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Abstract 

Excess  nitrogen  input  is  deleterious  to  coastal  waters,  resulting  in  deterioration  of 
the  water  quality,  increases  in  harmful  algal  blooms  and  disease  in  commercial  fish 
stocks.  A  significant  portion  of  this  nitrogen  enters  coastal  waters  through  groundwater 
systems.  Here  we  use  isotopic  and  molecular  biological  methods  to  identify  the 
populations  of  nitrifiers  and  denitrifiers,  where  they  occur,  and  what  levels  of  activity  are 
present  through  the  upper  four  meters  of  a  coastal  groundwater  system.  This  work  shows 
two  different  populations  of  putative  ammonia-oxidizing  archaea  (AOA)  based  on  the 
ammonia  monooxygenase  gene  (amoA),  one  shallow  population  most  closely  related  to 
open  ocean  water  column-like  sequences  and  a  deeper  population  that  is  more  closely 
related  to  estuarine-like  AOA.  Interestingly,  while  the  surface  population  has  a  potential 
nitrification  rates  (456  pmol  g'^  sediment  day'^)  similar  to  marine  sediments,  the  deeper 
population  does  not  show  detectable  evidence  of  nitrification.  Between  these  two 
archaeal  populations  resides  an  active  population  of  ammonia-oxidizing  bacteria  with 
similar  nitrification  rates  as  the  surface  AOA  population.  The  upper  meter  of  the  aquifer 
is  also  an  active  area  of  denitrification  as  evidenced  by  the  coincident  drop  in  nitrate 
concentration  and  increase  in  both  (up  to  -i-  20.1%o)  and  8^^0  (up  to  -i-  1 1.7%o), 
characteristic  of  groundwater  affected  by  denitrification.  16S  rRNA  gene  surveys  of  the 
organisms  present  in  the  upper  meter  also  are  similar  to  soil/sediment  type  environments 
including  many  potential  denitrifiers.  However,  nitrite  reductase,  nirS  and  nirK,  genes 
were  also  recovered  from  the  sediments  with  nirK  dominating  in  the  surface  sediments. 
This  contrasts  with  the  deep  salt  wedge,  where  the  microbial  community  16S  rRNA 
genes  appear  more  closely  related  to  marine  or  reducing  sediment/wastewater  type 
organisms,  and  nirS  genes  become  the  dominant  denitrification  gene. 
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Introduction 


Anthropogenic  nitrogen,  largely  in  the  form  of  nitrate  (NO3'),  is  a  deleterious 
pollutant  in  coastal  ecosystems.  Excess  nitrogen  addition  to  coastal  waters,  known  as 
eutrophication,  has  broad  deleterious  effects  on  coastal  ecosystems  including  increased 
hypoxia  and  anoxia,  altered  food  webs  and  degraded  sea  grass  beds,  kelp  beds  and  coral 
reefs.  These  impacts  result  in  decreased  biodiversity,  increased  frequency  and  duration  of 
harmful  algal  blooms,  and  compromised  finfish  and  shellfish  populations  (Valiela  et  al. 
1990;  Valiela  et  al.  1992;  Howarth  et  al.  2002).  Three  decades  after  the  Clean  Water  Act 
(1972),  the  nitrogen  flux  into  coastal  waters  has  increased  four  to  five-fold  and  as  much 
as  ten-fold  in  some  areas  as  a  result  of  human  activity  (Howarth  et  al.  2002).  Despite 
efforts  to  lower  aquatic  nitrogen  (N)  loading,  two  thirds  of  coastal  waterways  are 
moderately  to  severely  impacted  by  nutrient  pollution  (Bricker  et  al.  1999;  EPA  2001). 
With  rapid  population  growth  projected,  mainly  concentrated  within  50  miles  of  the 
coastline,  we  can  only  expect  N-loading  in  coastal  systems  to  increase.  To  better 
understand  the  processes  that  may  affect  the  eutrophication  of  coastal  waters  we  must 
understand  both  the  route  of  input  and  the  processes  that  may  remove  nitrogen  in  transit 
to  the  coastal  ocean.  The  work  presented  in  this  dissertation  examines,  in  a  coastal 
groundwater  aquifer,  (1)  the  balance  of  nitrification  and  denitrification,  (2)  which 
organism  are  involved  in  N-cycling  at  various  depths  and  (3)  what  are  the  potential  rates 
of  nitrogen  transformation  by  these  communities. 

It  has  long  been  known  that  much  of  the  terrestrially  derived  nitrogen  is 
transported  to  the  coastal  ocean  through  runoff  and  riverine  inputs.  However,  submarine 
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groundwater  discharge  (SGD),  defined  as  any  and  all  water  flowing  from  the  seabed  to 
the  continental  margin  regardless  of  composition,  has  been  recognized  as  a  mechanism  to 
transport  freshwater  into  coastal  waters  (Figure  1)  (Moore  1999;  Burnett  et  al.  2003; 
Moore  2006).  While  rivers  amount  to  a  point  source  of  constituents  to  the  coastal  ocean, 
and  thus  are  more  amenable  to  investigation,  SGD  could  be  viewed  as  a  spatially  broad 
signal  entering  the  coastal  ocean.  Due  to  the  large  spatial  and  temporal  variability  as  well 
the  remote  and  difficult  to  sample  nature  of  the  groundwater/sediment  matrix,  the  role  of 
SGD  in  coastal  biogeochemistry  has  been  understudied  and  underappreciated.  Despite  the 
inherent  difficulties  three  general  methods  have  been  used  to  quantify  the  contribution  of 
SGD  to  both  the  water  and  solutes  to  coastal  waters.  The  first  method  of  study  entails  the 
physical  measurements  of  discharge  rates  using  seepage  meters  or  flow  meters,  resulting 
in  a  calculation  of  the  local  water  flux  from  SGD  (Bugna  et  al.  1996;  Michael  et  al.  2003; 
ShoUcovitz  et  al.  2003;  Swarzenski  et  al.  2004).  As  a  second  approach,  discharge  rates 
can  be  estimated  using  tracer  techniques  usually  involving  the  tracking  of  radium  or 
radon  isotopes  (Moore  1996;  Charette  et  al.  2001;  Moore  et  al.  2008)  or  artificially 
injecting  compounds  such  as  sulfur  hexafluoride  (Cable  and  Martin  2008)  or  bromide 
(Hall  et  al.  1991)  into  the  aquifer  and  quantifying  the  flux  rate  of  the  tracer  to  the  coastal 
ocean.  Finally,  discharge  rates  can  be  estimated  using  various  modeling  approaches  that 
are  usually  based  on  the  mass  balance  of  water  and  the  assumption  of  steady  state  have 
been  used  to  estimate  groundwater  flux  into  coastal  waters  (Thompson  et  al.  2007;  Spiteri 
et  al.  2008).  Estimates  of  SGD  suggest  that  it  can  be  a  significant  contribution  of  water 
and  chemical  constituents  into  the  coastal  ocean,  with  water  fluxes  (freshwater  and 
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saltwater  recirculation)  as  high  as  40%  of  riverine  and  runoff  inputs  (Moore  1996). 
Freshwater  discharge  constitutes  6-10%  (Burnett  et  al.  2003)  of  the  total  discharge, 
roughly  3000-4000  km^  yr'^  (Zektser  and  Loaiciga  1993). 

The  water  flux  of  SGD  into  coastal  water  does  not  portray  the  full  impact  this 
process.  While  the  global  net  flux  of  water  from  SGD  is  relatively  small  when  compared 
to  river  inputs  on  the  estuary  scale  the  contribution  of  groundwater  to  the  water  flux  can 
vary  widely.  However,  water  flux  alone  is  not  the  only  impact  on  the  receiving  waters, 
the  composition  of  the  water  must  also  be  considered.  The  solute  load  of  SGD  can  differ 
greatly  from  the  solute  load  delivered  by  river  transport.  It  is  important  to  realize  that 
SGD  is  not  just  average  river  water  percolating  through  an  inert  system  of  permeable 
rock  and  debris  and  finally  empting  into  the  coastal  ocean  as  an  underground  river.  The 
transit  through  the  permeable  rock,  and  the  chemical  and  biological  reactions  occurring 
within  the  strata,  alter  the  composition  of  the  effluent  water.  Studies  examining  the 
leaching  of  elements  and  compounds  from  the  rock  and  the  ultimate  transport  to  the 
coastal  ocean  are  rare,  and  the  data  available  is  skewed  because  of  the  homogeneity  of  the 
rock  type  studied  (mostly  limestone  or  other  Karst  bedrocks  as  opposed  to  shale  or  other 
native  rocks),  however,  the  results  are  enlightening.  Calcium,  for  example,  is  a  major 
constituent  of  seawater  present  at  a  steady  state  concentration  of  10  mM.  In  order  to 
maintain  this  steady  state  concentration  given  the  calculated  inputs  of  calcium  to  the 
oceans,  groundwater  supply  of  calcium  to  the  ocean  must  be  20%  larger  than  the  largest 
known  source  of  the  calcium  to  the  ocean,  river  inputs  (Milliman  1993). 
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Calcium  is  not  the  only  element  or  compound  that  may  be  supplied  to  or  removed 
from  saline  groundwater  through  SGD  and  this  exchange  is  determined  by  the  rock 
involved  and  the  chemical  environment  at  the  location  of  groundwater/seawater  mixing 
analogous  to  the  particle  transport  and  flocculation  in  estuaries.  Moore  (1999)  called  this 
area  of  mixing  fresh  groundwater  and  seawater  the  subterranean  estuary  (STE).  Like 
traditional  estuaries,  STBs  are  areas  where  changes  in  the  chemical  nature  of  the  fluid 
results  in  flocculation,  sedimentation  and  ion  exchange  reactions.  However,  STBs  may 
also  have  very  different  chemical  characteristics  than  traditional  estuaries  because  of  a 
lack  of  communication  with  the  atmosphere  and  reaction  with  the  sediment  matrix  over 
the  relatively  long  transit  times  of  the  flowing  water.  The  chemical  reactions  modify  the 
rock  and  fluid  composition,  causing  certain  major  ions  to  be  lost  due  to  the  precipitation 
of  minerals  while  other  trace  elements  (Be,  Mn),  nutrients  (N,  P),  dissolved  organic  or 
inorganic  carbon  and  methane  may  be  released  to  the  coastal  ocean.  The  resulting  strong 
and  dynamic  chemical  gradients  in  a  subterranean  estuary  may  also  provide  overlapping 
habitats  for  microbial  populations  that  are  traditionally  thought  to  be  separate  due  to  the 
specific  environmental  controls  on  their  metabolism.  These  changes  in  both  chemical 
composition  and  biological  community  interactions  may  have  important  effects  on  the 
nutrient  load  of  advecting  groundwater  prior  to  discharge  in  to  the  coastal  ocean. 

Nutrients  delivered  to  the  coastal  ocean  by  groundwater  can  alter  the  ecology  of 
the  coastal  waters  as  described  above  in  the  balance  between  healthy  primary  production 
and  eutrophication.  The  transport  or  mobilization  of  nutrients,  trace  metals,  and  other 
compounds  may  be  due  to  the  natural  pressure  driven  flow  (including  tidal  flushing)  of 
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the  groundwater  aquifer  or  may  be  enhanced  due  to  anthropogenic  processes  such  as 
reduction  of  the  hydraulic  head,  thus  enhancing  salt  water  intrusion.  Many  reactions  can 
occur  during  the  intrusion  and  mixing  of  salt  water  into  the  groundwater  aquifer  due  to 
the  changes  in  the  ionic  strength  of  the  solution.  Desorption  of  some  bound  compounds 
from  the  substrate  or  flocculation  and  precipitation  of  other  compounds  is  expected  to 
occur  as  it  does  in  the  mixing  zones  of  surface  estuaries.  Radium  (Ra)  is  among  the 
compounds  that  desorb  during  the  mixing  process.  Seawater  naturally  contains  ^^^Ra  at 
an  activity  of  8  dpm/100  L  (Moore  1996),  any  excess  of  ^^^Ra  in  coastal  waters  must  be 
derived  from  other  sources.  Studies  have  shown,  by  eliminating  other  possible  sources, 
that  the  excess  ^^^Ra  is  an  indicator  of  groundwater  flow  into  the  coastal  ocean  (Krest  et 
al.  2000).  Furthermore,  the  activity  of  groundwater  ^^^Ra  roughly  correlates  with  salinity 
(more  mixing  desorbs  more  Ra)  (Moore  1996).  It  is  assumed  that  this  desorption  of  Ra 
coincides  with  desorption  of  other  compounds  including  dissolved  inorganic  nitrogen 
(DIN).  Unfortunately,  the  export  of  DIN  to  the  coastal  ocean  does  not  systematically 
follow  ^^^Ra  enrichments  due  to  other  sources  and  sinks  of  nitrogen  in  coastal  systems, 
including  nitrogen  fixation  and  nitrogen  uptake  by  vascular  plants.  Considering  these 
difficulties,  Krest  et  al.  (2000)  estimate  that  half  of  the  DIN  in  coastal  waters 
(extrapolated  from  a  marsh/estuary  system)  is  supplied  by  groundwater  flux.  However, 
physical  and  chemical  processes  are  not  the  only  factor  in  controlling  the  speciation  and 
amount  of  DIN  flux  to  the  coastal  ocean.  Microbiological  processes  can  also  influence 
DIN  export. 
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Nitrification,  Denitrification  and  Anaerobic  Ammonia  Oxidation: 

The  microbial  transformation  of  fixed  nitrogen,  mediated  by  members  of  the 
Bacteria  and  Archaea,  involves  the  processes  of  nitrogen  fixation,  nitrification, 
dissimilatory  nitrate  reduction  to  ammonia  (DNRA),  denitrification,  and  anaerobic 
ammonia  oxidation,  with  the  latter  two  processes  returning  bioavailable  nitrogen  back  to 
the  atmospheric  N2  reservoir  (Figure  2).  These  are  the  major  microbial  pathways  of 
inorganic  nitrogen  transformation  in  the  environment. 

Nitrification  is  the  microbial  transformation  of  ammonia  (NH3)  to  nitrate  (NO3') 
through  the  intermediates  hydroxylamine  (NH2OH)  and  nitrite  (NO2  ).  It  constitutes  an 
important  link  between  the  most  reduced  (NH3)  and  most  oxidized  (NO3')  forms  of  N  and 
generates  substrates  (NO3'  and  NO2  )  for  denitrification  and  anammox.  Nitrification  is  a 
multistep  process  which  involves  (at  least)  two  distinct  metabolic  groups  of  bacteria: 
ammonia-oxidizing  bacteria  (AOB)  and  nitrite-oxidizing  bacteria  (NOB),  though  archaea 
have  recently  been  implicated  in  ammonia  oxidation  as  well  (see  below).  AOB 
chemoautotrophically  oxidize  NH3  to  NO2  through  the  intermediate  NH2OH.  This 
pathway  requires  two  enzymes,  ammonia  monooxygenase  (AMO)  and  hydroxylamine 
oxidoreductase  (HAO),  to  conserve  energy  for  CO2  fixation  from  the  production  of 
nitrite.  With  the  exception  of  Nitrosococcus  sp.,  which  are  members  of  the  gamma- 
Proteobacteria,  the  majority  of  the  known  AOB  are  beta-Proteobacteria.  Thus,  until 
recently,  this  geochemically  important  group  of  microbes  was  believed  to  belong  to  a 
narrow  phylogenetic  group  (Rotthauwe  et  al.  1997).  Still,  with  the  low  abundance  of 
ammonia-oxidizing  bacteria  in  most  environments,  molecular  detection  and  diversity 
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studies  have  largely  relied  on  characterization  of  the  distribution  and  diversity  of  the  gene 
for  the  alpha  subunit  of  the  ammonia  monooxygenase  (amoA)  (Rotthauwe  et  al.  1997; 
Kowalchuk  et  al.  1999;  Purkhold  et  al.  2000;  Beman  and  Francis  2006). 

Estuaries  present  an  interesting  case  study  because  the  mixing  of  seawater  and 
freshwater  creates  significant  gradients  and  unique  chemical  environments.  Studies  in  the 
Chesapeake  Bay  and  Monterey  Bay  have  shown  that  the  diversity  of  amoA  genes  appears 
to  vary  with  salinity  along  the  estuaries  with  Nitrosospim-like  species  restricted  to  the 
mid-  to  high-salinity  regions  while  Nitrosomonas-Yike,  species  showing  no  trend 
(O’Mullan  and  Ward  2002;  Francis  et  al.  2003).  This  dependence  of  ammonia-oxidizer 
community  structure  on  salinity  has  implications  for  our  site,  where  strong  salinity 
gradients  co-occur  with  strong  gradients  in  nitrate  and  ammonium  concentrations. 

Much  is  known  about  the  physiology  and  biochemistry  of  bacterial  ammonia- 
oxidizers,  such  as  their  growth  substrates,  sensitivity  to  light,  and  ability  to  grow  at  low 
oxygen  tensions.  Indeed,  cultured  AOB  are  capable  of  growth  on  dissolved  oxygen  levels 
as  low  as  6  pM  (Carlucci  and  McNally  1969;  Goreau  et  al.  1980;  Lipschultz  et  al.  1981). 
These  fundamental  characteristics  of  AOB  have  shaped  our  view  of  when  and  where 
nitrification  may  be  an  important  player  in  nitrogen  transformations.  However,  AOB 
often  account  for  only  0.1%  of  the  community  in  pelagic  and  benthic  oceanic 
environments  (Ward  2000).  The  abundance  of  AOB  is  thus  not  representative  of  the 
importance  of  the  nitrogen  transformation  that  they  mediate  which  ultimately  leads  to 
oxidized  nitrogen  availability  for  removal  processes  such  as  denitrification. 
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Denitrification  is  the  stepwise  reduction  of  nitrate  to  dinitrogen  through  the 
intermediates  nitrite,  nitric  oxide  and  nitrous  oxide.  Because  of  the  requirement  for  a 
supply  of  oxidized  nitrogen  compounds,  denitrification  is  often  found  coupled  to 
nitrification  across  sediment  oxyclines.  Denitrifiers  are  facultative  anaerobic  heterotrophs 
and  chemoautotrophs  found  in  all  three  domains  of  life.  The  diverse  nature  of  denitrifying 
microbes  and  their  close  phylogenetic  affiliation  to  non-denitrifying  bacteria  and  archaea 
require  the  use  of  functional  gene  targets  to  analyze  the  diversity  of  denitrifying 
microbes.  Unlike  the  amoA  gene,  the  genes  that  encode  the  proteins  required  for 
denitrification  are  taxonomically  widespread  and  do  not  corroborate  the  phylogeny  built 
by  16S  rRNA  genes,  implying  horizontal  gene  transfer  (Jose  and  David  1999;  Song  and 
Ward  2003).  The  enzymes  of  interest  for  molecular  surveys  have  traditionally  been  nitrite 
reductase  (nirS,  nirK  genes)  (Braker  et  al.  2000;  Braker  et  al.  2001)  nitric  oxide  reductase 
{norB  gene)  (Braker  and  Tiedje  2003;  Casciotti  and  Ward  2005)  and  nitrous  oxide 
reductase  {nosZ  gene)  (Scala  and  Kerkhof  1998;  Scala  and  Kerkhof  1999).  Interestingly, 
the  diversity  of  the  denitrifying  population  has  been  shown  to  relate  to  the  carbon  content 
(quantity  and  age)  of  the  sediment,  rather  than  geochemical  gradients,  with  greater 
diversity  occurring  in  coastal  sediment  where  carbon  loads  are  both  greater  and  younger 
(Braker  et  al.  2000;  Braker  et  al.  2001). 

Anammox,  or  anaerobic  ammonia  oxidation,  involves  the  oxidation  of  ammonium 
with  nitrite  with  the  electron  acceptor  and  dinitrogen  as  the  end  product  (Van  de  Graaf  et 
al.  1995;  Jetten  et  al.  2001;  Dalsgaard  et  al.  2003;  Kuypers  et  al.  2003).  Anammox  is 
considered  a  shunt  in  the  nitrogen  cycle  because  the  pathway  is  capable  of  removing 
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reduced  forms  of  nitrogen  without  prior  oxidation  to  nitrate.  This  metabolism  was 
predicted  long  before  organisms  capable  of  it  were  discovered  (Broda  1977).  Since  then, 
a  representative  member  of  the  deeply  branching  group  within  the  Planctomycetales  has 
been  enriched  from  wastewater  treatment  facilities  and  sequenced  as  part  of  a  community 
genome  (Strous  et  al.  2006).  The  anammox  process  has  also  recently  been  shown  to  be 
important  for  nitrogen  removal  from  wastewater  treatment  systems  and  in  suboxic  marine 
sediments  and  water  columns  (Dalsgaard  and  Thamdrup  2002;  Dalsgaard  et  al.  2003; 
Kuypers  et  al.  2003;  Dalsgaard  et  al.  2005;  Schmid  et  al.  2007). 

Why  Archaea  may  be  important 

Archaea  have  been  traditionally  thought  of  as  the  prokaryotes  inhabiting  extreme 
environments  such  as  hot  springs,  hypersaline  lakes,  etc.  (Woese  et  al.  1978;  Barns  et  al. 
1994).  However,  Archaea  have  been  show  to  be  abundant  in  many  non-extreme 
environments  including  the  coastal  waters  (DeLong  1992)  and  open  ocean  waters 
(Fuhrman  1992),  Antarctic  water  (DeLong  et  al.  1994;  Murray  et  al.  1999),  terrestrial 
soils  (Buckley  et  al.  1998),  freshwater  sediments  (MacGregor  et  al.  1997)  and  even 
within  sponge  tissues  in  symbiosis  (Preston  et  al.  1996).  In  fact,  Archaea  may  account  for 
20-40%  of  the  marine  picoplankton  (DeLong  et  al.  1994;  Kamer  et  al.  2001).  The  non- 
thermophilic  crenarchaeota  are  the  most  widely  distributed  and  abundant  archaea  on  the 
planet  (DeLong  1998). 

The  distribution  of  archaea  in  the  marine  environment  has  been  a  focus  of  the 
early  inquiries  into  this  third  domain  of  life.  Open  ocean  16S  rRNA  gene  surveys  have 
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shown  that  below  the  euphotic  zone  (>150  m)  the  abundance  of  crenarchaeota  constitute 
the  major  fraction  of  the  picoplankton  and  by  1000  m  reach  number  equal  to  their 
bacterial  counterparts  (Massana  et  al.  1997;  Kamer  et  al.  2001).  Due  to  hybridization  of 
or  cellular  16S  rRNA  it  has  been  estimated  that  60-80%  of  the  microorganisms  present  in 
the  water  column  are  metabolically  active  (Karner  et  al.  2001).  However,  without 
cultured  representatives  of  non-extreme  organisms  it  was  difficult  to  ascertain  the 
ecological  role  of  the  archaeal  biomass  in  the  environment.  Interestingly,  positive 
correlations  between  crenarchaeal  distribution  and  nitrite  were  observed  in  the  Arabian 
Sea  (Damste  et  al.  2002),  coastal  ocean  (Murray  et  al.  1999),  and  in  low  oxygen 
conditions  in  the  Black  Sea  (Coolen  et  al.  2007).  Enrichment  culture  obtained  from  the 
North  Sea  showed  increases  in  crenarchaeota  and  amoA  gene  abundances  coinciding  with 
the  loss  of  ammonium  and  generation  of  nitrite  (Wuchter  et  al.  2006). 

More  recently,  genomic  surveys  have  shed  light  on  to  the  potential  ecology  of 
marine  crenarchaeota.  Shotgun  libraries  from  the  Sargasso  Sea  discovered  the  presence  of 
amoA  on  a  scaffold  containing  archaeal-linked  genes  (Venter  et  al.  2004).  This  finding 
suggested  that  some  pelagic,  marine  archaea  might  be  capable  of  nitrification. 
Subsequently,  Sargasso  Sea-like  amoA  genes  (archaeal  amoA  genes)  were  later  found  to 
occur  with  16S  rRNA  genes  from  crenarchaeota  in  a  large  genomic  fragment  isolated 
from  soils  (Schleper  et  al.  2005;  Leininger  et  al.  2006).  Archaeal  amoA  genes  have  also 
been  found  in  marine  waters  and  sediments,  including  Elkhorn  Slough,  CA  estuary 
sediments  (Erancis  et  al.  2005).  Additional  work  has  resulted  in  the  cultivation  of  the  first 
ammonia  oxidizing  archaea  (AO A)  related  to  the  abundant  Marine  Group  I 
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Crenarchaeota  (MGI)(Konneke  et  al.  2005),  demonstrating  unambiguously  that  a  marine 
archaeon  is  able  to  conserve  energy  chemolithoautotrophically  by  oxidizing  ammonia. 
Finally,  Wuchter  et  al.  (2006)  conclusively  showed  that  the  loss  of  nitrate  from  the  water 
column  of  the  North  Sea  is  correlated  with  archaeal  cell  and  archaeal  amoA  abundance 
while  bacterial  amoA  genes  occurred  at  much  lower  levels. 

Contrary  to  their  environmental  abundance  their  ecology  and  activities  are  less 
well  understood.  For  example,  it  is  known  that  AOB  are  inhibited  by  high  light  but  it  is 
unknown  if  AO  A  have  the  same  sensitivity.  There  is  also  some  evidence  that  AO  A  may 
not  be  strict  autotrophs  as  genomic  evidence  suggests  they  may  be  able  to  metabolize 
organic  carbon  as  well  as  carbon  dioxide  in  a  mixotrophic  lifestyle.  What  is  known  about 
their  physiologic  capacity  comes  from  only  a  handful  of  laboratory  cultures  and 
observations  of  their  distributions  in  a  range  of  environments.  Research  has  shown  AOA 
can  range  in  temperature  from  mesophilic  (Konneke  et  al.  2005)  to  thermophilic  (de  la 
Torre  et  al.  2008).  Their  ecological  significance  has  been  the  focus  of  studies  in  recent 
years,  almost  always  in  comparison  to  their  bacterial  counter  parts.  Relative  changes  in 
abundance  between  the  two  groups  has  been  correlated  with  salinity  and  oxygen  (Beman 
and  Francis  2006;  Mosier  and  Francis  2008;  Santoro  et  al.  2008),  C:N  (Mosier  and 
Francis  2008),  as  well  as  trace  metal  availability  (lead  and  nickel)  (Caffrey  et  al.  2007; 
Mosier  and  Francis  2008). However,  to  date  there  is  not  one  master  variable  that  seems  to 
determine  AOA  abundance. 


24 


Research  Approach 


Subterranean  estuaries  are  the  last  line  of  defense  against  nitrogen-polluted 
groundwater  seeping  into  the  coastal  ocean  through  SGD,  and  microbial  nitrogen  removal 
may  provide  an  important  environmental  service  in  these  settings.  Unfortunately,  little  is 
known  about  the  microbial  community  responsible  for  nitrogen  transformations  in 
subterranean  estuaries.  The  primary  objectives  of  this  study  are  to  determine  the 
diversity,  abundance  and  distribution  of  bacteria  and  archaea  involved  in  nitrification  and 
denitrification  in  the  Waquoit  Bay  subterranean  groundwater  estuary  and  to  determine  the 
potential  environmental  rates  of  these  processes.  Insight  into  these  processes,  the 
distribution  and  activity  of  the  microorganisms  involved,  will  allow  a  better  assessment 
of  nitrogen  loading  in  coastal  waters  and  the  natural  processes  that  may  attenuate  this 
load. 

To  begin  this  effort  the  microbial  communities  responsible  for  nitrification  and 
their  distribution  in  the  groundwater  system  must  be  described.  Chapter  2  surveys  the 
abundance  and  distribution  of  bacterial  and  archaeal  16S  rRNA  and  amoA  genes 
throughout  a  sediment  core  taken  at  Waquoit  Bay  in  June  2008.  The  phylogenetic 
distributions  of  the  archaeal  amoA  genes  are  also  presented.  The  results  of  this  work 
show  that  archaeal  16S  rRNA  gene  abundance  accounts  for  more  than  70%  of  the  cells 
deep  in  the  sediment  column,  as  the  pore  water  transition  from  fresh  groundwater  to  the 
deep  salt  wedge.  Looking  at  Marine  Group  I  Crenarchaeota,  a  subset  of  the  total  archaeal 
population  that  includes  the  AOA,  a  slightly  different  pattern  emerges.  MGI 16S  rRNA 
genes  show  two  peaks  in  abundance,  one  in  the  upper  0.5  meter  and  again  at  3  meters. 
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just  above  the  deep  salt  wedge.  The  distribution  of  archaeal  amoA  genes  is  consistent 
with  that  of  MGI  16S  rRNA  genes  with  a  roughly  1:1  ratio.  Bacterial  amoA  gene 
abundance  peaked  in  the  core  of  the  freshwater  plume  where  it  accounted  for  the  majority 
of  the  total  amoA  genes.  This  chapter  also  includes  the  first  experimental  comparison  of 
archaeal  amoA  primer  sets,  an  important  contribution  to  the  field  in  that  it  provides  a 
basis  for  the  interpretation  of  published  data  sets  with  differing  primer  sets  and  reaction 
conditions. 

Having  addressed  who  and  where  the  nitrifying  communities  are  in  the  aquifer, 
Chapter  3  focuses  on  the  role  of  denitrification  in  the  system.  To  do  this  we  describe  the 
chemistry  of  the  pore  waters  along  a  shore-normal  transect  as  well  as  the  abundance  of 
nitrite  reductase  genes  (nirS  and  nirK)  down  core.  Included  among  the  transect  data  are 
natural  abundance  isotopic  measurement  of  8^^Nno3,  8^^0no3  and  8^^0h2o-  The  data 
shows  a  loss  of  nitrate  occurs  in  the  upper  meter  of  the  sediment  column.  While  nitrate 
decreases  the  8^^Nno3  and  8^^0no3  increase  in  a  ratio  indicative  of  denitrification  in  a 
sandy  aquifer.  Unfortunately,  due  to  the  steep  nitrate  gradient  transitioning  from  the 
freshwater  plume  into  the  deep  salt  wedge  natural  abundance  isotopic  measurements  are 
not  available  over  this  transition.  NirK  is  the  dominant  denitrification  gene  in  the  surface 
meter  while  nirS  increases  rapidly  entering  into  the  deep  salt  wedge.  This  change  in  gene 
abundance  indicates  a  change  in  the  denitrification  community  from  the  surface  to  deep 
salt  wedge.  This  change  may  be  due  to  organic  carbon  or  some  other  environmental 
characteristic. 
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The  goal  of  Chapter  4  was  to  measure  potential  rates  of  nitrification  in  relation  to 
the  gene  profiles  reported  in  Chapter  2.  We  employ  the  use  of  tracers  to  track  the 
transfer  of  nitrogen  between  pools  and  calculate  the  rates  involved  based  on  the  transfer 
rates.  ^^N-labeled  ammonium  was  added  to  the  tracer  experiment  at  a  level  of  50  pM  and 
99.5-atom  %.  The  transfer  of  label  between  the  ammonium  and  nitrate  (no  nitrite  was 
evident)  pools  was  followed  along  with  transcription  of  specific  functional  genes  for 
nitrification.  We  show  that  the  surface  meter  is  a  zone  of  nitrification  mediated  by  AOA 
with  rates  significantly  faster  than  open  ocean  water  column  nitrification  but  slower  than 
nitrification  in  coastal  sediments.  In  the  core  of  the  freshwater  plume  (2  m)  we  find  AOB 
are  the  abundant  nitrifiers  with  potential  rates  equivalent  to  those  observed  in  the  surface 
sediments.  At  the  boundary  between  the  freshwater  and  deep  salt  wedge,  where  a  large 
population  of  AOA  had  been  described  in  Chapter  2,  we  find  no  chemical  evidence  for 
nitrification.  This  raises  some  interesting  questions  as  to  if  these  organisms  are  active  at 
this  depth  and  if  so  what  ecological  function  are  they  fulfilling? 

In  total,  the  data  presented  in  the  following  chapters  depicts  active  populations  of 
nitrifiers  and  denitrifiers  present  in  the  aquifer,  with  regions  of  potential  overlap  in  the 
upper  meter.  In  the  heart  of  the  large  nitrate  plume,  nitrification  appears  to  be  active 
while  denitrification  is  limited  in  this  portion  of  the  aquifer.  Of  course,  much  if  not  all  of 
this  fresh  groundwater  will  mix  with  seawater  prior  to  discharge.  If  the  communities  and 
activities  at  this  mixing  front  behave  similar  to  the  communities  in  the  upper  meter  of  the 
aquifer,  these  organisms  may  have  dramatic  effects  on  the  DIN  pools  discharged  into 
Waquoit  Bay. 
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A  B  C  D 


Figure  1.  Sampling  schematic  depicting  the  relative  position  of  the  piezometer  locations 
(A,  B,  C,  D)  in  relation  to  the  aquifer  structure.  The  aquifer  is  composed  of  three  distinct 
water  parcels,  (1)  a  surface  salt  zone  (stippled),  (2)  a  fresh  groundwater  plume  (white) 
and  (3)  a  deep  salt  wedge  (large  stipples).  The  mixing  zone  between  the  freshwater  and 
deep  salt  wedge  is  noted  with  a  gradient  in  stippling.  Pulse-coring  was  performed  at  site 
B  during  June  2008.  The  inset  shows  the  idealized  relative  concentration  variations  in 
nitrate,  oxygen,  salinity  and  ammonium  with  depth  in  the  aquifer  near  site  B. 
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Bacterial-Mediated  Nitrogen  Cycle 


Archaeal-Mediated  Nitrogen  Cycle 


Figure  2.  A  view  of  the  microbially-mediated  N-cycles.  Known  genes  for  each  pathway 
are  shown  in  blue,  anaerobic  N-evolving  pathways  are  shown  in  red.  Dashed  pathways  in 
the  archaeal-mediated  N-cycle  indicate  that  evidence  of  these  pathways  has  yet  to  be 
reported. 
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Abstract 


This  study  demonstrates  that  Archaea  are  abundant  members  of  the  microbial  community 
throughout  a  4  m  sediment  core  of  the  Waquoit  Bay  (Cape  Cod,  MA)  groundwater 
aquifer,  constituting  50%  or  more  of  the  total  16S  rRNA  gene  copies  at  each  depth 
examined  in  the  groundwater  aquifer.  16S  rRNA  gene  sequence  analyses  of  the  archaeal 
populations  showed  large  contributions  from  both  Euryarchaeota  (Methanosarcinales- 
like)  and  Crenarchaeota  (primarily  Cenarchaeales-like).  As  part  of  this  study,  the  result 
from  four  published  sets  of  archaeal  ammonia  monooxygenase  subunit  A  (amoA)- 
targeted  qPCR  primers  were  compared.  While  the  results  of  the  different  primer  sets  were 
generally  consistent,  the  Treusch  et  al.  primer  set  detected  fewer  amoA  copies  below  2  m 
depth  than  the  three  other  primer  sets.  Quantitative  PCR  results  for  both  16S  rRNA  and 
amoA  genes  revealed  that  ammonia-oxidizing  archaea  were  more  abundant  than  their 
bacterial  counterparts  at  most  depths  within  the  aquifer.  Ammonia-oxidizing  bacteria 
were  dominant  at  only  one  depth,  corresponding  to  the  core  of  the  freshwater  plume  and 
nitrate  peak,  where  they  accounted  for  approximately  80%  of  the  amoA  containing 
community.  From  archaeal  amoA  clone  libraries  we  found  that  the  ammonia-oxidizing 
population  from  the  surface  salt  zone  of  the  aquifer  was  characterized  by  amoA  gene 
sequences  most  closely  related  to  those  obtained  from  low-salinity  estuarine  sediments. 
On  the  other  hand,  the  archaeal  amoA  population  from  the  deep  salt  wedge  was 
dominated  by  sequences  most  closely  related  to  those  recovered  from  soils  and 
groundwater. 
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Introduction 


Delivery  of  excess  nitrogen  to  coastal  ecosystems  can  result  in  a  myriad  of 
deleterious  effects  ranging  from  excessive  or  harmful  algal  blooms  to  loss  of  sea  grass 
habitat,  hypoxia,  and  loss  of  fish  stocks  (Tomasky  et  al.  1999).  One  route  of  nitrogen 
input  into  coastal  water  is  through  submarine  groundwater  discharge,  which  can  also 
include  seawater  recirculated  through  coastal  sediments  (Moore  2006).  Mechanisms  of 
microbiological  transformation  of  the  nitrogen-containing  compounds  within  the 
groundwater  aquifer  prior  to  discharge  are  therefore  of  interest  because  they  can 
influence  both  the  quantity  and  composition  of  the  nitrogen  input  to  coastal  waters. 

Reduced  nitrogen  is  often  introduced  to  groundwater  as  ammonia  or  organic 
nitrogen  that  is  subsequently  catabolized  to  ammonia.  Nitrification  is  the  microbially- 
mediated  oxidation  of  reduced  nitrogen  (principally  ammonium)  to  nitrite  (NO2’)  and 
ultimately  to  nitrate  (NO3').  This  process  is  important  in  controlling  the  speciation  of 
nitrogen  between  the  reduced  and  the  oxidized  forms,  which  can  ultimately  be  removed 
from  the  environment  through  denitrification.  Nitrification  is  carried  out  by  at  least  three 
groups  of  microorganisms.  Ammonia  oxidation  to  nitrite  has  been  shown  to  be  the 
metabolic  function  of  specific  groups  of  ammonia-oxidizing  bacteria  (AOB)  (Bothe  et  al. 
2000)  and,  more  recently,  ammonia-oxidizing  archaea  (AO A)  (Konneke  et  al.  2005). 
Nitrite-oxidizing  bacteria  perform  the  latter  step  (NO2  to  NO3').  The  discovery  that 
certain  archaea  (specifically  Marine  Group  I  Crenarchaeota,  MGl)  that  are  ubiquitous 
and  abundant  in  the  environment  may  also  harbor  genes  responsible  for  ammonia 
oxidation  (Venter  et  al.  2004;  Treusch  et  al.  2005;  Hallam  et  al.  2006;  Leininger  et  al. 
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2006;  Wuchter  et  al.  2006;  Schauss  et  al.  2009)  has  prompted  reconsideration  of  the  roles 
of  bacteria  and  archaea  in  the  nitrification  process  and  where  (the  environmental 
contexts)  ammonia  oxidation  may  occur.  In  this  study,  we  apply  molecular  techniques  to 
ascertain  1)  who  the  most  abundant  ammonia-oxidizing  organisms  are  in  the  Waquoit 
Bay  subterranean  estuary  (Cape  Cod,  MA),  and  2)  where  in  the  gradients  of  nutrients, 
salinity,  and  oxygen  within  the  aquifer  the  communities  are  most  abundant  and  diverse. 

Materials  and  Methods 

Study  Area  and  Sample  Collection 

The  groundwater  aquifer  at  the  study  site  in  the  Waquoit  Bay  National  Estuarine 
Research  Reserve  (41°34’49”  N,  70°3r27”  W,  Waquoit,  MA  USA)  consists  of  outwash 
gravel,  sand  and  silt  deposited  during  the  retreat  of  the  Wisconsin  glaciation  (Oldale 
1976).  The  upper  10  m  of  sediment  are  homogeneous  and  highly  permeable,  providing 
conduits  for  transport  of  fresh  groundwater  toward  the  bay  and  tidally-forced  saltwater 
intrusions  (Charette  et  al.  2005).  The  Waquoit  Bay  groundwater  aquifer  has  been  the 
focus  of  several  studies  addressing  the  importance  of  the  groundwater  discharge  into  the 
bay  for  trace  metal,  nutrient,  and  pollutant  budgets  (Valiela  et  al.  1992;  Sham  et  al.  1995; 
Charette  et  al.  2001;  Testa  et  al.  2002;  Talbot  et  al.  2003),  and  a  long-term  monitoring 
program  for  nutrient  advection  into  the  bay  has  been  established  at  this  site  (Valiela  et  al. 
1990;  McClelland  and  Valiela  1998;  Michael  et  al.  2003). 
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Pore  water  and  Sediment  Sampling 

Pore  water  was  collected  in  June  2008  at  the  head  of  Waquoit  Bay  to  a  depth  of  4 
m  with  a  vertical  resolution  of  15-30  cm  using  a  piezometer  with  a  stainless  steel 
“Retract- a-tip”  (AMS,  Idaho).  After  flushing  with  at  least  one  liter  of  pore  water  (more 
than  6-fold  flushing  of  tubing  volume),  pH,  dissolved  oxygen,  salinity  and  temperature 
were  measured  using  a  YSI  600R  multi-probe  calibrated  following  the  manufacturer’s 
protocols  (Charette  and  Sholkovitz  2002).  These  indices  were  used  to  guide  pore  water 
sampling  and  sediment  coring  efforts. 

Triplicate  pore  water  samples  were  collected  and  filtered  through  a  0.2  pm  pore- 
size  Sterivex  filter  for  nutrient  concentration  (NO3',  NO2  ,  NH4'^)  analyses.  These  samples 
were  analyzed  within  12  hours  upon  return  to  the  laboratory.  Nitrite  concentrations  were 
measured  colorimetrically  using  sulfanilamide  and  N-(l-napthyl)  ethylenediamine  (Pai  et 
al.  1990).  Nitrate  plus  nitrite  (NOx)  concentrations  were  measured  by  chemiluminescence 
following  a  hot  vanadium  reduction  (Cox  1980;  Garside  1982;  Garside  1985). 
Ammonium  was  measured  using  the  phenolhypochlorite  method  (Solorzano  1969). 

After  pore  water  sampling,  a  sediment  core  was  collected  within  0.5  m  of  the 
piezometer  using  a  hand-operated  bailer-boring  auger  (Eijkelkamp,  The  Netherlands)  to  a 
depth  of  four  meters.  The  coring  effort  occurred  within  1  hr  of  low  tide.  The  recovered 
sediment  was  sampled  at  30-50  cm  resolution  using  aseptic  techniques,  and 
approximately  1-5  g  of  sediment  was  immediately  frozen  in  liquid  nitrogen  (on  site),  and 
then  stored  at  -80  °C  in  the  laboratory  for  later  DNA  extraction  and  analysis. 
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DNA  Extraction 


DNA  extractions  were  performed  in  triplicate  from  1  g  of  sediment  using 
lysozyme  (100  mg  mL'^)  and  proteinase  K  (20  mg  mL'^)  digestions  followed  by  physical 
disruption  by  successive  freeze-thaw  cycles  and  phenohchloroformdsoamyl  extraction 
(Edwards  et  al.  2000).  Extracted  DNA  was  precipitated  using  isopropyl  alcohol  and  3  M 
sodium  acetate  with  final  resuspension  in  50  pL  of  nuclease-free  water  or  Tris-EDTA 
buffer  (TE).  The  extracts  were  quantified  at  260  nm  using  a  Nanodrop  spectrophotometer 
(ThermoEisher  Scientific)  and  diluted  to  appropriate  concentrations  (see  below)  for 
amplification  and  quantification  by  polymerase  chain  reaction  (PCR)  methods  using 
nuclease-free  water. 

Quantitative  PCR 

Quantitative  PCR  (qPCR)  was  employed  to  determine  the  gene  copy  numbers  of 
bacterial  and  archaeal  16S  rRNA  and  ammonia  monooxygenase  subunit  A  (amoA)  genes 
at  30-50  cm  resolution  throughout  the  core  (Table  1).  Triplicate  DNA  extractions  of  each 
sample  were  diluted  to  10  ng  pl'^  and  amplified  using  either  the  iQ  SYBR  green 
Supermix  (Bio-Rad)  or  a  cocktail  consisting  of  1  X  colorless  master  mix  (Promega, 
Madison,  WI),  10  uM  dNTPs  (Promega),  0.5  X  SBYR-Green  I  (Invitrogen,  Carlsbad 
CA),  20nM  fluorescein  (Sigma),  1  mM  MgCl2  (Promega),  forward  and  reverse  primers 
(300-1000  uM;  Table  1)  and  1.25  units  of  GoTaq  Elexi  DNA  polymerase  (Promega). 
Pour  different  sets  of  qPCR  primers  (Prancis  et  al.  2005;  Treusch  et  al.  2005;  Wuchter  et 
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al.  2006;  Mincer  et  al.  2007)  were  used  to  analyze  the  abundance  of  archaeal  amoA 
(Table  1).  All  qPCR  amplifications  were  repeated  on  three  separate  days  with  optimized 
thermal  programs  using  an  iCycler  thermocycler  (Bio-Rad,  Hercules,  CA)  (Table  1). 
qPCR  efficiencies  ranged  from  88  -  105%  for  all  reactions  resulting  in  linear  standard 
curves  (r  =  0.90  -  0.99).  Melt  curves  were  performed  on  all  reactions  to  confirm 
specificity. 

Archaeal  amoA  Primer  Screening 

To  examine  the  target  specificity  of  the  four  different  archaeal  amoA  gene  primer 
sets  from  Francis  et  al.  (Francis  et  al.  2005),  Mincer  et  al.  (Mincer  et  al.  2007),  Wuchter 
et  al.  (Wuchter  et  al.  2006)  and  Treusch  et  al.  (Treusch  et  al.  2005),  the  primer  sequences 
were  aligned  to  the  amoA  gene  sequences  reported  from  the  Sargasso  Sea  WGS 
(AACY01435967),  Nitrosopumilus  maritimus  SCMl  (YP_00 1582834)  and  soil  fosmid 
54d9  (AJ627422)  using  the  Geneious  software  package  (Drummond  et  al.  2009).  These 
primer  sets  were  also  aligned  to  the  amoA  gene  sequences  recovered  from  the  Waquoit 
Bay  aquifer  sediments  to  determine  the  target  range  and  specificity  of  each  primer  set. 

PCR,  Cloning  and  Sequencing 

For  each  target  gene  (archaeal  16S  rRNA  and  archaeal  amoA),  three  replicate 
amplifications  from  0.5  m,  1  m,  2  m,  2.5  m,  3  m  and  4  m  were  conducted  according  to 
published  protocols  (Table  1  and  references  therein).  The  products  from  each  depth  were 
pooled  and  cleaned  with  the  SV  Wizard  PCR  Clean-up  (Promega),  ligated  following  the 
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manufacturer’s  protocols  (p-Gem  T-Easy  vector,  Promega),  and  transformed  into 
chemically  competent  E.  coli  JM109  cells.  Plasmid  purification  and  sequencing  were 
performed  at  the  Keck  Facility,  Josephine  Bay  Paul  Center  at  the  Marine  Biological 
Laboratories  (JBPC-MBL,  Woods  Hole,  MA).  Sequencing  chemistry  followed  a 
modified  Big-Dye  reaction  containing  0.5  pi  of  Big-Dye  Terminator,  0.4  pL  of  primer 
(M13R  5’-CAGGAAACAGCTATGAC,  15  pM),  0.1  pL  dimethyl  sulfoxide,  1/3  X 
reaction  buffer,  template  (200-400  ng  DNA)  and  water  to  a  final  volume  of  6  pL.  The 
recommended  thermal  protocol  of  60  cycles  of  96  °C  for  lOsec.,  50  °C  for  5sec  and  60  °C 
for  4  minutes  was  used  for  cycle  sequencing. 

Sequence  analysis 

Base  calling  and  sequence  quality  checks  were  performed  by  the  Phred  and  Phrap 
software  packages  (Ewing  and  Green  1998;  Ewing  et  al.  1998).  The  Lucy  software  (Chou 
and  Holmes  2001)  was  used  for  vector  trimming.  Processed  sequences  were  searched 
against  the  non-redundant  nucleotide  (nr/nt)  database  using  the  blastn  algorithm. 
Sequences  were  imported  into  ARB  (Ludwig  et  al.  2004),  aligned  with  the  ARB  aligner, 
manually  checked  and  corrected  then  added  by  parsimony  to  an  ARB  alignment  from 
Greengenes  (16S  rRNA  genes)  or  aligned  archaeal  amoA  gene  sequences  from  the 
GenBank  database  (Benson  et  al.  2000).  Conservative  column  filters  were  constructed 
masking  any  sites  missing  data  for  any  of  the  sequences.  Trees  were  constructed  on 
aligned  sequences  using  the  ARB  neighbor  joining  function. 
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Sequence  Accession  Numbers 

All  sequences  generated  and  used  in  this  study  have  been  reported  to  GenBank 
under  the  accession  numbers  HM 1 60 162-HM 160367. 

Results 

Environment  Context 

Long-term  chemical  monitoring  has  helped  to  elucidate  the  biogeochemical 
setting  within  the  Waquoit  Bay  subterranean  estuary  (Charette  et  al.  2001;  Talbot  et  al. 
2003;  Charette  et  al.  2005;  Charette  and  Sholkovitz  2006;  Mulligan  and  Charette  2006; 
Kroeger  and  Charette  2008).  The  most  prominent  features  observed  during  the  current 
study  include  (1)  a  surface  brackish  layer  (salinity  up  to  4)  with  low  dissolved  oxygen 
(<50  pM),  nitrate  (<50  pM)  and  ammonium  (<5  pM)  concentrations,  (2)  an  oxidized 
freshwater  parcel  with  high  levels  of  both  nitrate  (up  to  350  mM)  and  dissolved  oxygen 
(up  to  175  pM),  and  (3)  a  deep  reducing  salt  wedge  that  lacked  nitrate  but  contained  up  to 
53  pM  ammonium  (Figure  lA).  These  zones  provide  a  variety  of  niches  that  could  harbor 
distinct  microbial  communities  adapted  to  different  salinities,  oxygen  levels,  or  nutrient 
concentrations. 

16S  rRNA  gene  abundance 

The  abundance  of  bacterial  16S  rRNA  genes  ranged  from  3.9  x  10  to  1.3  x  10 
copies  per  gram  of  sediment  with  the  maximum  bacterial  16S  rRNA  gene  abundance 
observed  at  two  meters  depth  (Figure  IB).  Archaeal  16S  rRNA  gene  abundance  was 
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comparable  with  bacterial  16S  rRNA  genes  in  the  upper  1.7  meters  with  a  range  from  4.8 
X  10  to  7.3  X  10  copies  per  gram  of  sediment.  Below  1.7  meters,  the  archaeal  16S  rRNA 
gene  abundance  increased  gradually,  reaching  a  maximum  of  6.3  x  10^  copies  per  gram 
of  sediment  (13-fold  increase)  by  three  meters  depth. 

Taking  a  closer  look  at  the  archaeal  subgroup  MGl,  using  specific  qPCR  primers 
(Mincer  et  al.  2007),  we  found  a  two  order  of  magnitude  range  in  abundance  of  MGl  (4.5 
X  10"^  to  9.8  X  10^  copies  g'^  sediment )  over  the  studied  depth  interval.  Based  on  these 
qPCR  data,  MGl  apparently  accounted  for  less  than  one  percent  of  the  total  archaeal  16S 
rRNA  genes  at  all  but  the  shallowest  depths.  However,  this  minor  contribution  to  the 
community  abundance  was  not  observed  in  the  clone  libraries  (see  below),  perhaps 
indicating  an  overestimation  of  total  archaeal  abundance.  MGl  16S  rRNA  gene 
abundance  follows  archaeal  amoA  gene  abundance  in  a  nearly  1:1  ratio  except  for  the 
deepest  depth.  This  consistency  is  also  seen  in  other  studies  (for  example  Mincer  et  al. 
2007)  and  shows  our  estimation  of  these  genes  is  likely  robust.  Unlike  the  profile  of  total 
archaeal  abundance,  the  abundance  of  MGl  Crenarchaea  showed  two  peaks  in  the 
sediment  column  (Figure  IB).  The  first  peak  occurred  within  the  surface  meter,  reaching 
an  abundance  of  8.2  x  10^  copies  per  gram  of  sediment.  MGl  abundance  then  declined 
steadily  into  the  middle  portions  of  the  sediment  core,  reaching  a  minimum  of  4.5  x  lO"^ 
copies  per  gram  of  sediment  at  1.7  m.  The  second  MGl  peak,  at  a  depth  of  2.8  m, 
reached  an  abundance  of  9.8  x  10^  copies  per  gram  of  sediment.  Thus,  abundances  of 
MGl  rRNA  genes  increased  by  two  orders  of  magnitude  from  1.7  m  to  2.8  m,  whereas 
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the  abundance  of  total  archaeal  16S  rRNA  genes  increased  only  one  order  of  magnitude 
between  these  depths. 

Archaeal  amoA  gene  distributions 

The  results  from  all  four  qPCR  primer  sets  targeting  archaeal  amoA  genes  (Table 
1)  provided  reproducible  (self-consistent)  data  under  the  recommended  amplification 
conditions.  Moreover,  all  four  primer  sets  yield  the  same  range  of  abundances  for  the 
amoA  gene  throughout  the  upper  half  of  the  sediment  core  (Figure  1C).  However,  below 
2.8  m  the  Treusch  et  al.  (Treusch  et  al.  2005)  primer  set  yielded  archaeal  amoA  gene 
abundances  that  were  approximately  1  order  of  magnitude  lower  than  the  other  three 
primer  sets  (Figure  1C).  In  addition,  the  Francis  et  al.  (Francis  et  al.  2005)  primers  yield  a 
low  (order  of  magnitude)  estimate  compared  to  the  Mincer  et  al.  (Mincer  et  al.  2007)  and 
Wuchter  et  al.  (Wuchter  et  al.  2006)  primers  at  a  depth  of  2.7  m.  However,  given  the 
general  consistency  between  the  Wuchter  et  al.  (Wuchter  et  al.  2006),  Mincer  et  al. 
(Mincer  et  al.  2007),  and  Francis  et  al.  (Francis  et  al.  2005)  primer  sets  throughout  the 
sediment  core  (Figure  1C),  we  have  chosen  to  present  our  amoA  gene  data  as  obtained  by 
the  Francis  et  al.  (Francis  et  al.  2005)  primer  set  for  the  remainder  of  the  discussion. 

Archaeal  amoA  genes  occurred  in  high  abundance  in  the  surface  and  at  2.8  m  in 
the  core,  with  approximately  8.3  x  10  and  1.1  x  10  amoA  copies  per  gram  of  sediment 
respectively  (Figure  ID).  The  archaeal  amoA  gene  abundance  remained  low  until  2  m, 

n 

then  increased  to  a  second  maximum  of  1.1  x  10  copies  per  gram  sediment  at  2.8  m  and 
decreased  slightly  toward  the  bottom  of  the  sediment  core.  Each  primer  set  confirmed 
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these  main  features,  with  an  approximately  100-fold  decrease  in  amoA  gene  abundance 
between  the  surface  and  1.5  m  and  a  secondary  peak  at  2.8  m  (Figure  1C).  The  archaeal 
amoA  abundance  was  not  strongly  correlated  with  any  of  the  environmental  parameters 
measured.  The  best  observed  correlation  was  a  negative  correlation  between  archaeal 
amoA  gene  abundance  and  redox  potential  (Figure  2,  r^=  0.43). 

Normalizing  the  archaeal  amoA  gene  abundance  to  archaeal  16S  rRNA  gene 
abundance  revealed  that  the  AOA  component  accounted  for  a  maximum  of 
approximately  1%  of  the  total  archaeal  population  in  the  upper  0.5  m,  assuming  one  copy 
of  amoA  per  AOA  genome  (Mincer  et  al.  2007).  However,  comparison  of  archaeal  amoA 
gene  abundance  to  MGl  16S  rRNA  gene  abundance  showed  a  gene  ratio  close  to  1 
throughout  most  of  the  sediment  column,  abruptly  climbing  to  2  at  3.8  m.  Like  amoA,  a 
correlation  was  observed  between  MGl  16S  rRNA  gene  abundance  and  redox  potential 
(Figure  2,  r^  =  0.47). 

Bacterial  amoA  gene  distribution 

Bacterial  amoA  abundance  was  low  throughout  the  sediment  column,  with 
approximately  3.3  x  lO"^  copies  per  gram  of  sediment  at  the  surface  sediments,  increasing 
to  6.0  X  10^  copies  per  gram  of  sediment  at  2  m  (Figure  ID).  Below  2  m,  bacterial  amoA 
dropped  to  1.3  x  10  copies  per  gram  of  sediment.  Between  a  depth  of  1.2  and  2  m,  the 
bacterial  amoA  gene  reached  a  maximum  abundance  (6.0  x  10^  copies  per  gram  of 
sediment),  accounting  for  up  to  80%  of  the  total  amoA  gene  pool  at  2  m.  However,  except 
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for  this  interval,  the  bacterial  amoA  gene  was  1-2  orders  of  magnitude  less  abundant  than 
the  archaeal  version. 

Archaeal  16S  rRNA  gene  diversity 

The  archaeal  16S  rRNA  genes  recovered  from  the  aquifer  sediments  belonged  to 
both  the  Euryarchaeota  and  Crenarchaeota  subdivisions,  with  roughly  equal  contributions 
from  the  Euryarchaeota  (36  out  of  84,  43%),  and  the  Crenarchaeota  (48  out  of  84,  57%), 
for  which  MGl  Crenarchaeota  accounted  for  the  majority  of  the  sequences  (Eigure  4). 
Therefore,  in  contrast  to  the  qPCR  results,  which  showed  MGl  abundance  to  be 
approximately  1%  of  total  archaea,  the  archaeal  clone  libraries  were  approximately  50% 
MGl.  As  discussed  below,  this  discrepancy  may  be  the  result  of  biases  in  PCR 
amplification  and/or  cloning.  The  Crenarchaeota  sequences  were  recovered  from  all 
depths  in  the  aquifer,  and  they  were  most  closely  related  to  Cenarchaeales-Vike,  sequences 
(which  includes  Nitrosopumilus  maritimus),  although  sequences  were  also  recovered 
from  pSE  12-like  organisms  (4  out  of  48,  8%).  Both  MGl  and  PSE12  clades  have  been 
implicated  in  nitrification  (Konneke  et  al.  2005;  Hallam  et  al.  2006;  Mincer  et  al.  2007). 

A  third  group  of  organisms  related  to  SAGMA-like  (South  African  Gold  Mine  Archaea) 
subgroup  of  the  Crenarchaeota  were  also  recovered  throughout  the  aquifer  (7  of  48, 

15%).  It  is  unknown  whether  these  organism  harbor  amoA  genes.  Related  16S  rRNA 
genes  have  been  found  in  many  environments  including  surface  marine  waters,  Antarctic 
bathypelagic  and  other  deep  marine  sediments  (Gillan  and  Danis  2007),  uranium  mine 
tailings  (Takai  et  al.  2001),  and  agricultural  soils  (Mclnerney  1997). 
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Archaeal  amoA  Clone  Libraries 


The  archaeal  amoA  clone  libraries  (consisting  of  343  clones  total,  92  from  1  m,  96 
from  2  m,  86  from  3m  and  69  from  4  m)  from  the  aquifer  sediments  revealed  two  distinct 
communities  within  the  system.  The  first  community  was  mostly  found  in  the  upper  two 
meters  of  the  aquifer  sediments  and  a  second  community  between  two  and  four  meters 
within  the  sediments.  The  amoA  genes  recovered  from  the  upper  two  meters  clustered 
primarily  with  Water  Column-like  sequences  within  the  archaeal  amoA  phylogeny 
(Figure  5)  (Francis  et  al.  2005;  Beman  and  Francis  2006).  This  group  includes  sequences 
from  open  ocean,  coastal  waters  and  estuarine  systems  (Beman  and  Francis  2006;  Mosier 
and  Francis  2008).  The  sequences  from  Waquoit  Bay  were  most  similar  to  sequences 
from  the  low-salinity  sediments  of  North  San  Francisco  Bay  (21,  39).  Below  two  meters 
the  recovered  amoA  genes  revealed  a  different  community.  The  sequences  obtained  from 
this  deeper  community  were  distributed  within  a  cluster  of  archaeal  amoA  genes  obtained 
from  sediment,  soil  and  groundwater  (Hansel  et  al.  2008;  van  der  Widen  et  al.  2009). 

Discussion 

Many  studies  concerning  nitrogen  biogeochemistry  in  coastal  groundwater  focus 
on  the  role  of  denitrification  as  a  mechanism  of  nitrogen  loss  (LaMontagne  et  al.  2002; 
Addy  et  al.  2005)  with  less  consideration  for  how  this  process  fits  in  the  broader  context 
of  nitrogen  transformations  comprising  the  nitrogen  cycle.  For  example,  denitrification 
requires  a  source  of  nitrate,  which  can  be  transported  into  the  system  or  produced  locally. 
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Indeed,  a  great  deal  of  nitrogen  is  thought  to  be  lost  through  coupled  nitrification- 
denitrification  in  estuarine  systems.  A  few  previous  studies  have  used  modeled  and 
measured  chemical  fluxes  to  quantify  the  fluxes  of  nitrogen  within  the  Waquoit  Bay 
subterranean  estuary  (Kroeger  and  Charette  2008;  Spiteri  et  al.  2008).  Both  approaches 
have  suggested  that  nitrification  may  occur  in  both  at  the  surface  and  deep  salt  transition 
zone,  where  both  ammonia  and  oxygen  are  present.  Here  we  have  used  molecular 
techniques  to  describe  the  microbial  communities  harbored  by  the  aquifer  sediments  and 
examine  whether  these  communities  contain  the  genetic  potential  for  nitrification.  We 
also  report  the  congruity  of  the  major  published  primer  sets  for  the  archaeal  amoA  gene. 
Furthermore,  we  examine  the  relative  contribution  of  AOA  and  AOB  to  the  potential 
nitrifying  community,  focusing  on  the  potential  changes  in  their  abundance  and 
community  structure  in  the  different  zones  of  the  Waquoit  Bay  subterranean  estuary. 

16S  rRNA  gene  abundance 

The  observed  archaeal  16S  rRNA  gene  abundance  falls  within  the  range  of  values 
reported  for  near-surface  marine  sediments  (Schippers  and  Neretin  2006).  In  this 
subterranean  estuary,  the  increase  in  the  archaeal  16S  rRNA  gene  abundance  at  depth 
corresponds  to  a  region  that  is  dynamic  and  is  likely  exposed  to  both  saline  and  fresh  pore 
waters  over  the  course  of  the  annual  hydrologic  and  tidal  cycles.  As  the  intruding 
saltwater  wedge  advances  shoreward  during  the  drier  seasons  (lower  hydraulic  head),  this 
water  carries  with  it  a  greater  ionic  strength,  lower  oxygen  levels,  and  higher 
concentrations  of  ammonia.  Interestingly,  the  bacterial  16S  rRNA  gene  does  not  show 
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this  same  increase  at  depth.  Although  it  is  tempting  to  speculate  about  which 
environmental  factor(s)  may  contribute  to  the  observed  increase  in  archaea  at  depth,  16S 
rRNA  gene  abundance  does  not  show  a  strong  correlation  to  any  of  the  collected 
environmental  data  (NOa’,  salinity,  O2,  etc.). 

In  contrast  to  total  archaeal  16S  rRNA  abundance,  the  abundance  of  the  MGl  16S 
rRNA  genes  did  show  local  maxima  in  the  near  surface,  as  well  as  at  depth  in  the 
sediment  profile  (Figure  IB).  These  peaks  are  coincident  with  local  maxima  in  the 
abundance  of  archaeal  amoA  genes.  Moreover,  the  relative  abundance  of  MGl  16S  rRNA 
gene  copies  was  very  similar  to  the  archaeal  amoA  gene  copies  throughout  the  core, 
indicating  approximately  1  copy  of  amoA  per  16S  rRNA  gene  at  all  but  the  deepest  site, 
where  the  ratio  approached  2.  This  result  is  similar  to  that  obtained  by  Beman  et  al. 
(Beman  and  Francis  2006)  in  a  Mexican  estuary,  Wuchter  et  al.  (Wuchter  et  al.  2006)  in 
the  North  Sea,  and  by  Mincer  et  al.  (Mincer  et  al.  2007)  in  the  north  Pacific  subtropical 
gyre  (save  for  the  depth  interval  between  100  and  300  m  in  the  Pacific  data  set  where 
pSL  12-like  organisms  were  suspected  of  contributing  to  the  amoA  gene  pool).  An 
amoA:l6S  rRNA  gene  ratio  of  1  is  also  predicted  from  the  available  genomes  of 
ammonia-oxidizing  archaea  {Nitrosopumilus  maritimus  SCMl  and  Cenarchaeum 
symbiosum  A),  which  each  containing  one  copy  of  the  amoA  and  16S  rRNA  genes  per 
genome. 
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Archaeal  16S  rRNA  gene  libraries 

In  contrast  to  MGl  accounting  for  only  1%  of  the  total  archaeal  16S  rRNA  genes 
by  qPCR,  MGl  sequences  comprised  more  than  half  the  sampled  archaeal  16S  rRNA 
clone  library.  These  results  could  be  due  to  PCR  primer  bias,  cloning  or  library  sampling 
bias  and/or  changes  in  amplification  chemistries  and  efficiencies  between  the  two  PCR 
methods.  At  this  point,  we  do  not  fully  understand  the  reason(s)  for  this  discrepancy. 

The  presence  of  Euryarchaeota  in  this  environment  is  not  unexpected;  however, 
the  high  frequency  of  these  sequences  in  the  archaeal  16S  rRNA  clone  libraries  was 
surprising.  The  recovered  Euryarchaeota  sequences  all  fell  within  the  Methanosarcinales 
clade,  which  contains  anaerobic,  methanogenic  organisms.  This  is  an  interesting  finding 
given  the  bulk  of  the  sampled  water  is  oxic  and  should  inhibit  methanogenesis.  However, 
the  confluence  of  the  deep  salt  wedge  and  fresh  groundwater  plume  may  contain  low 
enough  levels  of  oxygen  and  to  favor  methanogenesis.  The  presence  of  these  autotrophs 
is  consistent  with  the  low  DOC  concentrations  typically  encountered  in  this  environment 
(Charette  et  al.  2005;  Charette  and  Sholkovitz  2006;  Bone  et  al.  2007).  However,  to  our 
knowledge,  neither  the  presence  of  methane  nor  the  activity  of  methanogens  have  been 
examined  in  this  aquifer. 

Archaeal  amoA  Gene  Libraries 

The  archaeal  amoA  gene  libraries  from  the  aquifer  sediments  revealed  two  distinct 
communities  within  the  system,  a  water  column  group  and  a  soils/sediments  group.  In  the 
Waquoit  Bay  groundwater  system,  the  water  column-like  sequences  were  recovered 
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exclusively  from  the  upper  two  meters  of  the  sediment  column.  Given  the  high 
abundance  and  subsurface  peak  of  archaeal  amoA  genes  in  this  upper  meter  it  is  unlikely 
that  these  sequences  originated  from  organisms  simply  trapped  in  the  sediment  column 
from  tidal  exchange  with  the  bay.  Moreover,  the  close  similarity  between  our  sequences 
and  those  recovered  from  low-salinity,  low  nutrient,  estuarine  sediments  (Francis  et  al. 
2005;  Mosier  and  Francis  2008)  suggests  the  possibility  of  a  niche  for  this  particular 
group  of  AO  A.  Below  two  meters,  the  phylogenetic  affiliation  of  the  recovered  amoA 
genes  shifted  from  water  column-like  to  sediment,  soil  and  groundwater-like  sequences. 
The  majority  of  the  sequences  recovered  from  the  deeper  section  were  most  closely 
related  to  sequences  recovered  from  groundwater  treatment  plants  (van  der  Widen  et  al. 
2009)  and  terrestrial  soils  (Hansel  et  al.  2008). 

Similar  “niche  separation”  has  been  reported  in  other  estuarine  environments, 
with  sequence  types  clustering  geographically  within  a  single  estuary  (Beman  and 
Francis  2006).  Unfortunately,  there  is  still  great  uncertainty  about  the  environmental 
conditions  or  physiological  properties  that  may  influence  the  distribution  of  different 
AOA  clades  (Erguder  et  al.  2009).  The  driving  force  behind  this  difference  in 
communities  is  not  evident  given  our  (or  other  published)  chemical  data  from  Waquoit 
Bay.  The  difference  in  community  composition  may  contribute  to  the  mismatch  in 
archaeal  amoA  abundance  determined  using  different  qPCR  primer  sets,  and  primer 
selection  should  be  considered  carefully  in  environmental  studies. 
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Primer  Sensitivity  in  amoA  qPCR 

This  work  represents  the  first  study  to  experimentally  test  congruity  of  the 
published  primer  sets  for  the  archaeal  amoA  gene  on  environmental  samples.  Despite  the 
different  target  regions  and  degeneracies  of  the  primer  sets,  they  yielded  remarkably 
similar  results  in  the  upper  2  meters  of  the  sediment  column  (Figure  1C).  The  main 
difference  between  the  primer  sets  was  observed  in  the  deeper  AOA  communities,  where 
the  Treusch  et  al.  (Treusch  et  al.  2005)  primers  underestimated  abundance  by  1  order  of 
magnitude  relative  to  the  other  primer  sets  (Figure  1C).  The  Treusch  et  al.  (Treusch  et  al. 
2005)  primers  were  initially  designed  for  use  in  a  TaqMan  assay  using  the  amoA  genes 
detected  in  a  soil  fosmid  54d9  (61)  as  well  as  the  Sargasso  Sea  Whole  Genome  Shotgun 
(WGS)  sequence  (Treusch  et  al.  2005),  and  it  should  be  noted  that  their  use  here  involved 
direct  application  in  a  SYBR  Green  I  assay,  rather  than  the  intended  TaqMan  assay. 
However,  given  the  larger  degree  of  degeneracy  of  these  forward  and  reverse  primers  (six 
degeneracies  in  the  forward  primer  and  five  degeneracies  in  the  reverse  primer)  (Figure 
3),  we  would  have  expected  a  broader  target  range  for  these  primers  relative  to  other  sets 
rather  than  a  more  restricted  target  range.  To  examine  the  possibility  that  the  Treusch  et 
al.  (Treusch  et  al.  2005)  primers  were  missing  a  subset  of  the  deep  Waquoit  Bay  AOA 
community,  these  primers  were  checked  in  silico  against  sequences  recovered  in  libraries 
constructed  using  the  Francis  et  al.  (2005)  primers.  The  alignment  showed  that  these 
primers  matched  only  5.7%  of  our  recovered  sequences  (4  of  69  clones  from  4  m). 
Furthermore,  matching  the  Treusch  et  al.  (Treusch  et  al.  2005)  primers  to  archaeal  amoA 
reference  sequences  (14  of  318  sequences  reported  in  Francis  et  al.  (Francis  et  al.  2005)) 
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and  a  survey  of  archaeal  amoA  genes  from  GenBank  (132  of  3530  sequences  from 


GenBank)  also  reveals  that  these  primers  match  less  than  4%  of  the  sequences.  Therefore, 
it  is  not  surprising  that  this  primer  set  returns  a  lower  abundance  of  amoA  than  data 
derived  from  the  other  amoA  gene  primers.  As  a  result  of  these  findings,  we  have  chosen 
to  report  our  abundances  as  determined  using  the  Francis  et  al.  (Francis  et  al.  2005) 
primer  set  for  the  following  discussion. 

Relative  abundances  ofAOB  andAOA 

For  many  years,  AOB  were  the  primary  organisms  believed  able  to  catalyze  the 
nitrification  processes  that  can  ultimately,  through  coupled  nitrification-denitrification, 
lead  to  nitrogen  removal  from  the  environment.  Studies  have  shown  AOB  abundance 
(Stehr  et  al.  1995;  Cebron  et  al.  2003)  and  nitrification  rates  (Rysgaard  et  al.  1999; 

Cebron  et  al.  2003)  increasing  in  the  freshwater  end  members  of  surficial  estuaries 
compared  with  abundances  and  rates  measured  in  the  seawater  end-members.  However, 
these  studies  did  not  include  AOA  cell  abundance,  or  their  contributions  to  nitrification 
rates.  Recent  reports  have  examined  spatial  patterns  in  the  distribution  of  the  AOA  and 
AOB  amoA  genes  in  relation  to  ammonia  and  organic  carbon  concentrations,  salinity 
fluctuations,  and  other  chemical  loads  (Bernhard  et  al.  2005;  Francis  et  al.  2005;  Freitag 
et  al.  2006;  Mosier  and  Francis  2008;  Bernhard  et  al.  2010).  Relative  changes  in 
abundance  between  AOB  and  AOA  have  been  found  to  be  correlated  with  C:N  (Mosier 
and  Francis  2008),  salinity  and  oxygen  (Beman  and  Francis  2006;  Mosier  and  Francis 
2008;  Santoro  et  al.  2008),  as  well  as  trace  metal  availability  (lead  and  nickel)  (Caffrey  et 
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al.  2007;  Mosier  and  Francis  2008).  The  recent  report  of  potential  mixotrophy  in  C. 
symbiosum  (Hallam  et  al.  2006)  and  planktonic  archaea  (Ouverney  and  Fuhrman  2000; 
Ingalls  et  al.  2006)  hint  at  the  potential  plasticity  of  the  metabolic  capabilities,  which  may 
allow  AOA  to  survive  in  dynamic  regions  of  the  environment  were  DOC  and  mobilized 
ammonia  may  periodically  be  injected.  Recent  culture  based  studies  have  shown  that  the 
archaea  are  adapted  to  life  in  nutrient  limited  conditions.  Nitrosopumilus  maritimus 
SCMl  has  been  shown  to  have  a  high  affinity  for  ammonium  and  a  low  substrate 
threshold  of  10  nM  NH3  or  less  (Martens -Habbena  et  al.  2009).  The  most  dynamic 
conditions  in  the  Waquoit  Bay  subterranean  estuary  are  found  in  the  upper  and  lower  salt 
transition  zones,  which  are  periodically  exposed  to  saline,  low  oxygen  pore  waters  with 
elevated  concentrations  of  ammonium,  and  in  some  cases  DOC  (Charette  and  Sholkovitz 
2006;  Bone  et  al.  2007).  However,  a  clear  consensus  on  the  interactions  between  these 
environmental  factors  and  AOA  distributions  has  not  been  reached  (Erguder  et  al.  2009). 

In  the  Waquoit  Bay  aquifer  system,  archaeal  amoA  dominated  the  total  amoA 
gene  pool,  accounting  for  >  90%  of  total  amoA  at  most  depths  in  the  core.  Peaks  in 
archaeal  amoA  distribution  occurred  in  the  surface  meter  and  in  the  lower  portion  of  the 
nitrate  peak  near  the  saltwater  interface,  while  the  bacterial  amoA  abundance  peaked  at  a 
shallower  depth  within  the  fresh  groundwater  plume.  In  the  present  study,  neither  the  raw 
nor  relative  abundance  of  sediment- associated  AOA  or  AOB  showed  a  strong  correlation 
with  any  of  the  simultaneously  measured  pore  water  parameters  (O2,  salinity,  NO3',  NO2’, 
NH4‘^,  pH,  or  redox  potential).  The  best  correlation  was  observed  between  archaeal  amoA 
and  redox  potential.  However,  fluctuations  in  pore  water  conditions  (the  depth  of  the 
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deep  salt- wedge,  the  presence  of  a  reduced  ammonium  bearing  plume  in  center  of  the 
groundwater  plume  and  tidal  pumping)  have  been  observed  or  modeled  at  nearby  sites 
(Charette  et  al.  2005;  Kroeger  and  Charette  2008;  Spiteri  et  al.  2008),  and  the  distribution 
of  sediment- associated  microbes  may  reflect  an  average,  rather  than  instantaneous, 
condition. 


Conclusions 

This  study  reports  the  diversity  and  distribution  of  bacterial  and  archaeal  16S 
rRNA  and  amoA  genes  within  a  subterranean  estuary.  We  also  report  on  the  congruity  of 
various  published  primer  sets  used  for  amplifying  archaeal  amoA  genes.  Closer 
inspection  of  the  abundance  and  diversity  of  nitrifying  organisms  inhabiting  the  aquifer 
sediments,  as  defined  by  the  amoA  gene,  indicated  a  community  dominated  by  AO  A 
rather  than  AOB  at  most  depths,  except  between  1.2  and  2  m  in  the  fresh  groundwater 
plume.  Sequencing  the  amoA  gene  from  the  two  peaks  in  archaeal  abundance  revealed 
that  these  represent  two  different  populations.  The  shallow  population  contained  amoA 
gene  sequences  within  the  Water  Column  A  cluster,  with  the  most  closely  related 
sequences  originating  from  low  salinity  estuarine  sediments.  The  deeper  aquifer 
population  is  more  closely  related  to  amoA  sequences  recovered  from  soil  and 
groundwater  systems.  This  difference  in  community  structure  suggests  an  environmental 
factor  (or  combination  of  factors)  may  control  the  distribution  of  these  populations  within 
the  aquifer.  However,  the  measured  environmental  conditions  in  this  study  do  not  provide 
an  adequate  answer  for  this  question.  Further  work,  both  in  the  environment  and  with 
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cultured  AOA,  is  needed  to  discover  what  environmental  eonditions  favor  different  AOA 


and  AOB  populations  and  what  other  metabolie  and  biogeoehemical  functions  AOA  may 
bring  to  the  environment. 
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Table  1.  Molecular  Methods 


Target 

Gene 

Primer  Names 

And 

Sequences 

Thermal  Program 

qPCR 

Efficiency 

Bacterial 

16S  rRNA 

(sequencing 

only) 

8F  (400nM)- 

AGAGTTTGATCCTGGCTCAG 

1492R  (400nM)- 
GGTTACCTTGTTACGACTT 

94°C  for  120sec. 

30  cycles  of 

94°C  for  30sec 

47°C  for  60sec 

72°C  for  60sec 

Final  extension  of 
72°C  for  600sec 

Bacterial 

16S  rRNA 

Bad 369  (lOOOnM)- 
CGGTGAATACGTTCYCGG 

Prokl541  (lOOOnM)- 
AAGGAGGTGATCCRGCCGCA 

50°C  for  120sec 

94°C  for  120sec 

40  cycles  of 

94°C  for  ISsec 

59°C  for  60sec 

88-90% 

Archaeal 

16S  rRNA 

(sequencing 

only) 

21Fa  (400nM)- 

TTCCGGTTGATCCYGCCGGA 

1492R  (400nM)- 
GGTTACCTTGTTACGACTT 

94°C  for  120sec. 

30  cycles  of 

94°C  for  30sec 

47°C  for  60sec 

72°C  for  60sec 

Final  extension  of 
72°C  for  600sec 

Archaeal 

16S  rRNA 

Arcmixl369F  (1:1  mixture) 

(SOOnM)- 

CGGTGAATACGTCCCTGC  and 
CGGTGAATATGCCCCTGC 

Prokl541  (lOOOnM)- 
AAGGAGGTGATCCRGCCGCA 

50°C  for  120sec 

94°C  for  120sec 

40  cycles  of 

94°C  for  ISsec 

59°C  for  60sec 

80°C  for  ISsec 
(read) 

98-10S% 

MGl 

Crenarchaeota 

GI-751F  (SOOnM)- 
GTCTACCAGAACAYGTTC 
GI_956R  (SOOnM)- 
HGGCGTTGACTCCAATG 

9S°C  for  900sec 

SO  cycles  of 

94°C  for  ISsec 

S8°C  for  30sec 

72°C  for  30sec 

80°C  for  ISsec 
(read) 

93-9S% 

Francis  et  al. 
amoA 
(sequencing 
and  qPCR) 

amoAF  (SOOnM)- 
STAATGGTCTGGCTTAGACG 

amoAR  (SOOnM)- 
GCGGCCATCCATCTGTATGT 

9S°C  for  300 

30  cycles  of 

94°C  for  30sec 

S3°C  for  60sec 

72°C  for  60sec 

80°C  for  ISsec 
(read) 

93-9S% 

Wuchter  et  al. 
amoA 

archamoAF  (SOOnM)- 
CTGAYTGGGCYTGGACATC 

archamoR  (SOOnM)- 
TTCTTCTTTGTTGCCCAGTA 

9S°C  for  300sec 

3S  cycles  of 

94°C  for  30sec 

61.S°Cfor  40sec 

9S-96% 

Positive  Control 

Reference 

Escherichia  coli. 

Pseudomonas 

stutzeri 

P.  stutzeri 

(Suzuki  et 
al.  2000) 

Ferroplama 

acidarmonas 

(Reysenbac 
h  et  al. 

1992) 

F.  acidarmonas,  N. 
maritimus 

(Suzuki  et 
al.  2000) 

N.  maritimus, 
environmental  clone 

(Mincer  et 
al.  2007) 

N.  maritimus, 
environmental  clone 

(Francis  et 
al.  200S) 

N.  maritimus, 
environmental  clone 

(Wuchter  et 
al.  2006) 

Target 

Gene 

Primer  Names 

And 

Sequences 

Thermal  Program 

qPCR 

Efficiency 

Positive  Control 

Reference 

72°C  for  40sec 

80°C  for  ISsec 
(read) 

Mincer  et  al. 
amoA 

crenamoAmodF  (500nM)- 
TGGCTAAGACGMTGTA 

crenamoA-Q-F  (500nM)- 
GCARGTMGGWAARTTCTAYAA 

crenamoAmodR  (500nM)- 
AAGCGGCCATCCATCTGTA 

95°C  for  900sec 

50  cycles  of 

94°C  for  15sec 

52°C  for  30sec 

72°C  for  30sec 

80°C  for  15sec 
(read) 

95-96% 

N.  maritimus, 
environmental  clone 

(Mincer  et 
al.  2007) 

Treusch  et  al. 
amoA 

Amol96F(1000nM)- 

GGWGTKCCRGGRACWGCMAC 

Amo277R(1000nM)- 

CRATGAAGTCRTAHGGRTADC 

95°C  for  ISOsec 

45  cycles  of 

95°C  for  15sec 

50°C  for  40sec 

80°C  for  15sec 
(read) 

90-91% 

N.  maritimus. 

Environmental 

clone 

(Treusch  et 
al.  2005) 

Bacterial 

amoA 

amol89F  (300nM)- 
GGHGACTGGGAYTTCTGG 

amoA2-R’  (900nM)- 
CCTCKGSAAAGCCTTCTTC 

50°C  for  120sec 

95°C  for  600sec 

40  cycles  of 

95°C  for  45  sec 

55°C  for  60sec 

72°C  for  45sec 

80°C  for  15sec 
(read) 

96-97% 

Nitrosomonas  sp. 

(Okano  et 
al.  2004) 
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Figure  1.  Pore  water  chemistry  at  study  site.  A)  Salinity,  nutrients  (NO3'  and  NH4"^),  and 
oxygen  concentrations  in  the  sediment  pore  water  at  the  time  of  coring.  The  aquifer 
contains  three  distinct  parcels  of  water:  i)  a  saltwater  surface  zone  (to  about  1  m),  with 
low  oxygen  and  nitrate  levels,  ii)  a  persistent  fresh  groundwater  plume  that  carries 
elevated  oxygen  and  nitrate  concentrations,  and  iii)  a  deep  anoxic  salt  wedge,  which 
introduces  ammonia  in  to  the  system.  B)  Depth  distribution  of  total  bacterial,  archaeal 
and  Marine  Group  I  (MGl)  archaeal  16S  rRNA  gene  abundance.  C)  Primer- specific 
abundance  of  archaeal  amoA  genes  using  the  four  published  primer  sets  (Table  1).  D) 
Depth  distribution  of  bacterial  and  archaeal  amoA  abundance  relative  to  MGl  16S  rDNA 
abundance. 
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Redox  potential  (mV) 

Figure  2.  The  abundance  of  MGl  16S  rRNA  and  archaeal  amoA  genes  are  weakly 
correlated  with  redox  potential.  No  other  measured  environmental  parameter  was 
strongly  correlated  with  the  distribution  of  these  genes. 
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Figure  3.  Alignment  of  each  of  the  primer  pairs  used  in  the  study  to  the  reference 
sequences:  Sargasso  Sea  WGS  (Venter  et  al.  2004),  Soil  fosmid  54d9  (Treusch  et  al. 
2005),  and  V.  maritimus  genome  (gene  1500)  (Konneke  et  al.  2005).  Numbers  beneath 
primer  sequence  represent  alignment  position  relative  to  Sargasso  Sea  sequence. 
Mismatches  to  each  of  the  three  reference  sequences  are  depicted  above  primer  sequence. 
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Figure  4.  Phylogenetic  tree,  based  on  1309  bp,  of  the  archaeal  16S  rRNA  gene.  Type 
strains  are  represented  with  bold  and  italicized  font.  Sequences  recovered  from  this  study 
are  in  bold.  Sequences  that  were  repeatedly  recovered  are  marked  with  *.  ARB  neighbor 
joining  methods,  with  Jukes-Cantor  corrections,  were  used  for  tree  construction. 
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Soils  and  Sediments  Water  Column 


Figure  5.  Phylogenetic  tree  based  on  456  bp  of  unambiguous  amoA  nucleotide  sequence. 
ARB  neighbor  joining  methods,  with  Jukes-Cantor  corrections,  were  used  for  tree 
construction.  Clone  libraries  were  constructed  at  six  depths  within  the  aquifer  using  the 
Francis  et  al.  (Francis  et  al.  2005)  primer  pair.  Recovered  sequences  fall  within  two  major 
grouping  of  the  archaeal  amoA  tree.  Sequences  from  above  a  depth  of  2  m  in  the  aquifer 
are  generally  most  similar  to  low  salinity  water  column  and  estuary  sequences  with 
closest  relatives  recovered  from  the  northern  section  of  San  Francisco  Bay.  Sequences 
from  2  m  and  below  are  more  closely  related  to  sequences  recovered  from  freshwater 
sediments,  rhizosphere,  groundwater  and  terrestrial  soils.  Numbers  within  the  wedges 
represent  the  number  of  sequences  included  in  that  branch. 
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Abstract 


Isotopic  and  molecular  techniques  were  combined  to  investigate  where 
denitrification  may  occur  within  a  transect  through  a  coastal  aquifer  prior  to  discharge  of 
the  groundwater  into  Waquoit  Bay,  MA.  The  sampled  aquifer  shows  three  distinct  parcels 
of  water  each  with  its  own  nutrient,  oxygen  and  salinity  characteristics.  The  surface  salt 
zone,  observed  in  the  upper  meter  of  the  sediment  column,  was  characterized  by  a  nitrate 
pool  enriched  in  both  (up  to  +  20.1%o)  and  (up  to  +  11.7%o),  characteristic  of 
groundwater  affected  by  denitrification.  The  near-surface  sediments  also  harbored  a 
denitrifying  community  containing  the  nitrite  reductase  genes  nirK  and  nirS.  This 
contrasts  with  the  deep  salt  wedge,  where  a  more  abundant  nirS  gene  characterized  the 
denitrifying  community.  Interestingly,  the  large  nitrate  pool  in  the  fresh  groundwater 
parcel  appears  to  advect  seaward  without  substantial  change  along  the  transect. 
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Introduction 


Denitrification  is  an  important  mechanism  for  nitrogen  removal  in  estuaries  and 
coastal  sediments  (Seitzinger  et  al.  2006).  More  than  half  of  the  inorganic  nitrogen 
entering  estuaries  and  coastal  ecosystems  may  be  removed  by  this  microbially-mediated 
process  (Seitzinger  1988).  Fresh  water  may  also  enter  the  coastal  ocean  through 
submarine  groundwater  discharge  (Moore  1999;  Burnett  and  Dulaiova  2003),  and  the 
amount  of  inorganic  nitrogen  introduced  through  SGD  could  be  influenced  by  microbial 
processes  occurring  within  the  subterranean  estuary.  Previous  studies  have  shown  that  the 
chemical  load  of  nutrients,  trace  metals,  organic  carbon,  and  methane  (Bugna  et  al.  1996; 
Paerl  1997;  Cai  et  al.  2003)  discharged  from  ground  waters  can  rival  that  of  riverine 
input,  although  the  total  water  flux  may  be  less  (Burnett  et  al.  2003;  Slomp  and  Van 
Cappellen  2004;  Kroeger  et  al.  2007).  Unlike  many  coastal  or  estuarine  sediments,  where 
water  flows  over  the  sediments  with  only  tidal  and  pressure  driven  ventilation, 
groundwater  flows  through  aquifer  sediments  resulting  in  a  much  larger  exposure  of  the 
water  parcel  to  reactive  sediment  surfaces  and  associated  microbial  communities. 
Groundwater  is  also  isolated  from  the  atmosphere,  which  allows  for  gradual  depletion  of 
oxygen  along  its  flow  path.  These  important  environmental  differences  between  overland 
and  groundwater  transport  may  contribute  to  different  microbial  consortia  and  nitrogen 
dynamics  within  aquifer  sediments  compared  with  more  traditional  riverine  or  estuarine 
sediments. 

Denitrifying  microbes  are  phylogentically  diverse,  and  are  found  in  all  three 
domains  of  life,  though  mostly  within  the  prokaryotes.  In  addition,  many  denitrifiers  are 


75 


closely  related  to  non-denitrifying  microbes,  making  16S  rRNA  genes  an  insufficient 
genetic  marker  for  this  metabolism  (Zumft  1992;  Braker  et  al.  2000).  Genetic  markers 
commonly  used  to  assess  the  abundance  and  diversity  of  denitrifying  bacteria  are  the 
genes  for  the  dissimilatory  nitrite  reductase  enzymes  that  catalyze  the  reduction  of  nitrite 
to  nitric  oxide,  the  first  gaseous  product  in  the  denitrification  pathway  (Ward  1996; 
Braker  et  al.  2000).  Nitrite  reductase  occurs  in  at  least  two  major  forms,  a  multi-heme 
cytochrome  containing  enzyme  (encoded  by  the  nirS  gene)  and  a  copper-containing 
protein  (encoded  by  the  nirK  gene)  (Zumft  1997).  Recent  work  has  used  these  genes  to 
describe  the  distribution  and  diversity  of  denitrifying  microbes  in  coastal  sediments  and 
groundwater  systems  (Braker  et  al.  1998;  Nowicki  et  al.  1999;  Braker  et  al.  2000; 

Nogales  et  al.  2002;  Santoro  et  al.  2006). 

The  presence  of  an  organism,  or  a  gene,  does  not  necessarily  indicate  activity, 
however,  particularly  for  a  facultative  process  such  as  denitrification.  Additional 
measures  are  needed  to  determine  where  and  when  denitrification  is  occurring.  Hence, 
stable  isotope  tracer  techniques  have  been  applied  to  many  different  environments  to 
detect  denitrification  (Fukada  et  al.  2004).  In  addition,  changes  to  the  natural  abundance 
stable  isotopic  signatures  of  nitrogen  oxides  can  indicate  potential  locations  of 
denitrification  (Aravena  and  Robertsona  1998).  Biological,  chemical  and  physical 
processes  can  fractionate  compounds  due  to  isotopic  differences  between  molecules.  This 
fractionation  can  yield  information  about  the  reaction  history  of  the  compound.  In 
particular,  with  nitrate  (NO3')  both  the  nitrogen  and  oxygen  atoms  can  be  used  to  track 
isotopic  fractionations  and  distinguish  removal  of  NO3'  via  denitrification  from  dilution 
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and  mixing  (Bottcher  et  al.  1990).  The  goal  of  this  study  was  to  investigate  the  removal 
of  nitrate  within  a  subterranean  estuary  on  Cape  Cod,  MA  by  identifying  regions  of 
denitrification  and  characterizing  the  microbial  communities  in  these  regions,  using  a 
combination  of  molecular  and  stable  isotopic  techniques.  Given  that  nutrient  input 
through  SGD  may  cause  increased  eutrophication  of  the  coastal  ocean  (Valiela  et  al. 
1992;  Paerl  1997;  Bricker  et  al.  1999),  it  is  important  to  understand  how  these  fluxes  are 
mediated  from  both  chemical  and  microbiological  perspectives. 

Methods 

Sample  Collection 

During  June  and  July  2008  pore  water  profiles  were  obtained  along  transects 
normal  to  the  shoreline  through  the  intertidal  zone.  The  transect  was  composed  of  four 
sites  beginning  with  a  site  1  m  above  mean  high-tide  and  extending  through  the  intertidal 
zone  (Figure  1).  Pore  water  was  collected  using  a  piezometer  with  a  stainless  steel 
“Retract-a-tip”  (AMS,  Idaho),  and  vertical  resolution  was  as  small  as  0.15  m.  After 
flushing  with  at  least  one  liter  of  pore  water  (more  than  6-fold  flushing  of  tubing 
volume),  pH,  dissolved  oxygen,  salinity  and  temperature  were  measured  using  a  YSI 
600R  multi-probe  (Charette  and  Sholkovitz  2002).  Oxygen,  pH,  and  redox  sensors  were 
calibrated  against  standard  solutions  following  the  manufacturers  protocols  prior  to 
deployment  in  the  field.  Triplicate  pore  water  samples  for  isotopic  (NO3',  NO2’)  and 
nutrient  analyses  (NO3',  NO2 ,  NH4‘^)  were  collected  at  a  resolution  of  0.15  m  to  0.30  m 
through  the  top  of  the  deep  salt  wedge  (3-4  m),  filtered  through  a  0.2  pm  pore-size 
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Sterivex  filter,  kept  on  ice  and  either  analyzed  upon  return  to  laboratory  or  frozen  for 
later  analysis. 

After  pore  water  sampling  (on  a  falling  tide),  a  sediment  core  (near  low  tide)  was 
collected  at  site  B  using  a  hand-operated  bailer-boring  auger  (Eijkelkamp,  The 
Netherlands).  The  recovered  sediment  was  sampled  using  aseptic  techniques.  Sediment 
samples  for  DNA  analyses  (ca.  1  g  and  5  g  of  sediment)  were  frozen  in  liquid  nitrogen  for 
later  processing. 

Nutrient  concentration  and  isotopic  analyses 

Nitrite  concentrations  were  measured  spectrophometrically  after  Pei  et  al.  (1990). 
Nitrate  plus  nitrite  (NOx)  concentrations  were  measured  by  chemiluminescence  following 
a  hot  vanadium  reduction  (Cox  1980;  Garside  1982;  Garside  1985).  Ammonium  was 
measured  using  the  phenolhypochlorite  method  (Solorzano  1969). 

The  nitrogen  and  oxygen  isotope  ratios  of  nitrate  (expressed  in  delta  notation  as 
6'^Nno3  and  8^^0no3,  respectively)  were  measured  after  conversion  to  N2O  using  the 
denitrifier  method  (Sigman  et  al.  2001;  Casciotti  et  al.  2002).  8^^Nno3  and  8^^0no3  values 
are  reported  in  %o  vs.  air  and  Vienna  Standard  Mean  Ocean  Water  (VSMOW), 
respectively,  and  are  calibrated  to  these  scales  by  analyzing  20  nanomole  (nmol)  aliquots 
of  the  potassium  nitrate  reference  materials  USGS32,  USGS34,  and  USGS35  after  every 
9  samples.  Samples  were  also  analyzed  at  the  20  nmol  level  to  match  the  standards.  All 
isotopic  measurements  were  made  using  a  Thermo-Finnigan  Delta^^^^  XP  isotope-ratio 
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mass  spectrometer  with  a  custom-built  continuous  flow  inlet  similar  to  that  described 
previously  (Casciotti  et  al.  2002). 

DNA  extraction 

DNA  was  extracted  from  the  aquifer  sediments  using  lysozyme  (5  mg  mL"',  Ihr, 
37°C)  and  proteinase  K  (2  mg  mL'\  2  hr,  60°C)  digestions  in  lysis  buffer  (200  mM  Tris- 
HCl,  4  %  sodium  dodecyl  sulfate,  100  mM  NaCl,  100  pg  mL'^  polyadenosine,  pH  8.0) 
followed  by  three  to  five  iterations  of  3  minutes  in  a  dry-ice/ethanol  bath  and  5  minutes 
60  °C  water  bath,  DNA  was  purified  using  three  phenohchloroformdsoamyl  alcohol 
separation  and  final  isopropanol  precipitation.  Recovered  DNA  was  resuspended  in  50 
pL  of  TE  for  all  subsequent  analyses. 

Quantitative  PCR 

Quantitative  PCR,  using  an  iCycler  thermocycler  (Bio-Rad,  Hercules,  CA),  was 
employed  to  determine  the  gene  copy  numbers  of  nirK  and  nirS  at  0.30  m  resolution 
throughout  the  depth  of  the  core.  Triplicate  DNA  extractions  of  each  sample  were  diluted 
to  lOng  pp'  and  amplified  using  either  the  iQ  SYBR  green  Supermix  (Bio-Rad)  or  a 
cocktail  consisting  of  1  X  colorless  master  mix  (Promega,  Madison,  WI)  amended  with 
10  nM  dNTPs  (Promega),  0.5  X  SBYR-Green  I  (Invitrogen,  Carlsbad  CA),  20  nM 
fluorescein  (Sigma),  1  mM  MgCl2  (Promega),  forward  and  reverse  primer  (Table  1,  300- 
1000  nM)  and  1.25  units  of  GoTaq  Flexi  polymerase  (Promega).  All  qPCR  amplifications 
were  repeated  on  three  separate  days  using  optimized  thermal  programs  and  resulted  in  a 
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PCR  efficiency  between  88-92  %  (Table  1  and  references  therein).  Standard  curves  were 
constructed  from  serial  dilutions  (10^  to  10^  copies  pL'^  plus  a  water  blank)  of  linearized 
p-Gem  T-Easy  plasmids  (Promega)  containing  inserts  of  the  functional  gene  of  interest 
(Table  1).  Plasmids  (0.5  pg  pL'^)  were  linearized  by  digestion  with  the  EcoRI  (New 
England  Biolabs)  restriction  enzyme  in  the  supplied  buffer  (2  hr  at  37  °C,  heat 
inactivation  20  min.  at  65  °C).  Restriction  digests  were  cleaned  using  the  SV  Wizard  PCR 
purification  kit  (Promega),  eluted  in  30  pE)  and  quantified  using  a  Nanodrop 
spectrophotometer  (Thermo-Eisher).  Melt  curves  were  performed  on  all  reactions  to 
confirm  specificity. 


PCR,  Cloning  and  Sequencing 

Eor  each  target  gene  (nirK  and  nirS),  three  replicate  PCR  amplifications  from  1  m 
(shallow  salt  zone),  2  m  (freshwater  peak)  and  3  m  (deep  salt  wedge)  were  conducted 
according  to  published  protocols  (Table  1  and  ref.  therein).  The  replicate  amplification 
products  were  pooled  and  cleaned  with  the  SV  Wizard  PCR  Clean-up  (Promega) 
following  the  manufacturers  protocols  and  eluted  in  50  pE  Tris-EDTA  (TE)  buffer. 
Cleaned  amplicons  were  ligated  using  the  p-Gem  T-Easy  vector  kit  (Promega)  using  the 
manufacturers  protocols  and  transformed  into  chemically  competent  E.  coli  JM109  using 
standard  procedures  (Seidman  et  al.  2002).  The  transformants  were  plated  on  Euria  Broth 
plates  containing  IPTG/X-gal/ampicillin  (10  pmols,  1  mg  and  50  pg  mE'^  respectively) 
and  incubated  at  37  °C  for  24  hr.  Blue/white  screening  was  used  to  pick  individual 
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colonies,  which  were  transferred  to  growth  blocks  containing  Superbroth  (Difco)  with 
ampicillin  (50  pg  mL'^)  and  grown  shaking  (250  rpm)  overnight.  The  overnight  cultures 
were  harvested  by  centrifugation  (2800  rpm  for  8  min.  at  4  °C).  Plasmid  purification  and 
sequencing  were  performed  at  the  Keck  Facility,  Josephine  Bay  Paul  Center  at  the 
Marine  Biological  Laboratories  (JBPC-MBL,  Woods  Hole,  MA)  using  a  BiomekFX 
workstation  (Beckman  Coulter)  and  a  3730XL  capillary  sequencer  (Applied  Biosystems). 

Sequencing  chemistry  was  a  modified  Big-Dye  reaction  containing  0.5  pi  of  Big- 
Dye  Terminator,  0.4  pL  of  primer  (M13R  5’-CAGGAAACAGCTATGAC-3’,  15  pM), 
0.1  pL  dimethyl  sulfoxide,  1/3  X  reaction  buffer,  template  (200-400  ng  DNA)  and  water 
to  a  final  volume  of  6  pL  per  reaction.  A  thermal  program  of  60  cycles  of  96  °C  for  10 
sec.,  50  °C  for  5  sec.  and  60  °C  for  4  min.  was  used.  Following  reaction,  the  DNA  was 
precipitated  by  a  series  of  75  %  and  70  %  isopropanol  washes,  air  dried  and  finally 
resuspended  in  7  pL  Hi-Di  formamide  and  frozen  until  sequencing. 

Sequence  analysis 

Successful  sequencing  reads  were  checked  and  the  vector  trimmed  at  the  JBPC- 
MBL  (Woods  Hole,  MA).  Base  calling  and  sequence  quality  was  performed  by  the  Phred 
and  Phrap  software  packages  (Ewing  and  Green  1998;  Ewing  et  al.  1998).  The  Eucy 
software  (Chou  and  Holmes  2001)  was  used  for  vector  trimming.  These  processed 
sequences  were  searched  against  the  non-redundant  nucleotide  (nr/nt)  database  using  the 
blastn  algorithm.  16S  rRNA  sequences,  and  related  GenBank  sequences,  were  imported 
into  ARB,  aligned  with  the  ARB  aligner  and  manually  checked  (Eudwig  et  al.  2004). 
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Conservative  column  filters  were  constructed  masking  any  sites  where  data  was  missing. 
Trees  were  constructed  on  the  aligned  sequences  using  the  ARB  neighbor  joining 
function.  NirS  and  nirK  gene  sequences  were  treated  in  the  same  manner  using  databases 
built  from  nirS  and  nirK  sequences  from  GenBank. 

GenBank  Submission 

All  sequencing  generated  and  used  in  this  study  are  reported  to  GenBank  under 
the  accession  numbers  HM160368-HM 160497. 

Results 

Nutrient  Chemistry 

The  distributions  of  salinity  and  nutrients  in  the  near-shore  aquifer  at  the  head  of 
Waquoit  Bay  have  been  described  in  detail  in  previous  studies  of  the  system  (Charette  et 
al.  2005;  Charette  and  Sholkovitz  2006;  Kroeger  and  Charette  2008).  In  short,  the  aquifer 
has  three  major  features:  1)  a  surface  salt  zone  (salinity  of  4-6)  driven  by  tidal  and/or 
wave  action,  2)  a  oxidized  freshwater  plume  and  3)  a  reduced  deep  salt  wedge.  The 
intrusion  of  deep  salt  wedge  beneath  the  seaward-flowing  freshwater  plume  results  in  a 
shoaling  of  the  ground  water  plume  and  a  deep  salt  wedge  in  the  seaward  stations  in  the 
transect  (Figure  2A).  A  reducing,  ammonium  and  iron  rich  freshwater  plume  was  never 
observed  at  the  current  study  site  but  has  been  reported  previously  at  study  sites  on  the 
eastern  portion  of  the  bay  studies  at  a  site  east  of  the  site  (see  Charette  et  al.  2005, 
Kroeger  and  Charette  2008). 


82 


Nitrate  is  the  dominant  species  of  inorganic  nitrogen  (>350  pM)  in  the  freshwater 
portion  of  the  aquifer.  Nitrate  decreases  to  concentrations  as  low  as  12.5  pM  within  the 
surface  brackish  layer  and  is  undetectable  in  the  deep  salt  wedge  (Figure  2C).  The 
oxygen  distribution  is  similar  to  the  nitrate  distribution,  with  elevated  concentrations 
within  the  freshwater  plume  and  decreasing  towards  the  two  saline  water  parcels.  In  the 
surface  salt  zone,  oxygen  decreases  to  ca.  20  pM.  A  similar  decrease  in  oxygen  is  seen 
between  the  freshwater  plume  and  the  deep  salt  wedge  (Figure  2B).  On  the  other  hand, 
ammonium  is  present  only  in  the  deep  salt  wedge  at  the  time  of  our  sampling  (Figure  2D) 
and  is  consistent  with  organic  matter  remineralization  in  anoxic  estuarine  sediments 
(Talbot  et  al.  2003;  Kroeger  and  Charette  2008). 

The  lowest  isotope  ratios  for  both  8'^Nno3  and  no3  are  found  within  the  core 
of  the  freshwater  plume,  with  8^^Nno3  values  of  4-5%o  (Figure  3A)  and  8^^0no3  values  of 
2-3%o  (Figure  3B).  Both  8^^Nno3  and  8^^0no3  increased  into  the  surface  salt  zone  as 
nitrate  concentrations  decreased  (Figure  3 A,  3B,  2C).  These  patterns  are  consistent  with 
nitrate  removal  by  a  process  with  normal  fractionation,  rather  than  through  advection  or 
dilution  that  ought  not  result  in  fractionation.  Nitrate  isotopic  measurements  within  the 
deep  salt  wedge  were  not  possible  because  sampling  resolution  within  the  steep  gradient 
in  nitrate  concentrations  did  not  yield  samples  within  our  analytical  limit  (Figure  2C). 

The  8^^0h2o  was  also  measured  throughout  the  transect  in  order  to  constrain  the 
d^^O  value  of  local  nitrate  production  (Figure  3C).  8^^0h2o  values  correlated  with  salinity 
(r^=  0.71)  with  the  most  depleted  values  (-8%o)  occurring  within  center  of  the  freshwater 
plume  (Figure  3D).  The  8^^0h2o  increased  by  2%o  moving  up  from  the  center  of  the 
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freshwater  plume,  and  approached  -6%o  in  the  surface  brackish  zone.  The  gradient  in 
6'^Oh2o  was  greater  moving  from  the  center  of  the  plume  toward  the  deep  salt  wedge, 
with  values  as  high  as  -1.8%o  in  the  salt  wedge,  at  salinities  of  26.2  (Figure  3C,  D). 

Nitrite  reductase  abundance  and  distribution 

In  the  Waquoit  Bay  aquifer,  both  nirS  and  nirK  were  detected  at  all  depths  (Figure 
5).  Of  the  two  genes,  nirK  displayed  less  variability  with  depth,  ranging  only  from  1.53  x 
lO"^  to  1.55  X  10^  copies  per  gram  sediment  throughout  the  core.  The  maximum 
abundance  for  nirK  occurred  at  a  depth  of  2  m  while  the  minimum  was  found  at  1  m.  The 
abundance  of  nirK  in  the  aquifer  sediments  is  similar  to  abundances  reported  for  soils, 
more  specifically  low  carbon  and  low  nitrogen  soils  (Henry  et  al.  2004).  In  contrast,  nirS 
abundances  varied  by  more  than  two  orders  of  magnitude,  from  7.83  x  10  copies  per 
gram  sediment  at  1  m  depth  to  6.8  x  10^  copies  per  gram  sediment  at  3.8  m  depth.  These 
abundances  are  similar  in  range  to  reports  for  the  Colne  estuary  sediments  (United 
Kingdom)  where  the  highest  abundances  were  observed  at  the  head  of  the  estuary  and 
copy  numbers  decreasing  moving  toward  the  estuary  mouth  (Smith  et  al.  2007). 

Bacterial  16S  rRNA  Gene  Diversity 

Phylogenetic  analysis  of  the  bacterial  16S  rRNA  genes  designated  11  operational 
taxonomic  units  using  a  cutoff  of  97%  sequence  identity.  The  Proteobacteria  division 
dominated  the  bacterial  community,  with  members  from  the  gamma-  and  alpha- 
subdivisions  accounting  for  43%  and  28%  of  the  population  respectively.  Lesser 


84 


contributions  were  from  the  Actinobacteria-,  Firmicute-  and  Planctomycetales-hke 
organisms. 

The  genes  recovered  from  the  gamma-Proteobacteria  division  were  recovered 
from  all  depths  of  the  core.  Alcanivorax-Vtke  organisms  are  often  observed  in 
environments  altered  by  hydrocarbon  pollution  and/or  haline  environments  (Kasai  et  al. 
2002).  The  Alcanivorax-Vtke  sequences  recovered  from  the  aquifer  are  most  closely 
related  to  genes  recovered  from  marine  sources  such  as  deep-sea  sediments  and  the  water 
column  of  the  North  Sea.  Alteromonad-Vtke  organisms  are  commonly  found  in  coastal 
and  open  ocean  waters,  often  in  the  deep-sea  (Mikhailov  and  Ivanova  1995).  These 
organisms  are  generally  aerobic  heterotrophs.  Sequences  closely  related  to  the  genes 
recovered  from  the  aquifer  have  been  obtained  from  deep-sea  sediments  (Reed  et  al. 
2002),  groundwater  (Gihring  et  al.  2006),  and  surface  seawater  (Bano  and  Hollibaugh 
2002)  environments.  Xanthomonas-ltke  organisms,  often  associated  with  plant 
degradation  and  disease,  are  commonly  in  soil  environments. 

The  alpha-Proteobacteria  are  represented  in  the  aquifer  by  Erythrobacter-, 
Mesorhizobium-  and  Rhodobartcr-hke  16S  rRNA  genes.  Erythrobacter-hke  sequences 
were  recovered  from  the  shallow  and  deep  salinity  zones  in  the  core.  Erythrobacteria  are 
common  photohetertrophs  in  nutrient-rich  coastal  waters  (Koblrzek  et  al.  2003).  The 
isolated  representatives  from  this  group  cannot  use  nitrate  as  a  nitrogen  source,  but 
instead  require  ammonium,  urea  or  amino  acids  (Koblrzek  et  al.  2003). 

Other  alpha-Proteobacteria  recovered  from  the  aquifer  are  closely  related  to 
known  denitrifying  organisms.  Mesorhizobium-hke  organisms  are  often  capable  of 
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nitrogen-fixation  and  are  typically  symbiotic  with  leguminous  plants.  However,  cultured 
Mesorhizobium  species  have  also  been  shown  to  denitrify  under  microaerobic  conditions 
(Okada  et  al.  2005;  Monza  et  al.  2006).  Related  sequences  have  also  been  recovered  from 
uranium  contaminated  soils  (Brodie  et  al.  2006),  ferromanganese  deposits  (Northup  et  al. 
2003)  and  aquiculture  biofilters.  Rhodobacter-like  organisms  are  common  in  both 
freshwater  and  marine  habitats.  These  organisms  are  known  for  their  very  plastic 
metabolism  ranging  from  phototrophy  to  lithotrophy,  aerobic  and  anaerobic  respiration, 
including  denitrification  (Schwintner  et  al.  1998;  Liu  et  al.  1999).  The  genes  recovered 
from  the  aquifer  are  most  closely  related  to  various  Stappia  sp.  within  the  Rhodobacter 
order.  Stappia  contains  members  that  are  known  denitrifiers  (Weber  and  King  2007). 
Interestingly,  the  species  with  the  most  closely  related  16S  rRNA  sequences  to  those 
recovered  from  the  aquifer  have  been  shown  to  denitrify  (Pujalte  et  al.  2005). 
Mesorhizobium-  and  Rhodobacter-like.  sequences  were  recovered  throughout  the 
sediment  core  and  represent  potential  denitrifying  organisms  in  the  nitrate-rich 
environment. 

Actinobacteria-like  organisms  are  common  soil  microorganisms  known  for  their 
ability  to  degrade  refractory  organic  matter  such  as  cellulose  and  chitin  but  also  for  the 
production  of  secondary  metabolites  such  as  antibiotics  (Bull  and  Stach  2007). 
Actinobacteria-like  sequences  were  only  recovered  from  the  upper  two  meters. 
Firmicute-like  organisms  are  another  common  soil  microorganism  that  are  classically 
divided  into  the  anaerobic  (Clostridia)  and  aerobic  (Bacilli)  groupings.  All  of  the 
Firmicute-like  sequences  obtained  from  the  aquifer  were  related  to  the  anaerobic 
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Clostridia  class  with  the  most  closely  related  sequences  recovered  from  anaerobic 
wastewater  digesters  (Tang  et  al.  2004),  groundwater  and  uranium  contaminated  soils 
(Brodie  et  al.  2006).  These  sequences  were  recovered  from  three  meters  and  deeper, 
Planctomycetales-like  organisms  include  the  anammox-like  organisms  as  well  as  other 
aquatic  (marine,  brackish  and  freshwater)  organisms.  Organisms  of  this  affiliation  were 
recovered  from  all  depths. 


Discussion 

Nitrogen  inputs  into  groundwater  on  Cape  Cod  are  primarily  derived  from 
wastewater  (48  %)  with  lesser  contributions  from  atmospheric  deposition  (29  %)  and 
fertilizer  (16  %)  (Valiela  et  al.  1997;  Valiela  et  al.  2000).  Nitrate  from  the  freshwater 
portion  of  the  aquifer  constitutes  the  bulk  of  the  DIN  transport  to  the  head  of  the  bay  and 
has  been  ealeulated  to  supply  6.2  kg  N  day'^  (80  %  of  total)  to  the  head  of  Waquoit  Bay 
(Talbot  et  al.  2003).  Model  results  suggest  limited  nitrate  removal  via  denitrification  is 
likely  to  oeeur  within  the  oxidized,  nitrate  rieh  freshwater  plume  (Spiteri  et  al.  2008). 
However,  below  this  freshwater  plume,  ehemical  data  from  a  similar  section  through  the 
groundwater  aquifer  show  an  increase  in  8^^Nno3  with  deereasing  nitrate  concentrations 
(Kroeger  and  Charette  2008).  Likewise,  a  similar  pattern  of  increase  8^^Nnh4  with 
deereasing  ammonium  concentration  is  also  observed  within  an  ammonium  plume  (>250 
pM  NH4'^)  of  the  freshwater  pareel  (Kroeger  and  Charette  2008).  These  authors  interpret 
these  data  to  refleet  a  loss  of  DIN  (nitrate  and  ammonia)  through  a  fraetionating 
biogeochemical  processes,  presumably  through  coupled  nitrification-denitrification  or 
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anammox  (Kroeger  and  Charette  2008).  Unfortunately,  to  date  there  is  not  enough 
evidenee  in  regard  to  the  isotopie  signatures  or  genetie  markers  to  adequately  examine 
the  role  of  anammox  in  the  eurrent  study.  The  goal  of  the  present  study  is  to  examine, 
using  isotopie  and  moleeular  teehniques,  where  denitrifieation  oeeurs  in  the  subterranean 
estuary  and  whieh  groups  of  organisms  might  be  involved  in  this  proeess. 

Isotopic  evidence  for  denitrification 

The  fresh  groundwater  pareel  eontained  low  8'^Nno3  and  8'^Ono3  values  of  5%c 
and  3%o,  respeetively  (Figure  3 A,  B).  These  values  are  similar  to  average  seawater  (Wada 
et  al.  1975;  Sigman  and  Caseiotti  2001;  Caseiotti  et  al.  2002)  and  do  not  ehange  as  the 
transeet  moves  toward  the  bay.  The  low  8'^Nno3  values  are  eonsistent  with  the  proposed 
anthropogenie  souree  of  nitrate  to  this  pareel.  Ammonium  from  wastewater  has  an 
expeeted  range  of  5%o  to  9%c  (Aravena  et  al.  1993).  Although  nitrifieation  has  a 
signifieant  fraetionation  faetor  (Mariotti  et  al.  1981;  Yoshida  1988;  Caseiotti  et  al.  2003), 
the  8^^Nno3  values  would  be  expeeted  to  approaeh  the  of  the  ammonium  souree 
when  the  ammonium  pool  is  eompletely  eonsumed  within  the  freshwater  plume  at  the 
time  of  our  sampling. 

The  low  values  of  8^^0no3  in  the  fresh  water  plume  suggest  little  alteration  of  the 
nitrate  pool  sinee  its  formation  via  nitrifieation.  Nitrifieation  adds  oxygen  atoms  from 
both  water  and  dissolved  oxygen  in  approximately  a  2: 1  ratio.  Exehange  of  nitrite  with 
water  during  nitrifieation  may  lead  to  a  greater  dependenee  on  water  (Buehwald  and 
Caseiotti  2010;  Caseiotti  et  al.  2010).  Due  to  this  dependeney  on  oxygen  atoms  from 
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water,  information  about  the  8^^0h2o  can  provide  constraints  on  whether  the  nitrate  is  of 
local  origin  or  derived  from  an  external  source.  8^^0h2o  variations  between  -8%o  and  -2%o 
were  well  correlated  (r  =  0.71)  with  salinity  (Figure  3D),  suggesting  mixing  between 
coherent  fresh  and  salt  water  sources.  The  depleted  values  for  8^^0h2o  within  the  core  of 
the  freshwater  plume  are  consistent  with  a  meteoric  source  (-8%o  to  -6%o)  (Kendall  and 
Coplen  2001),  and  are  similar  to  reported  d  Oh2o  values  for  nearby  river  water,  while  the 
high-salinity  end-member  (ca.  80%  seawater)  approached  the  average  d^^OH2o  for 
seawater  (0%o).  NO3'  produced  via  nitrification  with  a  source  of  water  between  -8%o  and  - 
2%o  is  expected  to  have  a  d^^ONos  value  between  -2.3%o  and  -i-0.2%o(Buchwald  and 
Casciotti  2010;  Casciotti  et  al.  2010).  This  range  is  largely  consistent  with  the  observed 
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d  Ono3  values  in  the  core  of  the  fresh  water  plume,  although  the  slightly  higher  d  Ono3 
values  observed  in  the  aquifer  could  be  the  result  of  mixing  with  nitrate  of  higher  d  Ono3 
or  partial  consumption  of  NO3'  within  the  plume. 

Consumption  of  NO3'  would  be  expected  to  leave  the  residual  NO3'  enriched  in 
both  O  and  N  (Mariotti  et  al  1981).  This  is  most  clearly  observed  in  moving  up 
through  the  sediment  into  the  shallow  salt  zone.  In  this  region  the  concentration  of  nitrate 
dropped  (Figure  2)  while  the  values  of  both  8^^Nno3  and  8^^0no3  increased  (Figure  3). 

The  gradients  in  the  isotopic  values  were  greatest  in  the  uppermost  meter  of  the  aquifer. 
Over  this  distance,  8'^Nno3  values  increased  from  +\0%c  to  maximum  values  of  +10%c 
near  the  surface  and  8^^0no3  values  increase  from  -i-5%o  to  +I07cc  (Figure  3).  The  increase 
in  8^^0no3  was  less  than  the  increase  in  8^^Nno3,  resulting  in  a  slope  for  8'^Ono3  vs. 
8'^Nno3  of  0.56  (Figure  4). 
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Alternative  explanations  for  the  increases  in  both  8^^Nno3  and  8^^0no3  towards 
the  surface  include:  1)  addition  of  nitrate  with  high  and  values  and  2)  isotopic 
fractionation  during  nitrate  consumption  within  the  aquifer.  A  detailed  study  of  the 
groundwater  in  a  different  section  of  Waquoit  Bay  revealed  a  shallow  nitrate  plume  with 
low  8^^Nno3  (+2.1%o),  which  the  authors  interpreted  as  nitrate  from  atmospheric 
deposition  (Kroeger  and  Charette  2008).  Atmospheric  deposition  tends  to  have  low 
8'^Nno3  (-4%o  to  +6%o)  (Spoelstra  et  al.  2001;  Burns  and  Kendall  2002)  and  high  8'^Ono3 
(+25%o  to  +75%o)  (Durka  et  al.  1994;  Spoelstra  et  al.  2001;  Burns  and  Kendall  2002) 
values.  This  could  explain  increases  in  8^^0no3,  but  is  not  consistent  with  an  increase  in 
d^^N  (up  to  +20%o).  Addition  of  NO3'  from  atmospheric  deposition  is  expected  to  result  in 
greater  increases  in  8^^0no3  than  in  8'^Nno3,  which  is  not  what  we  observe.  The 
enrichments  in  8^^0no3  and  8^^Nno3  suggest  that  the  bulk  of  this  nitrate  in  the  surface  salt 
zone  at  this  location  is  not  likely  to  arise  from  atmospheric  deposition,  unless  significant 
microbial  reprocessing  has  taken  place  since  deposition  (Barnes  et  al.  2008).  Another 
potential  source  of  nitrate  to  this  zone  could  be  through  infiltrating  seawater  driven  by 
wave  or  tidal  mixing.  However,  nitrate  introduced  through  this  path  would  be  expected  to 
be  isotopically  similar  to  bay  water,  which  has  lower  values  of  both  8'^Nno3  (+6.3  + 
0.06%o)  and  8^^0no3  (+3.8  +  0.05%o)  than  the  surface  aquifer  (Figure  3). 

The  increasing  values  of  both  8'^Nno3  and  8^^0no3  values  co-occuring  with 
decreasing  nitrate  concentrations  are  consistent  with  NO3'  removal  by  an  isotopically 
fractionating  process.  The  data  suggest  that  NO3'  in  the  surface  pore  waters  would  have 
undergone  partial  consumption  to  reach  the  high  8^^Nno3  and  8'^Ono3  values  observed  in 
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the  upper  meter  of  sediment.  The  ratio  of  8^^0no3  and  8^^Nno3  increase  within  this  zone 
(Figure  4)  was  also  consistent  with  reported  values  for  denitrification  in  groundwater 
(slope=0.47-0.77) ,  supporting  the  hypothesis  of  denitrification  occurring  in  the  upper 
sediments  of  the  aquifer.  Furthermore,  we  can  calculate  the  fractionation  factor  (s)  that 
would  result  from  the  consumption  of  nitrate  in  this  zone  from  the  freshwater  peak  (332 
pM  NOa')  to  the  surface  (46  pM  NO3')  using  a  Rayleigh  type  model: 

d„=d„-£ln(/)  (1) 

where  8n  is  the  isotope  value  at  time  n,  80  is  the  initial  isotope  value  and/ is  the  fraction 
of  nitrate  remaining.  This  calculation  provides  an  estimate  of  +7.6%o  for  s,  well  within 
the  reported  range  of  values  from  known  denitrifiers  (+5  to  +35%o).  The  natural 
abundance  stable  isotopes  thus  provide  compelling  evidence  for  active  denitrification. 
However,  additional  investigations  are  needed  to  constrain  the  instantaneous  rates. 

Below  the  fresh  groundwater  plume,  the  nitrate  concentration  drops  rapidly 
approaching  the  deep  salt  wedge,  resulting  in  a  steep  concentration  gradient. 
Unfortunately,  samples  along  this  gradient  with  sufficient  nitrate  concentrations  to 
measure  the  isotopic  values  were  limited.  However,  based  on  the  nitrate/salinity 
relationship  and  isotopic  fractionation  reported  by  others  ,  the  nitrate  pool  appears  to  mix 
non-conservatively  within  this  lower  portion  of  the  aquifer.  This  suggests  nitrate  is 
removed  between  the  freshwater  and  deep  salt  wedge. 

We  interpret  these  data  collectively  as  representing  a  nitrate  pool  produced 
somewhere  upstream  of  the  transect  zone  that  remains  largely  unchanged  at  its  core 
during  its  transit  through  the  area  covered  by  our  transect.  At  both  the  surface  salt  zone 
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and  deep  salt  wedge  transitions,  nitrate  removal  appears  to  occur  under  lowered  O2 
conditions.  In  the  surface  salt  zone  at  the  landward  end  of  the  transect,  both  8^^Nno3  and 
8'*Ono3  values  increase  with  decreasing  nitrate  concentrations,  which  is  consistent  with 
denitrification.  The  question  now  becomes  which  organisms  might  be  responsible  for  this 
consumption  of  nitrate. 

Molecular  evidence  for  denitrification 

Modern  molecular  techniques  can  be  used  to  shed  light  on  which  organisms  and 
perhaps  which  metabolisms  could  be  functioning  at  different  depths  within  the  aquifer 
sediments.  In  a  general  overview  the  bacterial  16S  rRNA  gene  clone  libraries  suggest  that 
the  Waquoit  Bay  groundwater  system  harbors  organisms  that  are  likely  adapted  to 
nutrient- amended  environments.  Many  of  the  sequences  recovered  in  this  study  such  as 
the  gamma-Proteobacteria,  Planctomycetales-  and  Firmicute-Wke  organisms  indicate 
adaptation  to  anaerobic  sediments  and  soils.  Given  that  many  of  these  sequences  have 
been  recovered  from  wastewater  or  sludge  bioreactors  it  may  be  that  these  organisms  are 
also  adapted  to  high  nitrate,  low  oxygen  environments.  The  study  site  at  Waquoit  Bay 
contains  steep  gradients  in  nutrient  and  oxygen  concentrations  that  may  set  boundaries  on 
the  distribution  of  organisms  within  this  environment.  To  evaluate  this,  we  combined 
phylogenetic  analysis  of  the  16S  rRNA,  nirK  and  nirS  genes  as  well  as  the  quantitative 
depth  distribution  of  the  two  functional  genes  through  the  sediment  column. 
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The  upper  meter 

The  upper  meter  of  the  sediment  column  represents  an  area  that  likely  harbors 
denitrifying  organisms.  Analysis  of  the  16S  rRNA  genes  recovered  from  the  upper  meter 
shows  that  it  is  dominated  by  alpha-and  gamma-Proteobacteria.  Many  of  these  organisms 
are  related  to  aerobic  heterotrophs  from  environments  such  as  marine  and  estuary 
sediments,  marine  water  column  and  groundwater  environments  among  others.  Notably, 
Mesorhizobium-like,  and  Stappia-like  sequences  were  recovered  from  this  depth  range  but 
not  in  the  deeper  regions  of  the  sediment  column.  The  Stappia-like  sequences  may 
represent  a  very  versatile  metabolic  lifestyle  that  potentially  includes  denitrification 
(Kroeger  and  Charette  2008).  Mesorhizobium-like  organisms  include  known  nitrogen- 
fixers  as  well  as  some  aero  tolerant  denitrifying  strains  (Weber  and  King  2007). 
Mesorhizobium-like  organisms  have  been  implicated  in  the  nirK-mediated  denitrification 
in  rhizosphere  environments  (Okada  et  al.  2005),  wastewater  (Yoshie  et  al.  2004)  and 
soils  (Sharma  et  al.  2005).  This  finding  is  congruent  with  the  abundance  data  that  shows 
nirK  genes  are  more  abundant  than  nirS  in  the  upper  meter  of  the  sediment  column.  The 
nirK  sequences  recovered  from  this  region  of  the  aquifer,  as  well  as  from  other  sections, 
are  also  most  closely  related  to  terrestrial  organisms  such  as  Mesorhizobium-, 
Rhodobacter-,  and  Hyphomicrobium-like  organisms.  Consistent  with  low  abundance 
levels,  only  one  sequence  type  of  nirS  was  retrieved  from  the  upper  meter  of  the  sediment 
column  most  closely  related  to  the  beta-Proteobacterium  Ralstonia  sp.  Taken  together 
with  the  chemical  data  above  the  surface  salt  zone  appears  to  be  a  zone  of  active 
denitrification  carried  out  by  a  soil  or  sediment  type  microbial  community. 
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The  Freshwater  Plume 


The  freshwater  plume  contains  elevated  nitrate  and  oxygen  concentrations  and 
depleted  isotopic  values  for  both  8^^Nno3  and  8^^0no3-  As  described  above,  the  chemical 
signatures  do  not  indicate  appreciable  denitrification  in  this  zone.  Many  of  the  16S  rRNA 
genes  recovered  from  this  depth  are  also  present  in  the  1  m  and  the  3  m  samples. 
Sequences  related  to  various  alpha-  and  gamma- Proteobacteria  including  Alteromonas- 
like  organisms  were  most  common.  Alteromanades  are  common  marine  heterotrophs, 
ranging  from  coastal  waters  to  open  ocean  and  deep-sea  habitats  that  have  also  been 
implicated  in  denitrification  (Bremer  et  al.  2007).  At  this  depth  nirS  gene  abundance  has 
increased  to  a  level  equal  to  nirK  gene  abundance.  Without  a  chemical  signature  for 
denitrification,  and  given  that  these  are  gene  abundance  values  and  not  transcript 
(mRNA)  abundance,  we  interpret  the  presence  of  these  genes  at  this  depth  as  indicating 
the  presence  of  potential  denitrifiers  that  are  utilizing  another  metabolic  pathway.  Within 
the  freshwater  plume,  the  high  oxygen  concentrations  are  expected  to  inhibit  traditional 
denitrifying  organisms.  Studies  have  shown  denitrification  can  be  inhibited  by  oxygen 
concentrations  above  10  pM  for  environmental  studies  (Enger  et  al.  1987),  though  some 
cultured  isolates  have  displayed  greater  aerotolerance.  For  example,  cultures  of 
Pseudomonas  stuzeri,  a  well-studied  denitrifying  bacterium,  has  been  shown  to  denitrify 
in  up  to  150  pM  oxygen  (Nielsen  et  al.  1990).  However,  our  isotopic  and  chemical  data 
suggest  that  the  nitrate  in  the  well-oxidized  fresh  groundwater  is  largely  unaltered  by 
denitrifcation,  regardless  of  the  abundance  of  denitrification  genes  at  this  depth. 
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Deep  Salt  Wedge  Transition 

The  transition  between  the  freshwater  plume  and  deep  salt  wedge  results  in  a 
steep  gradient  in  nitrate  and  oxygen  concentrations.  Recent  reports  suggest  that 
denitrification  does  occur  within  this  zone  (Korner  and  Zumft  1989).  Interestingly,  there 
is  a  transition  observed  in  the  16S  rRNA  gene  libraries  at  this  depth.  The  sequences 
recovered  from  three  meters  and  below  are  most  commonly  represented  by  sequences 
related  to  Erythrobacter  (also  seen  in  upper  meter),  Clostridia  and  Planctomycetales. 
Clostridia-  and  Planctomycetales-like  sequences  were  not  recovered  in  the  surface  salt 
zone.  Similar  sequences  have  been  recovered  from  coastal  marine  waters  (Erythrobacter- 
and  Planctomycetes-  but  not  Clostridia-Yike),  anaerobic  soils  and  marine  sediments 
(Kroeger  and  Charette  2008),  groundwater,  contaminated  soils  (Tang  et  al.  2004)  and 
wastewater  reactors  {Clostridia-  and  Planctomycetales-  but  not  Erythrobacter-like). 
Plancotmycetales  are  organisms  of  special  interest  in  nitrogen  cycling  because  some 
members  of  this  phylum  have  also  been  implicated  in  anammox,  although  the  sequences 
recovered  from  the  aquifer  are  more  closely  related  to  sequences  obtained  from  marine 
sediments  and  soils.  Planctomycetes  are  common  constituents  of  soil  and  marine 
sediment  microbial  communities  and  some  have  very  broad  metabolic  capabilities  as  well 
as  complex  life  cycles. 

Nitrite  reductase  genes  from  the  deep  salt  wedge  are  different  than  found  in  the 
upper  meter.  NirS  increases  by  3  orders  of  magnitude  compared  to  the  surface  salt  zone, 
becoming  dominant  over  the  nirK  gene  (Figure  5).  Recovered  nirS  genes  are  most  closely 
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related  to  Marinobacter-,  Pseudomonas-  and  Paracoccus-Vike,  organisms.  The  presence 
of  these  genes  and  organisms  suggests  that  denitrification  is  possible  as  suggested  by 
earlier  chemical  studies  (Brodie  et  al.  2006).  The  isotopic  and  chemical  signature  of 
groundwater  denitrification  is  not  observed  in  the  present  study,  most  likely  because  of 
the  steep  nitrate  gradient  at  depth,  which  prevented  our  natural  abundance  isotopic 
measurements  due  to  detection  limitations  (0.5  mM  NO3'  required). 

Another  pathway  for  nitrogen  removal  is  anammox.  Anammox  traditionally 
combines  ammonia  and  nitrite  to  form  dinitrogen,  in  contrast  to  classical  denitrification 
that  uses  only  the  oxidized  nitrogen  species.  Previous  (Kroeger  and  Charette, 
unpublished)  and  concurrent  incubation  studies  (Rogers  and  Casciotti,  unpublished)  using 
isotopic  tracers  did  not  detect  significant  rates  of  anammox  in  the  subterranean  estuary, 
although  both  denitrification  (Kroeger  and  Charette,  unpublished)  and  nitrification 
(Chapter  4)  were  observed.  Interestingly,  ladderane  lipids  that  are  thought  to  be  specific 
to  anammox  were  recovered  from  the  aquifer  sediments  (Saenz,  unpublished).  Recent 
reports  have  suggested  that  anammox  organisms  may  contain  a  specific  type  of  nirS  gene 
that  could  be  quantified  separately  from  the  classical  denitrifiers  in  the  environment 
(Kroeger  and  Charette  2008).  This  was  not  pursued  in  the  present  study  because  the  nirS 
sequences  from  anammox-like  organism  were  not  available  at  the  time  this  study  was 
conducted.  Given  the  presence  of  lipid  markers  and  the  reported  non-conservative 
behavior  of  ammonium  in  the  mixing  zone  between  the  freshwater  plume  and  deep  salt 
wedge  (Lam  et  al.  2009)  anammox  processes  are  expected  to  be  active  in  the  aquifer. 
There  are  several  reasons  the  present  study  may  have  missed  anammox  in  the  aquifer. 
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The  steep  gradient  in  nitrate  within  this  deep  mixing  zone  limits  our  analysis  in  terms  of 
sample  number  for  both  concentration  and  natural  abundance  isotopes  and  we  do  not 
observe  a  nitrite  pool.  Even  if  natural  abundance  isotope  values  were  available,  it  is  not 
known  what  signal  is  produced  by  biological  or  manganese-oxide  mediated  anammox 
processes.  In  our  tracer  level  studies  handling  of  the  samples,  especially  exposure  to 
elevated  oxygen  concentrations,  may  temporarily  or  permanently  inhibit  the  anaerobic 
microorganisms.  The  reaction  rate  of  bacterially  mediated  anammox  could  also  be  slow 
enough  that  it  cannot  be  observed  over  the  time  course  of  the  incubation  experiment. 
Finally,  the  requirement  of  both  reduced  ammonium  and  oxidized  nitrite  and  an 
anaerobic  zone  for  the  anammox  process  is  likely  to  also  require  the  steep  gradients  and 
associations  with  aerobic  organisms  for  the  production  of  nitrite.  The  incubations  used 
here  were  designed  to  eliminate  gradients  that  may  obscure  nitrification  or  denitrification 
processes.  Clearly,  this  area  of  research  will  benefit  from  future  targeted  studies  and 
technique  development  in  both  isotopic  methods  and  the  development  of  molecular 
markers. 

Implications 

The  results  of  the  chemical  interrogation  of  the  study  site  reveal  three  distinct 
zones:  1)  a  shallow  salt  zone  with  low  pore  water  O2  and  nitrate  concentrations,  2)  a 
fresh,  oxidized  nitrate  plume  (up  to  350pM  nitrate)  and  3)  a  deep  intruding  salt  wedge 
that  is  rich  in  ammonium,  but  lacking  oxygen  and  nitrate.  While  each  of  these 
environments  harbors  communities  that  are  capable  of  cycling  and  removing  nitrogen 
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before  the  aquifer  waters  are  discharged  into  coastal  waters,  our  data  suggest  that  the 
fringes,  rather  than  the  central  plume,  may  be  the  most  active  denitrification  zones.  The 
surface  salt  zone  displays  a  pattern  of  decreasing  nitrate  and  elevated  d^^NNos  and 
8'*Ono3  that  are  consistent  with  previously  observed  groundwater  denitrification. 
Moreover,  the  organisms  present  are  representative  of  soil/sediment  type  organisms 
including  many  potential  denitrifiers.  In  contrast,  the  community  as  assayed  by  the  nir 
genes  reveals  a  low  diversity  of  denitrifiers,  with  nirK  playing  an  important  role. 
Removal  of  nitrogen  from  the  transition  between  the  freshwater  plume  and  deep  salt 
wedge  is  also  likely  where  the  community  shifts  to  a  nirS  dominated  denitrifier 
community.  This  microbial  community  appears  more  closely  related  to  marine  or 
reducing  sediment/wastewater  type  organisms.  The  elevated  oxygen  levels  within  the 
freshwater  nitrate  plume  most  likely  inhibit  denitrification,  although  available  carbon  or 
electron  donors  may  also  be  limiting.  Of  course,  the  processes  that  occur  where  the 
freshwater  plume  discharges  into  a  saline  environment  may  be  similar  to  the  salinity 
transition  zones  observed  in  this  study,  and  may  yet  have  the  potential  to  attenuate  NO3' 
concentrations  prior  to  discharge. 
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Table  1.  qPCR  primers  and  thermal  programs. 


Target 

Gene 

Primer  Names 

And 

Sequences 

Thermal  Program 

Positive 

Control 

Reference 

nirS 

nirSchSAF  (lOOOnM) 

GTSAACGTSAAGGARACS 

GG 

nirSR3cd  (lOOOnM) 

GASTTCGGRTGSGTCTTG 

A 

50°C  for  120sec 

95°C  for  600sec 

40  cycles  of 

95°C  for  60sec 

50°C  for  60sec 

60°C  for  60sec 

80°C  for  15sec  (read) 

Psuedomonas 

stuzeri 

(Throback  et  al. 

2004;  Kroeger  and 
Charette  2008) 

nirK 

nirK876  (SOOnM)- 
ATYGGCGGVAYGGCGA 
nirK1040  (SOOnM)- 
GCCTCGATCAGRTTRTGG 
TT 

50°C  for  120sec 

95°C  for  900sec 

6  cycles  of 

95°C  for  15  sec 

63-58°C  for  60sec  (- 
1°C /cycle) 

72°C  for  30sec 

80°C  for  15sec  (read) 

40  cycles  of 

95°C  for  15sec 

58°C  for  30sec 

72°C  for  30sec 

80°C  for  15sec  (read) 

Alcaligenes 
faecalis  ATCC 
8750 

(Geets  et  al.  2007) 

Bacterial 

16S  rRNA 

8F  (400nM)- 

AGAGTTTGATCCTGGCTC 

AG 

1492R  (400nM)- 

GGTTACCTTGTTACGACT 

T 

94°C  for  120sec. 

30  cycles  of 

94°C  for  30sec 

47°C  for  60sec 

72°C  for  60sec 

Final  extension  of 

72°C  for  600sec 

Escherichia 

coli. 

Pseudomonas 

stutzeri 
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A  B  C  D 


Figure  1.  Sampling  schematic  depicting  the  relative  position  of  the  piezometer  locations 
(A,  B,  C,  D)  in  relation  to  the  aquifer  structure.  The  aquifer  is  composed  of  three  distinct 
water  parcels,  (1)  a  surface  salt  zone  (stippled),  (2)  a  fresh  groundwater  plume  (white) 
and  (3)  a  deep  salt  wedge  (large  stipples).  The  mixing  zone  between  the  freshwater  and 
deep  salt  wedge  is  noted  with  a  gradient  in  stippling.  Pulse-coring  was  performed  at  site 
B  during  June  2008. 
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Figure  2.  Chemical  profiles  along  a  transect  normal  to  the  shore  line  at  the  head  of 
Waquoit  Bay,  Cape  Cod,  MA.  Transect  moves  from  on-shore  to  off,  left  to  right  with 
positive  distance  above  the  mean  high  tide.  Black  circles  represent  sampling  points  along 
each  of  the  four  piezometer  locations.  Sediment  coring  was  done  at  station  B.  A)  Salinity 
is  most  pronounced  in  the  deep  salt  wedge  and  shoals  moving  off-shore.  B)  Dissolved 
oxygen  (pM)  peaks  in  the  fresh  groundwater  and  shoals  off-shore.  C)  Nitrate  (pM)  peaks 
in  the  fresh  groundwater  but  is  also  present  in  the  surface  salt  zone.  Nitrate  is  absent  from 
the  deep  salt  wedge.  D)  Ammonium  is  not  observed  in  the  fresh  groundwater  parcel  but 
routinely  found  in  the  deep  salt  wedge  likely  due  to  remineralization  of  organic  carbon  or 
cation  exchange. 
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Figure  3.  Natural  abundance  isotope  profiles  along  a  transect  normal  to  the  shore  line  at 
the  head  of  Waquoit  Bay,  Cape  Cod,  MA.  Transect  moves  from  on-shore  to  off,  left  to 
right  with  positive  distance  above  the  mean  high  tide.  Black  circles  represent  sampling 
points  along  each  of  the  four  piezometer  locations.  A)  The  8^^Nno3  across  the  transect 
reveals  low  values  similar  to  bay  water  in  the  core  of  the  freshwater  plume  and  increasing 
values  in  the  surface  pore  water  of  the  two  most  shoreward  piezometers.  B)  8^^0no3 
reveals  a  similar  distribution  as  d^^NN03-  C)  8^^0h2o  is  most  depleted  in  the  core  of  the 
fresh  groundwater  plume.  D)  Plotting  8^^0h2o  against  salinity  reveals  a  steady  mixing 
line  at  all  but  the  lowest  values. 
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Delta  '^N^o3 

Figure  4.  The  relationship  between  8'^Nno3  and  8^^0no3  in  the  upper  meter  of 
groundwater  at  Waquoit  bay  is  consistent  with  previous  reports  of  groundwater 
denitrification. 
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Figure  5.  The  nitrite  reductase  occurs  in  two  types,  a  copper  based  nirK  and  a 
cytochrome  based  nirS.  Abundance  of  nirS  genes  (grey  circles)  increase  with  depth, 
peaking  at  3.8  m.  NirK  gene  abundance  (black  squares)  varies  less  than  an  order  of 
magnitude  throughout  the  core.  Error  bars  represent  one  standard  deviation  based  on 
triplicate  analysis. 
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A 


Pseudomonas  stutzeri,  AM230893 
Pseudomonas  chloritidismutans.  AJ664572 
Alcaligenes  faecatis.  AJ224913 

uncultured  bacterium.  AY195909 
“  urKultured  bacterium,  AY19S906 
uncultured  bacterium,  AY195904 
uncultured  bactenum,  AY195902 
uncultured  bacterium,  AY195900 
uncultured  bacterium.  AY195905 
Azoarcus  evartsli.  AY07B269 
Azoarcus  loluvorans.  AY07d270 
Thauera  sefenatis.  AY078264 
Comamonas  denliriricans.  D0865926 
Comamonas  sp  R-28452,  AM230922 
Acidovorax  sp  R-2S212,  AM230905 
—  Comamonas  sp  R-2M66,  AM230e97 
Comamonas  sp  R-24447,  AM230695 
Comarrwnas  sp  R-26227,  AM230921 
Pseudomonas  sp  R-25061 ,  AM230904 
Comamonas  sp.  R-25014.  AM230898 
bacterium  K,  DQ106865 

Azoarcus  totulyticus,  AY078271 
uncultured  bacterium.  AB456900 
WBnirS42 
WBnlrSSS 

uncultured  bacterium,  AB456673 
uncultured  bactenum.  AB456871 
uncultured  bacterium,  AB456929 

uncultured  bacterium,  AB456885 
•  WBnirS43 
1  WBnirSSa 

'  uncultured  organism,  AM419579 
uncultured  organism,  AM419556 

j  uncultured  organism,  AM419630 
I  I  uncultured  organism,  AM419566 
I  J  uncultured  organism,  AM4195S3 
L_r^  uncultured  organism,  AM419562 
uncultured  organism.  AM419557 

WBnirS39 

WBnirSS 

WBnlrS49 

uncultured  organism.  AM419565 
WBnlrS38 
WBnirS54 

uncultured  bactenum,  AB456942 
WBnlrSH 
WBnirSlO 
WBnirS12 

ur>cultured  bactenum.  D0676040 


_l  uncultured  bacterium.  DQ676190 

^  WBnirS31 


I7wb8 


Azoarcus  lolulytkxis.  AY078272 
I  uncultured  bacterium.  EF556472 
'  WBnirS21 

uncultured  bacterium.  EF556416 
I  uncultured  bacterium,  DO676057 
'  WBnirS63 

uncultured  bactenum.  DQ676127 
I  WBnlrS32 
I  WBnirS57 

'  uncultured  bacterium.  EF55848S 

uncultured  bacterium,  EF558391 


— I  7  WB  sequences 


-c 


uncultured  bacterium,  EF558379 
WBnirS22 

uncultured  bacterium,  EF55S407 
uncultured  bacterium.  EFS58378 
1  WBnirS27 

P  uncultured  organism,  DQ767847 
L  WBnirS29 

~  uncultured  organism,  DQ767853 
uncultured  bactenum.  EF558532 


— I  5  WB  sequences 


WBnirSi 
—  WBnlrS68 


_|  WBnirS28 
^  WBnlrS67 


PseudorT}onas  aeruginosa.  FJ976652 

Pseudo/TTonas  fluorescens.  AF197466 
Marinobacter  sp.  NBC35.  AJ626837 
Marinobacter  sp  R53.  AJ626634 
Marinobacter  sp  R171,  AJ626836 
uncultured  bacterium.  EU325639 


j  WBnirS41 
1  WBnirSS2 

'  uncultured  bacterium,  EU32S642 

H  12  /  WB  strains  10  sequences 
WB  sequences 

uncultured  bacterium.  AJ440476 

uncultured  bacterium,  00085422 
uncultured  bacterium.  DQ085424 

uncultured  organism,  00094218 
uncultured  organism,  DQ094195 


WBnirS23 

Marinobacter  sp 


Paracoccus  pantotrophus,  AJ401462 

Ralstonia  eutropha  H16.  X91394 
Brachymonas  denitrificans.  D0865925 
—  Paracoccus  denitrificans.  U7S413 

Paracoccvsdenitrifican$PD122.  U05002 
WBnlrS45 


S64.  A  J 626835 


0.10 
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17 


Figure  6.  A)  Phylogenetic  relationship  of  recovered  nirS  (A)  and  nirK  (B)  sequences  to 
sequences  present  in  Genbank.  Sequences  were  imported  in  to  the  ARB  phylogenetic 
software  and  aligned  using  the  ARB  aligner  followed  by  manual  editing.  Trees  were 
constructed  using  the  neighbor  joining  function  with  a  conservative  mask  of  552  bp 
(nirS)  and  427  bp  (nirK).  Sequences  recovered  from  this  study  are  in  bold  and  known 
cultures  are  italicized. 
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Abstract 


Submarine  groundwater  discharge  can  be  an  important  source  of  nitrogen  to 
coastal  waters.  However,  little  is  known  about  the  microbial  nitrogen  cycle  in 
groundwater  settings.  Previous  work  has  shown  the  presence  three  populations  of 
nitrifiers  in  the  Waquoit  Bay  groundwater  system,  two  different  archaeal  populations,  one 
in  the  upper  meter  and  one  at  the  transition  between  the  freshwater  and  deep  salt  wedge 
as  well  as  a  bacterial  population  within  the  freshwater  plume.  The  aim  of  the  current 
study  is  to  estimate  the  potential  nitrification  rates  at  three  depths  in  the  groundwater 
system  using  a  ^^N-labeled  ammonium  technique.  Allylthiourea  (ATU),  an  inhibitor  for 
bacterial  nitrification,  was  used  to  differentiate  between  archaeal  and  bacterial 
nitrification  activity.  The  expression  of  bacterial  and  archaeal  ammonia  monooxygenase 
genes  (amoA)  was  also  followed  over  the  time  course  of  the  incubations.  Potential 
nitrification  rates  in  the  surface  salt  zone  and  freshwater  plume  were  similar  to  estuarine 
sediments.  Incubations  from  the  surface  salt  zone  contained  archaeal  amoA  expression 
that  was  not  suppressed  in  the  presence  of  ATU.  This  contrasts  with  the  incubations  from 
the  freshwater  plume  where  bacterial  amoA  was  expressed  but  showed  sensitivity  to 
ATU.  The  deep  salt  wedge  contained  a  population  of  putative  archaeal  ammonia- 
oxidizers,  as  determined  from  the  presence  of  archaeal  amoA  genes;  however,  evidence 
of  active  nitrification  was  lacking  in  these  deep  sediments. 
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Introduction 


Submarine  groundwater  discharge  (SGD),  the  transport  of  terrestrial  groundwater 
and  recirculated  seawater  through  sediments  and  into  coastal  waters,  is  increasingly 
recognized  as  an  important  route  for  chemical  discharge  to  the  coastal  ocean  (Moore 
1999).  Though  the  global  freshwater  flux  from  SGD  is  small  (ca.  6%,  Burnett  et  al.  2003) 
compared  to  riverine  water  flux,  the  impact  on  the  local  environment  may  be  important 
due  to  the  large  chemical  load  these  waters  can  carry.  For  example,  with  one-third  of  the 
world’s  population  living  within  100  km  of  the  coast  (Cohen  et  al.  1997)  many 
groundwater  systems  carry  large,  anthropogenically  derived,  nutrient  loads.  This  nutrient 
addition  to  coastal  waters,  known  as  eutrophication,  has  broad  deleterious  effects  on 
coastal  ecosystems  including  increasing  the  frequency  and  duration  of  harmful  algal 
blooms  and  expanding  regions  of  hypoxia  and  anoxia.  Nutrient  addition  also  alters  local 
food  webs,  degrades  sea  grass  beds,  decreases  biodiversity,  and  compromises  fin-fish  and 
shellfish  populations  for  human  consumption  (Valiela  et  al.  1990;  Valiela  et  al.  1992; 
Howarth  et  al.  2002). 

The  biogeochemical  processes  that  occur  within  the  groundwater  system  largely 
control  the  scale  of  nutrient  input,  principally  nitrogen,  and  are  therefore  important  to 
understand.  Unlike  surface  waters,  which  flow  rapidly  over  a  sediment  matrix  in 
communication  with  the  atmosphere,  groundwaters  typically  have  slow  advection  rates 
flowing  through,  rather  than  over,  the  sediment  matrix  and  are  isolated  from  the 
atmosphere.  Thus  the  geochemical  environment  of  groundwaters  can  be  very  different 
than  that  found  in  riverine  settings.  This  contrast  in  physical  environment  can  have  a 
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significant  effect  on  the  biogeochemical  cycling  of  nitrogen  in  groundwaters.  In 
particular,  many  of  the  transformations  nitrogen  can  undergo  are  microbially  controlled 
and  sensitive  to  oxygen  availability,  which  can  vary  widely  in  groundwater  systems. 
Therefore,  the  suite  of  nitrogen  transformations  that  is  possible  in  a  subterranean  estuary 
is  very  broad. 

Oxidation  of  nitrogen  in  the  environment  is  observed  in  the  conversion  of 
ammonium  through  nitrite  to  nitrate.  Specific  bacteria  (ammonia-oxidizing  bacteria, 
AOB)  and  archaea  (ammonia-oxidizing  archaea,  AOA)  catalyze  the  oxidation  of 
ammonia  for  the  conservation  of  energy.  A  separate  group  of  bacteria,  nitrite-oxidizing 
bacteria,  conserve  energy  through  the  conversion  of  nitrite  to  nitrate.  Both  the  AOB  and 
AOA  utilize  a  (putatively  related)  protein  called  ammonia  monooxygenase  to  catalyze  the 
first  step  in  this  pathway,  converting  ammonia  to  the  first  intermediate,  hydroxylamine. 
The  gene  coding  for  one  of  the  subunits  of  this  enzyme,  ammonia  monooxygenase 
subunit  A  (amoA),  has  been  used  as  a  marker  for  the  presence  (Francis  et  al.  2003; 
Leininger  et  al.  2006;  Wuchter  et  al.  2006)  and  activity  of  these  organisms  (Treusch  et  al. 
2005;  Lam  et  al.  2009) 

AOB  were  discovered  over  one  hundred  years  ago  by  Winigradsky  (Winogradsky 
1890)  and  were,  until  recently,  thought  to  be  the  only  organisms  capable  of  carrying  out 
this  portion  of  the  nitrogen  cycle.  In  the  last  decade,  however,  archaea  also  have  been 
demonstrated  to  participate  in  ammonia  oxidation,  first  from  molecular  evidence  from 
various  environments  (Venter  et  al.  2004;  Francis  et  al.  2005;  Treusch  et  al.  2005) 
followed  by  further  genomic  evidence  of  carbon  dioxide  fixation  and  ammonia-oxidation 
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pathways  (Hallam  et  al.  2006)  and  activity  in  cultured  representatives  (Konneke  et  al. 
2005;  Santoro  et  al.  2008).  AOB  and  AO  A  have  been  found  to  co-occur  in  many 
environments,  presumably  competing  for  the  ammonia  substrate.  It  is  still  unclear  what 
environment  factor(s)  may  favor  one  population  of  the  other.  Possible  roles  for  oxygen 
tension  (Santoro  et  al.  2008),  carbon  concentration  (Hallam  et  al.  2006),  salinity  (Francis 
et  al.  2005;  Caffrey  et  al.  2007),  pH  (de  la  Torre  et  al.  2008)  and  differences  in 
ammonium  concentration  (Martens-Habbena  et  al.  2009)  have  all  been  put  forth  as  key 
environmental  factors  controlling  the  balance  between  AOBs  and  AO  As. 

Here  we  combine  molecular  techniques  to  determine  which  genes  are  present  and 
expressed  with  isotope  tracer  incubations  to  build  our  knowledge  of  the  players  and 
processes  active  in  a  groundwater  system  and  how  their  activity  may  alter  the  input  of 
nitrogen  into  the  coastal  waters. 


Methods 

Sample  Collection 

Groundwater  samples  were  collected  on  the  ebbing  tide  using  a  stainless  steal 
piezometer  equipped  with  a  ”Retract-a-tip”  (AMS,  Idaho)  with  a  vertical  resolution  of 
0.15  to  0.30  m.  The  tubing  was  flushed  with  at  least  one  liter  of  pore  water  (more  than  6- 
fold  flushing  of  tubing  volume)  before  measuring  pH,  redox  potential,  dissolved  oxygen, 
salinity  and  temperature  using  a  YSI  600R  multi-probe  (Charette  and  Sholkovitz  2002). 
Oxygen,  pH,  and  redox  sensors  were  calibrated  against  standard  solutions  following  the 
manufacturers  protocols  prior  to  deployment  in  the  field.  Porewater  from  from  target 
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depths,  as  guided  by  previous  molecular  studies  (Chapter  2),  was  collected  in  4  L 
cubitainers  and  capped  while  overflowing  to  minimize  alteration  of  the  in  situ  oxygen 
concentrations. 

After  pore  water  sampling,  sediment  was  collected  near  low-tide  for  nitrification 
rate  measurements  from  within  0.5  m  of  the  piezometer,  using  a  hand-operated  bailer¬ 
boring  auger  (Eijkelkamp,  The  Netherlands).  The  interfaces  between  the  shallow  salt 
zone/fresh  groundwater  (-0.6  m)  and  the  fresh  groundwater/deep  salt  zone  (-4  m)  were 
targeted,  as  well  as  the  center  of  the  groundwater  plume  (2  m). 

Nutrient  measurement 

Samples  were  collected  for  dissolved  inorganic  nitrogen  analyses  (NO3',  NO2’, 
NH4‘^),  filtered  through  a  0.2  pm  pore-size  Sterivex  filter  and  frozen.  Nutrient  samples 
were  frozen  until  analysis.  Nitrite,  nitrate  and  ammonium  were  determined  using  a  Lachet 
QuickChem  8000  continuous  flow  injection  system  at  the  WHOI  Nutrient  Analytical 
Facility. 

^^N-Ammonium  Tracer  Experiment 

In  September  2009  a  pulse-core,  reaching  a  depth  of  4  m,  was  recovered  from 
near  piezometer  B  (Figure  1).  Three  depths  (0.6,  2  and  4  m)  were  targeted  based  on 
previous  observations  of  microbial  populations  within  and  the  prevailing  groundwater 
chemistry  of  the  aquifer  (Chapter  2).  Groundwater  from  each  of  these  depths  was 
collected  using  a  stainless  steel  piezometer,  filtered  (0.2  pm)  and  amended  with  ^^N- 


120 


labeled  ammonium  (^^NH4‘^,  99.8  atom  %,  50  )j,M  final  concentration)  or  '^N-labeled 
ammonium  plus  allylthiourea  (ATU,  50  pM  final  concentration).  Incubations  were  done 
in  triplicate  and  contain  approximately  100  g  sediment  (quantified  post-incubation)  in  a 
160  mL  glass  serum  bottle  with  the  remaining  volume  filled  with  the  amended  water 
from  the  same  depth.  Oxygen  concentrations  in  the  collected  water  were  determined 
using  the  YSI  meter  prior  to  ammonium  addition.  Incubations  were  stored  at  12  °C  and 
shaken  continually  (ca.  100  rpm)  during  the  course  of  the  experiment.  Incubations  were 
destructively  sampled  in  triplicate  at  24  hr  intervals  up  to  96  hr.  Nutrient  concentration 
and  isotope  (NO2  and  NO3'  at  every  time  point,  NH4’^  initial  atom  %)  samples,  taken  in 
separate  20  mL  aliquots  (200  mL  for  initial  NH4’^),  were  filter  through  a  0.2  pm  syringe 
filter  (Millipore)  into  and  HDPE  bottle  and  frozen  (-20  °C)  until  analysis.  Sediment  was 
also  collected  and  frozen  (-80  °C)  for  DNA  (2  g)  and  RNA  (10  g)  extraction  and 
subsequent  molecular  analysis  (see  below). 

Isotopic  Analysis 

of  the  initial  ammonium  pool  was  measured  in  the  ammonia-oxidation 
experiments  using  the  ammonia  diffusion  method  (Sigman  et  al.  1997;  Holmes  et  al. 
1998).  In  this  procedure,  ammonia  is  diffused  from  a  100  mL  sample  and  trapped  on  pre¬ 
combusted  and  acidified  glass  fiber  filter  (GF-D,  Millipore)  that  is  encapsulated  between 
two  Teflon  membranes  (Millipore).  After  14  days  of  diffusion  on  a  shaking  (ca.  100 
rpm),  60  °C  incubator  the  filter  was  dried  and  packed  into  a  tin  cup  and  analyzed  by 
combustion  on  a  PDZ  Europa  ANCA-GSE  elemental  analyzer  interfaced  to  a  PDZ 
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Europa  20-20  isotope  ratio  mass  spectrometer  (Sercon  Ltd.,  Cheshire,  UK)  by  the 
University  of  California  at  Davis,  Stable  Isotope  Facility  (UCDavis-SIF).  Standards  were 
made  from  dilutions  a  400  pM,  99.5-atom  %  NH4'^  (Cambridge  Isotope  Laboratories) 
with  a  400  pM  stock  NH4"^  (Thermo-Fisher)  to  make  additional  stocks  of  75,  50  and  25- 
atom  %.  Standards  (12.5  mL)  were  diluted  with  DI  (87.5  mL)  for  diffusion  and  trapping 
with  the  sample  set.  Standards  and  samples  were  run  using  5  pmols  NH4'^. 

of  the  nitrite  plus  nitrate  pools  was  analyzed  following  reduction  to  nitrous 
oxide  for  d^^N  analysis  using  the  denitrifier  method  (Sigman  et  al.  2001).  However, 
nitrite  was  commonly  present  at  insufficient  concentrations  (rarely  higher  than  1  pM)  for 
isotopic  analysis.  The  nitrous  oxide  produced  by  these  methods  was  analyzed  on  a  PDZ 
Europa  20-20  Isotope  Ratio  Mass  Spectrometer  at  the  UCDavis-SIF.  Standards  were 
made  by  equimolar  (5  mM)  dilutions  of  labeled  stocks  with  natural  abundance  nitrate 
(Thermo-Fisher,  ACS  grade)  resulting  in  standards  98.5,  75,  50,  25-atom  %  '^N.  These  5 
mM  stocks  were  diluted  1  mL  into  99  mL  of  DI  to  make  50  pM  working  solutions. 
Working  solutions  (400  pL)  were  injected  into  purged  denitrifier  cultures  to  generate  10 
nmol  of  N2O  gas.  All  samples  and  standards  were  run  at  20  nmols  N. 

Rate  Calculations 

Nitrification  rate  calculations  based  on  traditional  concentration  measurements  are 
sensitive  to  processes  that  consume  nitrate,  thus  measuring  a  net  rate.  The  advantage  of 
using  the  tracer  technique  to  calculate  nitrification  rates  is  that  the  technique  is  not  as 
sensitive  to  consumption  processes  because  despite  possible  isotopic  fractionation,  these 
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processes  do  not  significantly  alter  the  atom  percent  of  nitrogen  pool  at  tracer  levels. 
Processes  such  as  denitrification  or  ammonium  adsorption  to  sediments  have  an  isotope 
effect  in  the  per  mil  range,  which  is  much  smaller  than  the  percent  level  change  observed 
in  the  transfer  of  label-N  between  the  ammonium  and  nitrate  pools.  Therefore,  the  ratio  of 
in  nitrate  can  be  used  to  calculate  the  nitrification  rate.  The  ratio,  is 

calculated  from  atom  %  by: 
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Where  and  are  the  concentrations  of  heavy  and  light  N  isotopes.  These  ratios  can 
be  used  to  calculate  the  flux  of  N  from  the  labeled  ammonium  pool  to  the  nitrate  pool  by: 
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Where  and  are  the  initial  concentrations  of  labeled  and  unlabeled  nitrate,  Fnii  is 
the  flux  of  N  from  the  ammonium  pool  into  the  nitrate  pool.  At%NH4  is  the  atom  percent 
in  the  ammonium  pool  and  is  assumed  to  be  stable  (99.5-atom  %)  over  the  course  of 
the  experiment  (i.e.  ammonium  regeneration  from  other  unlabeled  sources  is  not 
significant),  k  denotes  the  first-order  rate  constant  for  the  uptake  of  ^"^N-containing 
nitrate,  and  a  =  1.005  is  the  ratio  of  rate  constants  ('"k/'^k)  for  consumption  of  unlabeled 
and  labeled  nitrate  (Sigman  et  al.  1999).  Equation  1  was  used  to  calculate  for 
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each  sample.  Equation  2  was  used  to  fit  the  data  using  Igor  Pro  (Wavemetrics) 

solving  for  Fnii  and  k.  Initial  values  for  Fnii  and  k  were  between  1  and  0.001  pmol  F'^  day' 
^  and  0.1  and  0.000001  day'\  respectively.  Fit  sensitivity  to  the  initial  guess  of  each 
parameter  was  tested  by  holding  one  variable  at  a  best  guess  and  varying  the  intial  value 
of  the  second  variable  over  many  orders  of  magnitude  (1,  0.1,  0.001  etc.)  until  the  fit 
degraded.  Rates  (Fnii)  in  from  the  model  fit,  in  moles  F''  hr',  were  converted  to  moles  g"' 
day"'  by  multiplying  by  the  volume  of  water  in  the  slurry  incubation  and  normalizing  to 
the  dry  weight  of  the  sediment. 

Molecular  Techniques 

DNA  was  extracted  from  1  g  (wet  weight)  of  sediment  using  lysozyme  and 
proteinase  K  digestions  followed  by  physical  disruption  by  successive  freeze-thaw  cycles 
and  phenolxhloroformdsoamyl  extraction  (Edwards  et  al.  2000).  Extracted  DNA  was 
precipitated  using  isopropyl  alcohol  and  3  M  sodium  acetate  with  final  resuspension  in  50 
pF  of  nuclease-free  water  or  Tris-EDTA  buffer  (TE).  The  extracts  were  quantified  at  260 
nm  using  a  Nanodrop  spectrophotometer  (Thermo-Fisher  Scientific)  and  diluted  using 
nuclease-free  water  to  appropriate  concentrations  for  amplification  and  quantification  by 
polymerase  chain  reaction  (PCR)  based  methods  (see  below). 

RNA  was  isolated  from  2  g  of  sediment  using  the  PowerSoil  Total  RNA  Isolation 
kit  (MoBio,  USA).  Extracts  were  treated  with  RQl  RNase  free  DNase  (Promega)  to 
remove  DNA.  cDNA  was  constructed  using  specific  antisense  primers  (1541R  for  16S 
rRNA,  and  the  antisense  primer  for  the  functional  genes.  Table  1)  and  the  Superscript  III 
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First-Strand  Synthesis  System  (Invitrogen)  following  the  manufacturer’s  instructions. 
Samples  without  the  reverse  transcription  were  used  as  negative  controls. 

Quantitative  RT-PCR  (qRT-PCR)  was  employed  to  determine  the  gene  transcript 
copy  numbers  of  bacterial  and  archaeal  16S  rRNA  and  ammonia  monooxygenase  subunit 
A  (amoA).  cDNA  was  constructed  as  detailed  above  and  quantified 
spectrophotometrically  using  a  Nanodrop  spectrophotometer  (Thermo-Fisher).  The 
cDNA  samples  were  diluted  to  100  ng  pL'^  and  amplified  using  an  iCycler  thermocycler 
(Bio-Rad,  Hercules,  CA)  and  a  cocktail  consisting  of  1  X  colorless  master  mix  (Promega, 
Madison,  WI),  10  uM  dNTPs  (Promega),  0.5  X  SBYR-Green  I  (Invitrogen,  Carlsbad 
CA),  20  uM  fluorescein  (Sigma),  1-3  mM  MgCl2  (Promega),  forward  and  reverse  primers 
(300-1000  uM;  Table  1)  and  1.25  units  of  CoTaq  Flexi  DNA  polymerase  (Promega).  All 
amplifications  used  optimized  magnesium  concentrations  and  thermal  programs  using 
(Table  1)  and  had  efficiencies  ranging  from  91  -  111%  for  all  reactions  resulting  in  linear 
standard  curves  (r  =  0.90  -  0.99).  Standard  curves  were  constructed  from  serial  dilutions 
(10^  to  10^  copies  pL'^  plus  a  water  blank)  of  linearized  p-Cem  T-Easy  plasmids 
(Promega)  containing  inserts  of  the  functional  gene  of  interest  (Table  1).  Plasmids  (0.5 
pg  pL"')  were  linearized  by  digestion  with  the  EcoRI  (New  England  Biolabs)  restriction 
enzyme  in  the  supplied  buffer  (2  hr  at  37  °C,  heat  inactivation  20  min.  at  65  °C). 
Restriction  digests  were  cleaned  using  the  SV  Wizard  PCR  purification  kit  (Promega), 
eluted  in  30  pE)  and  quantified  using  a  Nanodrop  spectrophotometer  (Thermo-Eisher). 
Melt  curves  were  performed  on  all  reactions  to  confirm  specificity. 
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Results 


Field  Chemistry 

The  groundwater  profiles  from  the  field  effort  confirmed  the  geochemical  context 
of  the  system,  detailed  in  previous  reports  (Chapter  2  and  3).  In  short,  the  shallow 
groundwater  system  is  composed  of  three  different  parcels  of  water  and  the  interfaces 
between  them.  1)  The  surface  parcel  extends  to  about  1  m  in  depth  and  is  brackish 
(salinity  less  than  5),  being  composed  of  a  mixture  of  groundwater  and  seawater  from 
overtopping  of  the  bay  water.  Oxygen  and  nitrate  are  present  in  concentrations  ranging 
from  30  to  90  pM  O2  and  0  to  30  pM  NO3'  with  the  highest  concentration  occurring  in 
the  deeper  portions  of  the  parcel.  2)  In  the  fresh  groundwater  plume,  nitrate  and  oxygen 
rise  to  maxima  (97.5  and  156  pM  respectively)  before  decreasing  approaching  the  deep 
salt  wedge.  3)  The  deep  salt  wedge  is  composed  of  recirculated  seawater  and  is  high  in 
NH4'^  but  low  in  O2  and  NO3'  (Figure  1). 

^^N-Ammonium  Incubations 

Sediment  from  three  depths  (0.6  m,  2  m,  and  4  m)  was  incubated  in  a  12  °C 
shaker  amended  with  ^^N-labeled  ammonium  with  and  without  allylthiourea  (ATU),  an 
inhibitor  of  bacterial  ammonia-oxidation  (Hooper  and  Terry  1973).  The  shallow 
sediments  from  0.6  m  depth  showed  a  linear  decrease  of  ammonium  concentration  (ca.  5 
pM)  over  4  days  while  nitrate  concentration  increased  by  3  pM  (Figure  2).  The  loss  of 
labeled  ammonium  was  accompanied  by  an  increase  in  the  ^^N-labeling  of  the  nitrate 
pool  up  to  33-atom  %.  This  increase  in  isotopic  value  is  consistent  with  the  observed 
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increase  in  the  nitrate  pool  resulting  from  nitrification  of  the  labeled  ammonium  pool. 

The  nitrite  pool  was  negligible  throughout  the  experiment.  Similar  patterns  were 
observed  in  the  ATU  treatment,  with  approximately  10  pM  of  ammonium  consumed  and 
3  pM  nitrate  produced  while  the  ^^N-labeling  of  the  nitrate  pool  increased  to  ca.  30-atom 
%.  Nitrification  rates  were  calculated  as  0.019  ±  0.006  pmol  L'^  hr'  or  456  pmol  g"' 
sediment  day"'  given  the  100  mL  aqueous  and  100  g  sediment  volumes  for  incubation 
(Figure  5).  This  calculation  was  not  sensitive  to  changes  in  consumption  rates  (k  in 
equation  2)  over  seven  orders  of  magnitude  start  values  (10  to  1  x  10'^).  A  linear  fit  was 
also  tested  yielding  a  rate  370  pmol  g"'  sediment  day’'.  Nitrification  rates  in  the  presence 
of  ATU  reached  0.018  ±  0.001  pmol  U'  hr'  or  432  pmol  g"'  sediment  day"'  (218  pmol  g"' 
sediment  day"',  linear  fit).  Published  N-flux  rates  from  the  environment  are  often  reported 
in  mol  m'  time'  .  Therefore  to  compare  these  area-normalized  estimates  to  our  volume 
derived  estimates  some  assumptions  must  be  made.  First  to  convert  from  an  area  to  a 
volume  we  calculate  a  diffusion  length  using  Fick’s  law: 

L  =  (3) 

where  D  is  the  diffusion  coefficient  and  t  is  time.  This  calculation  yields  a  diffusion 
length  of  1.59  cm  day"'  for  nitrate  in  a  12  °C  seawater  solution.  Converting  from  a 
volume  to  weight  of  sediment  requires  an  assumption  for  the  sediment  density  (2.2  g  mL' 
')  and  the  porosity  (0.9). 

The  2  m  incubations  showed  a  loss  of  ammonium  (ca.  9  pM)  and  production  of 
nitrate  (ca.  7  pM).  This  increase  in  nitrate  was  associated  with  a  corresponding  increase 
in  '^N  up  to  29-atom  %,  an  increase  that  only  accounts  for  ca.  4  pM  of  the  observed 
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increase  assuming  all  N  originates  in  the  labeled  ammonium  pool  (Figure  3).  Unlike  the 
shallow  incubations,  the  2  m  incubation  with  ATU  showed  only  an  initial  loss  of 
ammonium  in  the  first  24  hr  (ca.  3  pM)  and  remained  stable  throughout  the  remainder  of 
the  incubation.  Nitrate  increased  by  ca.  4  pM  in  the  first  24  hr  and  remained  stable 
thereafter  (ca.  10  pM).  The  isotopic  signature  of  the  nitrate  pool  increased  only  slightly 
during  the  incubation  from  <l-atom  %  to  2.8-atom  %.  This  increase  in  accounts  for 
0.2  pM  transfer  of  N  from  the  ammonium  to  nitrate  pools.  Nitrification  rates  at  this  depth 
were  calculated  as  0.016  ±  0.009  nmol  L'^  hr'  or  384  pmol  g'^  day’'  (327  pmol  g"' 
sediment  day"',  linear  fit).  The  nitrification  rate  in  the  ATU  addition  0.002±  0.0007  nmol 
U'  hr'  or  22  pmol  g'^  day"'  (22  pmol  g"'  sediment  day"',  linear  fit). 

The  4  m  tracer  incubation  revealed  a  different  chemical  pattern  than  the  two 
shallower  incubations.  Ammonium,  at  this  depth,  decreased  by  ca.  2  pM  over  the  first  24 
hr  then  remained  steady  for  the  duration.  Nitrate  increased  by  ca.  2  pM  over  the  first  24 
hr  remaining  at  ca.  4  pM  for  the  rest  of  the  time  course.  This  increase  in  nitrate  was  not 
associated  with  an  increase  in  '^N-atom  %  of  the  nitrate  pool.  In  fact,  the  '^N-atom  %  of 
the  nitrate  pool  at  this  depth  was  less  1-atom  %  for  the  duration  of  the  experiment. 

Gene  Expression 

Quantitative  RT-PCR  was  used  to  determine  the  expression  level  of  specific 
functional  genes  (bacterial  amoA  or  archaeal  amoA)  as  well  as  bacterial  and  archaeal  16S 
rRNA  genes.  mRNA  was  extracted  from  the  sediment  and  (reverse)  transcribed  into 
cDNA  for  each  incubation  in  the  '^NH4'^  tracer  experiment.  It  must  be  noted  that  our 
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results  will  be  discussed  as  relative  changes  in  transcript  abundance.  Normalization  of 
mRNA  expression  is  challenging  (Bustin  2000)  and  some  authors  ignore  normalization 
entirely  in  environmental  studies  (eg.  Lam  et  al.  2009).  Typically,  mRNA  expression  is 
normalized  to  expression  levels  of  a  gene  or  suite  of  genes  that  is  continually  expressed  at 
a  constant  level  without  sensitivity  to  culture  conditions,  such  genes  are  commonly  called 
housekeeping  genes.  Unfortunately,  choosing  an  appropriate  housekeeping  gene  in  an 
environmental  setting,  as  opposed  to  a  pure  culture  or  tissue  culture,  is  difficult.  We 
initially  amplified  the  rRNA  genes  with  the  idea  that  these  genes,  which  can  constitute 
85-90%  of  the  total  cellular  mRNA  (Bustin  2000),  may  provide  an  appropriate  basis  for 
normalization.  Unfortunately,  16S  rRNA  transcripts  were  not  constant  over  the  course  of 
the  incubation  (Figure  2,  3,  4),  though  total  RNA  yield  from  extractions  were  relatively 
constant.  The  changing  abundance  of  16S  rRNA  transcripts  may  be  due  to  the  disruption 
of  sampling  (a  bottle  effect)  or  the  simulation  of  cells  through  the  addition  of  substrate 
(NH4'^).  Therefore,  the  patterns  discussed  below  are  internally  comparable  as  relative 
abundance  but  extrapolation  to  other  environments  or  other  experiments  is  not  advisable. 

The  0.6  m  incubations  showed  a  nearly  3-fold  increase  in  archaeal  and  bacterial 
16S  cDNA  abundance  (Figure  2).  The  archaeal  16S  cDNA  increased  over  the  first  24  hr 
of  the  experiment  from  1.4  x  10^  to  4.9  x  10^  copies  g'\  and  remained  elevated 
throughout  the  remainder  of  the  time  course.  Bacterial  16S  cDNA  showed  a  similar 
increase  from  1.6  x  10^  to  4.6  x  10^  copies  g'^  by  72  hr.  Archaeal  amoA  cDNA  increased 
in  the  first  24  hr,  from  9. lx  10  to  5.4  x  10  copies  g'  ,  then  steadily  decreased  over  the 
remainder  of  the  incubation.  Bacterial  amoA  cDNA  displayed  a  similar  pattern. 
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increasing  from  2.3  x  10^  to  9.4  x  copies  g"'  within  24  hr,  5-fold  lower  than  the 
archaeal  gene.  Addition  of  the  ATU  inhibitor  resulted  in  little  change  in  the  expression 
pattern  of  these  genes.  Archaeal  16S  rRNA  and  amoA  cDNA  abundance  was  not  affected 
by  the  presence  of  ATU.  Notably,  the  bacterial  16S  rRNA  cDNA  maintained  levels  about 
half  of  the  uninhibited  experiment.  Likewise,  bacterial  amoA  expression  in  the 
ammonium  plus  ATU  incubations  remained  less  than  5x10  copies  g'  throughout  the 
time  course. 

The  incubations  from  2  m  revealed  a  different  expression  pattern  than  the  shallow 
incubations  (Figure  3).  Bacterial  16S  rRNA  cDNA  abundance  ranged  from  1.4  x  10^  to 
3.9  X  10^  copies  g'\  Archaeal  16S  rRNA  cDNA  abundance  showed  a  similar  range,  1.3  x 
10^  to  3.0  X  10^  copies  g"'.  Bacterial  amoA  cDNA  ranged  from  4.3  x  10^  to  6  x  10^  copies 
g'\  peaking  at  48  hr  before  decreasing  to  1.2  x  10^  copies  g’'  by  96  hr.  Archaeal  amoA 
abundance  was  not  as  dynamic  as  its  bacterial  counterpart,  ranging  1.4  x  10  to  3.2  x  10 
copies  g In  the  presence  of  ATU,  bacterial  16S  rRNA  cDNA  abundance  ranges  from 
9.1  X  10^  to  5.0  X  10^  copies  g'^  but  tends  to  decrease  over  the  course  of  the  experiment. 
Archaeal  16S  rRNA  cDNA  abundance  increases  from  1.4  x  10^  to  4.4  x  10^  copies  g"'  in 
the  first  24  hr  and  remains  at  this  abundance  for  the  remaining  incubation.  Interestingly, 
in  the  presence  of  ATU  the  abundance  of  bacterial  amoA  cDNA  is  less  than  5.0  x  10 
copies  g'^  Archaeal  amoA  cDNA  abundance  is  also  low. 

Gene  expression  in  the  4  m  incubations  is  unlike  the  two  shallower  incubations. 
Bacterial  16S  rRNA  cDNA  abundance  is  lower  than  found  in  the  other  incubations, 
ranging  from  5.1  x  10^  to  8.3  x  10^  copies  g'^  while  archaeal  16S  rRNA  cDNA  abundance 
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was  comparable  to  other  incubations,  6.7  x  10^  to  9,1  x  10^  copies  g'^  (Figure  4).  Unlike 
the  shallow  incubation  both  amoA  cDNA  abundance  remain  low  during  the  incubation. 
Bacterial  amoA  abundance  varies  from  5.4  x  10  to  8.4  x  10  copies  g  while  archaeal 
amoA  ranges  from  5.6  x  10  to  8.9  x  10  copies  g'  under  the  same  conditions.  Gene 
express  of  all  four  genes  was  consistent  between  the  ammonium  and  ammonium  plus 
ATU  incubations. 


Discussion 

The  three  groundwater  parcels,  based  on  salinity,  of  the  Waquoit  Bay  aquifer  have 
previously  been  described  in  detail  (Charette  et  al.  2005;  Kroeger  and  Charette  2008). 
Briefly  they  include:  (1)  a  surface  salt  zone  that  is  the  result  of  mixing  between  meteoric 
waters  and  overtopping  bay  water  from  tidal  and  wind  driven  circulation,  (2)  a  fresh 
groundwater  and  (3)  a  tidally  driven  deep  salt  wedge.  These  parcels  represent  three 
distinct  geochemical  environments  within  which  nitrogen  transformations  can  occur. 
Within  these  three  zones  previous  studies  have  shown  the  presence  of  different  microbial 
communities  (Chapter  2)  and  potentially  different  biogeochemical  processes  in  regard  to 
nitrogen  cycling  (Chapter  3;  (Kroeger  and  Charette  2008;  Spiteri  et  al.  2008).  This  study 
aims  to  estimate  the  rate  of  nitrification  within  each  parcel  by  following  the  transfer  of 
from  the  NH4'^  to  NO3'  pools.  Furthermore,  gene  expression  of  relevant  functional 
genes  are  quantified  over  the  time  course  of  the  experiment  with  the  goal  of  elucidating 
which  groups  of  organisms  may  be  active. 
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Surface  Meter 

In  earlier  work,  a  large  community  of  AOAs  was  observed  in  the  surface  meter 
(Chapter  2).  Here  we  show  through  quantification  of  archaeal  amoA  transcripts  that  these 
AOA  are  active.  Furthermore,  the  transfer  of  the  tracer  from  the  ammonium  to  the 
nitrate  pools  confirms  nitrification  at  this  depth.  ATU  is  believed  to  specifically  inhibit 
ammonia  oxidation  by  chelating  the  copper  atom,  obstructing  the  active  site  of  the 
ammonia  monooxygenase  enzyme  (Ginestet  et  al.  1998),  however,  in  complex 
community  and  environment  such  as  the  sediments  ATU  may  have  unforeseen  reactions. 
For  example,  ATU  may  be  titrated  by  the  presence  of  aqueous  or  sedimentary  copper  or 
perhaps  by  interactions  with  organisms  and  molecules  other  than  those  intended.  These 
interactions  are  difficult  to  measure  or  predict  in  complex  environmental  settings. 

The  resulting  nitrification  rates  (456  pmol  g"'  sediment  day'^)  are  comparable  to 
those  reported  for  estuarine  sediments  (400  to  9000  nmol  m'  day'  or  1 1  pmol  g' 
sediment  day'^  to  257  pmol  g"'  sediment  day'^)  (Henriksen  et  al.  1981;  Caffrey  et  al. 
2003).  Compared  to  carbon  rich  environment  like  salt  marsh  sediment  our  nitrification 
rate  estimates  are  100  to  1000  times  slower  (Dollhopf  et  al.  2005).  The  activity  observed 
in  this  zone  is  important  because  much  of  the  groundwater  discharge  may  transit  through 
a  brackish  zone  similar  to  this  prior  to  discharge.  However,  given  the  short  residence  time 
of  water  in  the  upper  meter  (Michael  et  al.  2005),  and  the  observed  1-5  pM  NH4'^ 
concentrations  (Chapter  3)  it  is  possible  that  the  biological  rates  are  too  slow  to 
completely  consume  ammonium  either  transported  in  from  the  bay,  remineralized  from 
organic  matter  or  release  through  cation  exchange  as  saline  water  invades.  Oxidation  of 
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ammonium  in  this  zone  would  result  in  a  lighter  (smaller  8^^N)  nitrate  pool,  which  is  not 
the  pattern  that  has  been  observed  (Chapter  3).  However,  it  is  likely  that  the  observed 
ratio  of  8^^0/8^^N  (Chapter  3)  is  not  solely  a  signal  of  denitrification  in  this  or  other  - 
groundwater  environments  (Bottcher  et  al.  1990).  Likely,  these  environments  harbor  both 
nitrifying  and  denitrifying  communities  whose  coupled  metabolism  produce  the  observed 
ratio  of  N  and  O  isotopes  in  the  nitrate  pool. 

Fresh  Groundwater  Plume 

The  fresh  groundwater  plume  has  been  described  as  a  location  where  AOBs  may 
be  more  important  than  AOA  due  to  the  relative  abundance  of  their  respective  amoA 
genes  (Chapter  2).  Our  data  here  corroborates  this  hypothesis  by  showing  both  an 
abundance  of  active  AOBs  (through  bacterial  amoA  transcripts)  and  transfer  of  the 
isotopic  label  from  the  reduced  to  oxidized  nitrogen  pools.  The  potential  rates  of 
nitrification  at  2  m  are  equal  to  that  of  shallow  aquifer  sediments.  Interestingly,  the 
presence  of  ATU  completely  blocks  nitrification  at  this  depth,  suggesting  that  nitrifier 
community  at  this  depth  is  sensitive  to  the  inhibitor.  We  interpret  this  sensitivity  as  an 
indication  that  the  active  nitrifier  community  is  mainly  bacterial  as  opposed  to  archaeal  in 
the  freshwater  plume.  Interestingly,  a  drop  in  the  expression  level  of  bacterial  amoA 
genes  also  occurred  in  the  presence  of  ATU.  This  is  unexpected  because  ATU  is  a  copper 
chelator,  binding  to  and  physically  blocking  the  active  site  of  the  ammonia 
monooxygenase  enzyme.  If  the  enzyme  is  blocked  it  does  not  function  and  it  would  be 
expected  that  the  expression  levels  of  the  gene  encoding  for  the  needed  metabolic 
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enzyme  would  remain  elevated.  However,  ATU  could  have  other  effects  in  the  complex 
environment  represented  by  the  bottle  incubations.  One  possible  explanation  is  that  ATU 
altered  a  factor  that  promotes  the  expression  of  the  functional  gene.  Another  explanation 
could  be  that  the  ATU  altered  some  other  required  enzyme  system  in  the  organisms  at 
this  depth  that  resulted  in  the  population  no  longer  metabolizing  ammonia.  A  third 
possibility  could  be  that  ATU  reacted  with  the  sediment  matrix  at  this  depth,  perhaps 
cleaving  the  inhibitor  or  otherwise  altering  the  chemical  environment  that  resulted  in  the 
organisms  no  longer  processing  the  added  ammonia.  At  this  point,  it  is  unclear  why  the 
addition  of  ATU  resulted  in  the  loss  of  the  bacterial  amoA  gene  expression  signal.  The 
archaeal  amoA  expression  at  this  depth  was  low  both  in  the  presence  and  absence  of  the 
inhibitor  indicating  that  archaea  are  probably  not  the  major  ammonia-oxidizers  within  the 
core  of  the  freshwater  nitrate  plume. 

The  presence  of  an  active  population  of  AOB  at  this  depth  is  interesting  because 
ammonium  is  not  observed  here  as  it  is  in  the  surface  and  deep  salt  wedge.  However, 
there  are  possible  sources  of  ammonia  in  this  zone  such  as  in  situ  remineralization  of 
organic  matter  or  perhaps  the  advection  of  ammonium-bearing  water  into  the  fresh 
groundwater  plume.  The  sands  as  this  depth  in  the  aquifer  are  carbon  poor,  thus  ammonia 
generation  ought  to  be  small.  Likewise,  an  advecting  ammonia  plume  has  not  been 
observed  at  the  current  study  site,  though  such  a  plume  has  been  mapped  about  0.5  km 
east.  If  an  ammonium  plume  similar  to  the  one  in  the  east  were  to  enter  this  zone  it  is 
unlikely,  given  the  reaction  rates  and  high  concentrations  of  ammonia  (>50  pM),  that  the 
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ammonia  would  be  efficiently  attenuated  in  this  zone  though  this  is  dependent  on  the  rate 
of  advection. 

Deep  Salt  Wedge 

The  deep  salt  wedge  has  been  shown  to  harbor  a  large  abundance  of  archaeal 
amoA  genes  (Chapter  2).  This  has  been  interpreted  as  a  population  of  AOA.  Curiously, 
our  examination  has  not  demonstrated  the  presence  of  an  active  ammonia-oxidizing 
population.  Gene  expression  for  both  amoA  versions  was  low  and  there  was  not  a  transfer 
of  label  between  the  ammonium  and  nitrate  pools.  The  ammonium  from  the  deep  salt 
wedge  has  been  described  as  mixing  nearly  conservatively  entering  the  freshwater  plume 
from  below  with  some  evidence  of  ammonium  loss  (Kroeger  and  Charette  2008).  The 
data  presented  above  is  not  consistent  with  an  AOA  population  ready  to  consume 
ammonium  as  soon  as  it  is  available.  If  this  population  was  tuned  to  ammonium 
metabolism,  the  addition  of  the  tracer  ought  to  have  stimulated  consumption.  Of  course, 
our  experimental  procedure  may  not  mimic  the  environment  accurately  enough  to 
maintain  a  viable  population  of  AOA  from  this  depth.  There  are  a  few  possible 
explanations  for  this  inactivity  in  the  deep  AOA  population,  first  perhaps  these  organisms 
are  limited  by  some  other  required  growth  element  that  is  not  limiting  in  the  surface 
sediments,  perhaps  this  population  is  more  sensitive  to  the  handling  required  to  set  up 
cultures,  especially  changing  oxygen  concentrations  or  perhaps  the  lag  time  of  this 
population  is  greater  than  the  96  hr  of  the  experiment.  Alternatively,  anammox  processes, 
either  biological  or  by  oxidation  with  Mn-oxides  could  also  account  for  the  loss  of 
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ammonium  at  this  depth  (Kroeger  and  Charette  2008).  Biological  anammox  processes 
may  or  may  not  be  linked  to  archaeal  nitrification  as  has  been  suggest  for  the  Black  Sea 
(Lam  et  al.  2007)  and  Peruvian  upwelling  (Lam  et  al.  2009).  Future  work  and  technique 
development  will  be  able  to  elucidate  this  relationship. 

Conclusions 

This  work  confirms  the  presence  of  different  communities  of  ammonia-oxidizers 
throughout  the  upper  meters  of  the  Waquoit  Bay  groundwater  system.  The  upper  meter  of 
the  aquifer  harbors  an  active  population  of  archaeal  ammonia-oxidizers,  which  appear  to 
be  responsible  for  the  production  of  oxidized  nitrogen  species  at  potential  rates 
comparable  to  estuarine  sediments.  In  contrast,  bacteria  appeared  to  be  the  active 
nitrifiers  in  the  freshwater  plume.  Further  work  is  needed  to  determine  the  geochemical 
role  of  the  deep  archaeal  population,  which  did  not  appear  active  in  nitrification  at  the 
time  of  our  sampling.  Prior  evidence  suggests  that  a  fraction  of  these  deep  archaea 
possess  the  amoA  gene;  however,  the  gene  was  not  expressed  under  the  incubation 
conditions  in  this  study.  In  fact,  no  nitrification  was  evident  within  the  deep  salt  wedge 
even  when  NH4^  was  added.  Closer  examination  of  this  deep  population  of  archaea, 
including  cultivation  could  shed  new  light  on  the  metabolism  of  these  organisms. 
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A 


Figure  1.  (A)  Head  of  the  bay,  Waquoit  Bay,  Cape  Cod,  MA.  Contours  represent  10  ft. 
intervals  and  hatehed  area  denotes  wetlands.  (B)  The  eurrent  and  previous  study  locations 
are  depicted.  Cape  Cod  Commission  wells  (CCCl,  CCC2),  transect  from  Charette  et  al. 
2005,  2006  (A  to  A’)  and  Kroeger  et  al.  (2008)  (A  to  A’,  B  to  B’  and  C  to  C’).  Mulligan 
and  Charette  (2006)  extended  the  B  to  B’  transect  east  and  west. 
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A)  B) 


C) 


D) 


Figure  2.  Shallow  (0.6  m)  incubation  time  course.  A)  Ammonium  ()j,M,  circles),  nitrate 
()j,M,  squares)  and  '^N03  atom  %  (triangles)  are  plotted  against  time  during  an  50  pM, 
99.5  atom  %  slurry  incubation  .  Ammonium  decreases  while  nitrate  and 

increase  over  the  coarse  of  the  experiment.  B)  Ammonium  (pM,  circles),  nitrate  (pM, 
squares)  and  atom  %  (triangles)  development  in  the  ammonium  plus  ATU 
incubations.  C)  Gene  expression  (qRT-PCR)  data  for  archaeal  and  bacterial  amoA  and 
16S  rRNA  cDNA  during  the  time  course  of  the  labeled  ammonium  incubation.  D)  Gene 
expression  (qRT-PCR)  data  for  archaeal  and  bacterial  amoA  and  16S  rRNA  cDNA  during 
the  time  course  of  the  labeled  ammonium  plus  ATU  incubation. 


138 
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Figure  3.  Tracer  incubation  using  sediments  and  water  collected  from  2  m.  A) 
Ammonium  (gM,  circles),  nitrate  (gM,  squares)  and  atom  %  (triangles)  are  plotted 
against  time  during  an  50  pM,  99.5  atom  %  slurry  incubation  .  Ammonium 

decreases  while  nitrate  and  '^NOa  increase  over  the  coarse  of  the  experiment.  B) 
Ammonium  (pM,  circles),  nitrate  (pM,  squares)  and  ^^NOa  atom  %  (triangles) 
development  in  the  ammonium  plus  ATU  incubations.  C)  Gene  expression  (qRT-PCR) 
data  for  archaeal  and  bacterial  amoA  and  16S  rRNA  cDNA  during  the  time  course  of  the 
labeled  ammonium  incubation.  D)  Gene  expression  (qRT-PCR)  data  for  archaeal  and 
bacterial  amoA  and  16S  rRNA  cDNA  during  the  time  course  of  the  labeled  ammonium 
plus  ATU  incubation. 
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Figure  4.  Tracer  incubation  using  sediments  and  water  collected  from  4  m.  A) 
Ammonium  (gM,  circles),  nitrate  (gM,  squares)  and  atom  %  (triangles)  are  plotted 
against  time  during  an  50  pM,  99.5  atom  %  slurry  incubation.  Ammonium 

decreases  while  nitrate  and  '^NOa  increase  over  the  coarse  of  the  experiment.  B) 
Ammonium  (pM,  circles),  nitrate  (pM,  squares)  and  ^^NOa  atom  %  (triangles) 
development  in  the  ammonium  plus  ATU  incubations.  C)  Gene  expression  (qRT-PCR) 
data  for  archaeal  and  bacterial  amoA  and  16S  rRNA  cDNA  during  the  time  course  of  the 
labeled  ammonium  incubation.  D)  Gene  expression  (qRT-PCR)  data  for  archaeal  and 
bacterial  amoA  and  16S  rRNA  cDNA  during  the  time  course  of  the  labeled  ammonium 
plus  ATU  incubation. 
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A)  B) 


Time  (hr?)  Tim«  (hr?) 


Figure  5.  Nitrification  rates  as  calculated  by  15N-label  transfer  from  ammonium  to 
nitrate  pools.  A  simple  box  model  was  fit  to  the  ratio  of  15N/14N  over  time  using  Igor 
Pro.  A)  0.6  m  incubations,  B)  2  m  incubations,  C)  0.6  m  plus  ATU  incubations  and  D)  2 
m  plus  ATU  incubations. 
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Introduction 


Throughout  coastal  Massachusetts  and  on  Cape  Cod  in  particular,  groundwater 
quality  is  heavily  impacted  by  population  density  (as  measured  by  individuals  or  housing 
per  acre)  with  the  principle  impact  on  nitrogen  loading  (Persky  1986;  Valiela  et  al.  1992; 
Sham  et  al.  1995;  McClelland  et  al.  1997;  McClelland  and  Valiela  1998;  Cole  et  al. 

2006).  This  nitrogen,  initially  in  reduced  forms,  is  transformed  primarily  to  nitrate  and 
enters  coastal  waters  as  new  nitrogen  through  surface  waters  and  groundwater  discharge. 
Interest  in  groundwater  as  a  freshwater  input  to  coastal  waters  has  increased  in  recent 
years  due  to  of  the  projected  magnitude  of  the  water  flux,  chemical  constituents 
associated  with  this  flux,  and  potential  impact  of  those  constituents  compared  to  overland 
runoff  (rivers  and  storm  water). 

The  subterranean  estuary  at  Waquoit  Bay  is  an  appropriate  and  tractable  model 
system  in  which  to  study  microbe-mediated  transformations  of  anthropogenic  nitrogen 
for  several  reasons  related  to  both  the  physical  and  chemical  environment.  First,  the 
presence  of  a  sandy  aquifer  is  common,  especially  in  the  areas  where  glaciation  has 
altered  the  underlying  strata.  Second,  the  large  population  burden  on  natural  resources 
and  anthropogenic  influence  on  the  aquifer  system  is  common.  This  setting  should  serve 
as  an  analog  to  many  temperate  coastal  watersheds  including  places  such  as  Long  Island, 
NY.  In  addition,  there  is  a  good  deal  of  relevant  previous  and  ongoing  research  conducted 
at  this  site  to  provide  context,  although  very  little  microbiology  has  been  done  in  this 
system. 
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Our  work  has  shown  the  presence  of  the  archaeal  ammonia  monooxygenase  gene 
(amoA)  at  all  sampling  depths  while  the  bacterial  version  of  the  gene  is  only  abundant  in 
the  core  of  the  freshwater  plume  (Chapter  2).  Interestingly,  our  data  show  that  archaeal 
amoA  expression  and  the  resulting  nitrification  dominate  in  the  surface  meter  while  the 
bacterial  counterparts  are  more  abundant  in  the  core  of  the  freshwater  plume  (Chapter  4). 

Denitrification  within  the  groundwater  system  has  also  been  demonstrated  using 
chemical  and  gene  abundance  data  (Chapter  3).  The  isotopic  ratios  of  the  nitrate  pool, 
S'^Nno3  and  8^^0no3,  increase  while  the  nitrate  concentration  decreases  moving  from  1  m 
toward  the  surface.  Over  the  same  interval  denitrification  genes  have  been  detected,  with 
the  gene  encoding  for  the  copper-based  nitrite  reductase  enzyme  dominant  at  this  depth 
(Chapter  3).  Interestingly,  tracer  experiments  coupled  with  gene  expression  data 
demonstrate  active  nitrification  may  also  occur  over  the  same  depth  range  as 
denitrification  (Chapter  4). 

Ladderane  lipids,  a  biomarker  for  anammox  organisms,  and  crenarchaeol,  a 
biomarker  associated  with  crenarchaeota,  have  also  been  detected  within  the  aquifer  with 
peaks  in  both  lipids  associated  with  the  mixing  zones  between  the  freshwater  plume  and 
the  two  saline  waters  (Saenz,  personal  communication).  Thus,  this  subterranean  estuary 
has  biological  markers  indicating  that  a  variety  of  nitrogen  transformation  pathways 
within  the  system,  many  of  which  may  be  overlapping  in  depth  or  co-occurring.  The 
discussion  that  follows  is  an  effort  to  place  our  geomicrobiological  data  and 
interpretations  in  the  context  of  previous  geochemical  and  modeling  work  at  Waquoit 
Bay. 
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The  Bay  Proper 


Waquoit  Bay  is  a  semi-enclosed  estuary  located  on  the  Southwest  coast  of  Cape 
Cod,  in  the  town  of  Falmouth,  MA.  The  bay  is  contained  by  the  landmasses  of  Cape  Cod 
to  the  north  and  Washburn  Island  to  the  south.  The  bay  has  direct  communication  with 
Vineyard  Sound  to  the  south  through  a  small  (ca  200  m)  straight  bay  between  Washburn 
Island  and  South  Cape  Beach  State  Park  (Cape  Cod,  MA)  to  the  east.  Indirect 
communication  with  Vineyard  Sound  also  occurs  through  a  small  but  navigable 
waterway  in  the  northwestern  part  of  the  bay,  traveling  along  the  western  edge  of 
Washburn  Island  through  the  narrow  Eel  Pond. 

The  geologic  setting  of  the  western  end  of  Cape  Cod,  and  thus  Waquoit  Bay, 
consists  of  sand  and  silt  deposited  on  a  gravel  bed  during  the  retreat  of  the  Wisconsin 
glaciation  (Oldale  1976;  Oldale  and  O'Hara  1984).  These  highly  permeable  soils  allow 
for  high  groundwater  recharge  from  precipitation  (ca.  1 14  cm  y'^)  (Walter  et  al.  2004) 
leading  to  groundwater  fed  streams  and  rivers  with  little  overland  runoff  (Cambareri  and 
Eichner  1998).  Water  table  height  in  the  surrounding  area  shows  some  seasonality  with 
peak  water  table  heights  in  April,  following  a  net  recharge  from  late  autumn  to  early 
spring,  and  water  table  lows  occurring  in  December,  following  net  evapotranspiration 
from  late  spring  to  early  autumn  (Michael  et  al.  2005).  Direct  groundwater  discharge  is 
estimated  to  account  for  34%  of  the  freshwater  flux  into  Waquoit  Bay,  while  atmospheric 
deposition  is  responsible  for  11%  of  the  freshwater  (Cambareri  and  Eichner  1998).  The 
majority  of  freshwater  supply  to  Waquoit  Bay  comes  by  way  of  two  rivers,  the  Quashnet 
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and  Child’s  rivers,  which  account  for  up  to  55%  of  the  freshwater  flux  (Cambareri  and 
Eichner  1998).  Precipitation-  and  radium-derived  recharge  estimates  for  total 
groundwater  flux  in  to  Waquoit  bay  range  from  23,700  to  37,000  m  d'  (Cambareri  and 
Eichner  1998;  Charette  et  al.  2001).  This  flux  includes  both  a  freshwater  flux  ranging 
from  0.14  -  3.93  m  day'  m'  shore  line  and  a  recirculated  seawater  flux  of  0.56  -  6.88  m 
day'^  m'^  shore  line,  calculated  by  Darcy’s  law  or  specific  seepage  measurements 
(Cambareri  and  Eichner  1998;  Michael  et  al.  2003;  Michael  et  al.  2005;  Mulligan  and 
Charette  2006).  The  residence  time  of  water  in  Waquoit  Bay  is  9.4  days,  as  determined 
by  radium  measurements  (Charette  et  al.  2001).  The  tidal  range  within  the  bay  is  1.1  m. 
Waquoit  Bay  is  a  shallow  with  an  average  depth  of  1  m  over  it  3  km  area  (Michael  et  al. 
2003)  draining  a  watershed  of  approximately  65  km^  (Charette  et  al.  2001). 

The  watershed  is  densely  populated  with  190  people  km'  (Kroeger  et  al.  2006; 
Spiteri  et  al.  2008),  and  therefore  a  significant  nutrient  input  derived  from  private  septic 
systems  within  the  residential  area  (Valiela  et  al.  1992).  Nutrient  input  in  the  bay  has 
caused  a  decrease  in  the  coverage  area  of  sea  grass  as  well  as  harming  the  population 
density  of  shellfish  within  the  bay  (Valiela  et  al.  1990;  Bowen  and  Valiela  2001; 
Evgenidou  and  Valiela  2002).  The  majority  of  this  nutrient  input  into  the  bay  is  due  to 
wastewater  (10,998  kg  yr'\  48%)  followed  by  atmospheric  deposition  (6923  kg  yr', 
30%)  and  fertilizer  (3442  kg  yr',  15%)  (Valiela  et  al.  1997). 
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Groundwater  Hydrology 


On  a  broad  scale,  groundwater  input  into  Waquoit  Bay  supplies  a  significant 
amount  of  water  and  dissolved  constituents  to  the  shallow  bay.  Interestingly,  from  the 
current  study  and  recent  studies  at  different  locations  within  the  bay,  the  composition  and 
flux  of  the  groundwater  constituents  varies  across  the  bay  (Figure  1).  Much  of  the 
previous  work  has  focused  on  an  area  in  the  northern  section  of  the  bay  where  the  Cape 
Cod  Commission  has  installed  wells  for  groundwater  sampling  (CCC  line)  and 
subsequent  transect  surveys  have  extended  this  line  through  the  intertidal  zone.  These 
wells  are  placed  in  a  topographic  low,  bound  to  the  east  and  west  by  6  m  tall  ridges  and  to 
the  north  by  Bog  pond  and  Bourne  pond.  The  coring  and  transect  sites  of  the  current 
study  are  approximately  50  m  west  of  the  CCC  line  but  within  the  bounds  of  previous 
work  and  mapping  projects  in  the  head  of  the  bay.  In  general,  the  groundwater  system  at 
Waquoit  Bay  is  composed  of  three  main  parcels  of  water  (Figure  2),  1)  the  surface  salt 
zone  (also  referred  to  as  a  brackish  zone),  2)  the  freshwater  plume  and  3)  the  deep  salt 
wedge.  The  major  difference  between  the  current  study  site  (west  of  the  CCC  line)  and 
previous  work  done  at  the  CCC  line  and  east,  lies  in  the  in  the  composition  of  the 
freshwater  plume  and  the  absence  of  an  additional,  reduced  freshwater  plume  at  our  site. 
The  freshwater  plume  at  the  current  site,  west  of  the  CCC  line,  has  similar  oxygen 
concentrations  compared  to  previous  studies  (average  86  pM  O2  at  core  B  vs.  78  pM, 
Kroeger  and  Charette  2008),  however,  ammonium  is  largely  absent.  West  of  the  CCC 
line,  ammonium  is  only  observed  in  the  surface  salt  zone  (<5  pM)  and  in  the  deep  salt 
wedge  (up  to  53  pM).  Ammonium  is  not  observed  within  the  freshwater  parcel  along  the 
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transect  but  could  exist  more  landward.  The  source  of  this  ammonium  in  the  two  salt 
zones  is  likely  the  aquifer  sands  due  to  ion  exchange  caused  by  the  intruding  saltwater 
(Charette  and  Sholkovitz  2006)  or  possibly  the  remineralization  of  organic  matter,  though 
organic  carbon  concentrations  are  low  within  the  freshwater  zone  (Bone  et  al.  2007).  At 
the  CCC  line  and  east,  a  reducing  zone  within  the  freshwater  plume  is  consistently 
observed  within  the  larger  freshwater  plume.  At  the  CCC  line  this  reducing  zone  contains 
large  concentrations  of  ammonium  (>50  pM)  (Charette  et  al.  2005;  Kroeger  and  Charette 
2008),  manganese  (up  to  15  pM)  (Charette  et  al.  2005)  and  ferrous  iron  (up  to  175  pM) 
(Charette  et  al.  2005).  In  areas  where  this  plume  occurs,  sediment  coring  reveals  a  zone 
of  iron-oxide  coated  aquifer  sands  (Charette  and  Sholkovitz  2002).  However,  previous 
studies  have  reported  that  nitrate  is  the  major  dissolved  nitrogen  species  entering  Cape 
Cod  coastal  waters  with  ammonium  accounting  for  only  a  small  percentage  (Valiela  et  al. 
2000;  Kroeger  et  al.  2006).  The  microbiology  and  geochemistry  of  three  major  parcels  of 
water  in  the  Waquoit  Bay  groundwater  aquifer  will  be  discussed  below. 

Surface  Salt  Zone 

The  surface  salt  zone  originates  from  tidal  or  wave  induced  overtopping  of  bay 
water  on  to  the  surface  sediments  and  mixing  with  the  advecting  freshwater  plume.  This 
water  is  recharged  at  high  tide  through  clean  coarse  quartz  sands  and  has  a  short 
residence  time  on  the  order  of  hours  or  days  (Michael  et  al.  2005).  Furthermore,  tracer 
injections  have  shown  that  this  saline  water  never  reaches  deeper  than  1.2  m  into  the 
sediment  column  (Michael  et  al.  2005).  The  shallow  salt  zone  contains  nitrate  mixed  in 
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from  the  freshwater  plume,  oxygen,  and  on  occasion  ammonium  resulting  from  saltwater 
intrusion.  Atmospheric  deposition  has  been  implicated  as  the  source  of  the  shallow  nitrate 
east  of  our  study  site  (Kroeger  and  Charette  2008),  however,  the  and  of  the 
nitrate  pool  at  our  study  site  is  not  consistent  with  an  atmospheric  source.  Alternatively, 
the  shallow  nitrate  could  be  the  result  of  mixing  the  freshwater  plume  up  into  the 
overtopping  seawater.  Plotting  nitrate  against  salinity  for  the  upper  meter  is  consistent 
with  mixing,  however,  such  mixing  does  not  explain  the  enrichment  in  both  N  and  O 
isotopes  over  this  range.  This  parcel  of  water  was  only  observed  in  the  two  most 
shoreward  piezometer  locations  (locations  A  and  B,  Figure  2)  in  our  study  area. 

Nitrogen  biogeochemistry  in  this  zone  is  complex.  The  ammonium  pool  is  small 
in  this  setting  and  is  only  detectable  in  a  few  samples  in  the  most  shoreward  site  and  even 
then  only  in  the  upper  half  meter  with  a  maximum  of  5  pM.  The  concentration  of  nitrate 
in  the  two  most  shoreward  piezometer  sites  shows  a  decrease  in  concentration  (ca.  40  to 
12  pM)  moving  from  1  m  up  to  the  surface.  With  this  decrease  in  nitrate  concentration, 
the  8^^Nno3  and  8^^0no3  increase  (4  to  20%o  and  2  to  10%c  respectively)  from  1  m  to  the 
surface.  This  concurrent  decrease  in  concentration  and  increase  in  both  isotope  values  is 
consistent  with  a  fractionating  consumption  process.  The  ratio  (0.56,  r^=  0.83)  of  8'*Ono3 
and  8^^Nno3  during  this  fractionating  process  is  consistent  with  previous  studies  of 
groundwater  and  sandy  aquifer  denitrification  .  Molecular  queries  show  a  microbial 
community  similar  to  those  found  in  marine  and  estuarine  sediments.  The  recovered  16S 
rRNA  gene  sequences  include  representatives  that  are  most  closely  related  to  Stappia- 
a.nd  Mesorhizobium-like.  organisms.  These  genera  encompass  a  large  variety  of  metabolic 
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lifestyles,  including  known  denitrifying  organisms  though  it  is  tenuous  to  assign 
ecological  function  based  on  16S  rRNA  phylogeny.  However,  nitrite  reductase  (nirK  and 
nirS)  gene  sequences  were  also  recovered  throughout  the  sediment  core.  In  the  upper 
meter  of  sediment  the  copper-based  nirK  gene  was  dominant  over  the  cytochrome-based 
nirS  gene. 

Despite  this  in  situ  chemical  and  molecular  evidence  that  denitrification  may  be 
important  in  the  upper  meter  of  the  aquifer,  incubation  studies  with  labeled  nitrate  did  not 
result  in  a  production  of  labeled  dinitrogen  (not  shown).  These  incubations  were  done 
approximating  in  situ  chemistry  at  the  time  of  sampling,  including  oxygen 
concentrations.  Elevated  oxygen  concentrations  (50-100  pM),  such  as  those  reported 
during  the  chemical  mapping  in  Chapter  3,  have  been  shown  to  limit  denitrification  in  the 
environment .  These  high  oxygen  concentrations  may  explain  the  observation  that  while 
denitrifiers  may  be  present  and  though  there  is  evidence  of  their  activity  in  the 
environment,  denitrification  was  not  measureable  in  our  bottle  incubation  over  the  time 
span  of  the  experiments.  On  the  other  hand,  we  did  find  substantial  evidence  of  active 
nitrification  in  the  surface  meter. 

Incubation  with  labeled  ammonium  resulted  in  the  transfer  of  the  label  to  the 
nitrate  pool.  This  transfer  of  nitrogen  from  the  reduced  to  oxidized  pool  is  mirrored  by 
the  up-regulation  of  ammonia  monooxygenase  subunit  A  (amoA).  Incubations  with 
labeled  ammonium  and  the  specific  inhibitor  allylthiourea  (ATU)  had  little  effect  on  the 
rate  of  isotope  transfer  in  these  incubations  or  the  expression  levels  of  the  archaeal  amoA 
gene.  The  presence  of  an  active  (defined  by  gene  expression  and  transfer  of  nitrogen) 
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archaeal  nitrifying  community  is  consistent  with  the  quantitative  study  of  the  arehaeal 
population,  which  demonstrated  a  peak  in  archaeal  amoA  genes  in  the  upper  half-meter  of 
sediments  (Chapter  2).  Recovered  amoA  genes  from  this  depth  are  most  closely  related  to 
sequences  obtained  from  water  column  and  low  salinity  estuarine  environments  including 
Monterey  Bay  and  northern  San  Franeisco  Bay. 

Taken  together,  this  zone  appears  to  have  both  denitrifying  and  nitrifying 
eapacity.  Given  the  short  residenee  time  of  the  water  in  this  area,  and  the  development  of 
fractionation  in  the  natural  abundance  isotopes,  denitrification  ought  happen  at  significant 
rates  to  maintain  the  isotope  signal.  Alternatively,  an  enriched  source  nitrate  mixing  with 
the  groundwater  could  produce  the  same  observed  patterns.  The  surface  bay  water  cannot 
be  the  souree  of  this  enriched  nitrate  as  that  pool  has  a  low  isotopic  signature  for  both 
atoms  (8^^N  =  6%c,  =  3%o).  Furthermore,  presenee  of  an  enriehed  nitrate  source  is 

less  likely  than  our  inability  to  maintain  the  denitrifier  community  in  bottles.  Perhaps 
microenvironments,  such  as  grain  surfaces  or  grain  boundaries,  where  oxygen 
concentrations  are  limited  by  the  balance  of  consumption  processes  like  aerobic 
heterotrophy  and  diffusion  from  the  surrounding  waters  are  important  in  the 
biogeoehemical  eycling  of  nitrogen  in  the  aquifer.  The  slurry  approaeh  used  in  the 
isotope  tracer  (^^NH4'^  addition.  Chapter  4)  would  necessarily  eliminate  such  zones.  In 
fact,  the  slurry  ineubations  were  shaken  so  that  the  diffusion  of  reactants  would  not  limit 
metabolism.  Of  course  in  the  case  of  an  inhibiting  substance  such  as  oxygen  to 
denitrifiers,  agitation  prevents  the  development  of  a  proteeted  harbor  for  the  organisms. 
Nitrification  is  not  limited  in  the  same  way  by  this  experimental  design.  However,  it  is 
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important  to  note  that  nitrification  in  these  surface  sediments  is  not  affected  by  the 
addition  of  the  traditional  bacterial  nitrification  inhibitor  ATU.  Taken  with  the  dominance 
of  the  archaeal  amoA  gene,  the  data  is  interpreted  as  the  upper  meter  of  the  aquifer 
harbors  an  active  population  of  archaeal  nitrifying  organisms  that  are  not  inhibited  by  the 
presence  of  substrate  level  ATU. 


Freshwater  Plume 

The  freshwater  plume  is  characterized  by  high  concentrations  of  nitrate  and 
oxygen,  and  to  the  east  of  the  CCC  line  high  concentrations  of  ammonium  in  a  reducing 
plume.  At  the  current  study  site  the  concentration  contours  of  nitrate  and  oxygen  in  the 
freshwater  plume  are  consistent  throughout  the  transect,  only  shoaling  as  the  transect 
progresses  seaward  and  the  deep  salt  wedge  shoals.  Contrary  to  enriched  8^^Nno3 
reported  near  the  CCC  line  (Nielsen  et  al.  1990),  both  8'^Nno3  and  8'^Ono3  remain  steady 
throughout  the  transect  (Chapter  3).  Addition  of  label  to  the  nitrate  or  ammonium  in  the 
isotope  pairing  experiments  does  not  result  in  a  production  of  labeled  dinitrogen.  At  this 
depth  the  denitrification  genes,  nirK  and  nirS,  show  similar  levels  of  abundance.  Similar 
to  the  surface  sediments  a  transfer  of  the  label  from  the  ammonium  pool  to  the  nitrate 
pool  is  observed  at  mid-depth,  with  potential  rates  similar  to  those  of  the  surface  meter. 
Interestingly,  the  addition  of  ATU  to  the  incubation  at  2  m  prevented  the  transfer  of  label 
between  the  reduced  and  oxidized  pools.  The  sensitivity  of  nitrification  to  ATU  is 
consistent  with  a  bacterially  mediated  process.  This  is  congruent  with  the  observation  that 
the  nitrate  peak  in  the  freshwater  plume  is  the  only  zone  where  the  bacterial  version  of 
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the  amoA  gene  is  dominant  over  the  archaeal  version.  Examining  the  expression  data  is 
consistent  with  an  active  bacterial  nitrifying  population  with  little  contribution  from  their 
archaeal  counterparts  as  seen  by  the  decline  in  both  rate  of  transfer  and  bacterial 
amoA  expression  in  the  presence  of  ATU. 

Deep  Salt  wedge 

The  deep  salt  wedge  is  characterized  by  a  sharp  transition  in  salinity  from 
freshwater  to  near  bay  water  salinity  in  a  half-meter.  This  study  and  previous  studies  have 
had  difficulty  obtaining  samples  with  enough  nitrate  for  natural  abundance  isotopic 
measurement  of  the  nitrate  pool  due  to  the  rapid  loss  of  nitrate  between  the  nitrate  rich 
freshwater  and  nitrate  lacking  salt  wedge.  Molecular  evidence  suggests  a  peak  in  archaeal 
amoA  genes  in  the  bottom  of  the  nitrate  plume  (2.8  m)  and  an  increase  in  nirS  genes 
going  into  the  deep  salt  wedge.  Also  of  interest  is  the  phylogenetic  relationship  between 
this  deep  archaeal  amoA  population  and  the  shallow  population.  While  the  shallow 
population  is  closely  related  to  water  column  type  sequences,  genes  recovered  from  the 
deeper  in  the  aquifer  are  more  closely  related  to  sediment  and  soil  sequences.  Like  the 
shallower  depths  the  paired  isotope  experiments  did  not  show  a  transfer  of  tracer  to  the 
dinitrogen  pool.  Unlike  the  shallower  sediments,  however,  the  labeled  ammonium 
experiments  did  not  show  a  transfer  of  label  in  the  nitrate  pool.  In  addition,  the  gene 
expression  data  does  not  show  significant  transcript  activity,  compared  to  the  shallow 
depths. 
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What  about  Anammox? 


Recent  years  have  seen  a  new  appreciation  of  the  potential  role  of  anammox,  the 
concurrent  reaction  of  ammonium  and  nitrite  to  form  dinitrogen  in  anaerobic  zones,  in 
many  environments.  If  present  at  the  site  of  our  study,  the  isotope  pairing  experiment 
ought  to  have  revealed  anammox  activity  through  the  production  of  N2  from  labeled 
^^N-ammonium  and  unlabeled  nitrate  or  nitrite  (as  well  as  similar  production  from 
unlabeled  ammonium  and  labeled  nitrate).  However,  our  experiments  did  not  show  an 
active  anammox  process  at  any  depth  tested.  Like  denitrification,  this  lack  of  chemical 
evidence  for  anammox  may  have  been  a  result  of  our  incubation  method  involving 
sediment  slurries.  Anammox  organisms  may  require  protected  microenvironments  in  the 
aquifer  that  were  not  available  in  the  agitated  slurry  incubations. 

Other  evidence  points  towards  the  presence  of  anammox-capable  microbes  in  the 
STE  sediments.  Planctomycetes-like  sequences  were  recovered  from  the  16S  rRNA  gene 
clone  libraries;  although  they  were  minor  a  contribution.  In  addition,  Planctomycetes-like 
sequences  were  not  closely  related  to  the  known  anammox  organisms,  and  therefore  the 
presence  of  these  particular  clones  is  not  a  strong  indication  of  the  potential  for  anammox 
processes.  However,  concurrent  study  of  the  recovered  sediment  core  has  revealed  the 
presence  of  ladderane  lipids  throughout  the  recovered  core  (Saenz  personal 
communication).  Ladderane  lipids  are  thought  to  be  specific  indicators  of  anammox 
organisms,  which  use  the  lipid  to  construct  the  anammoxosome  compartment  where  the 
anammox  reaction  occurs.  Interestingly,  a  second  lipid  crenarchaeol,  specific  for 
crenarchaeota,  was  shown  to  be  closely  associated  with  the  distribution  of  ladderanes. 
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which  may  indicate  a  relationship  or  shared  eeology  between  the  erenarehaeota  and 
anammox  organisms  (Saenz  personal  communieation).  Kroeger  and  Charette  (2008)  have 
also  reported  the  concurrent  loss  of  ammonia  and  nitrate  with  increasing  8^^Nno3  in  the 
deep  salt  wedge  as  well  as  the  freshwater  plume  (on  the  boundary  of  an  intruding  reduced 
ammonia  plume).  However  as  the  authors  note,  this  concurrent  loss  and  potential 
fractionation  could  also  be  the  result  of  eoupled  nitrification/denitrifieation,  or  reaction 
with  Mn-oxides  on  the  sand  grains. 

The  current  studies  show  an  abundant  archaeal  population  in  the  surfaee  meter 
and  are  likely  the  main  nitrifiers  in  the  shallow  sands,  however,  the  eeologieal  funetion 
the  deep  arehaea  population  is  uncertain.  Recent  work  in  the  Black  Sea  has  shown  a 
relationship  in  the  abundanee  of  arehaeal  amoA  genes  (marker  for  erenarehaeal 
ammonium  oxidizing  organisms)  and  anammox  speeific  nirS  gene  abundance,  and  linked 
the  two  proeesses  by  observing  nitrifieation  using  ^^N-label  in  the  suboxie  zone  (van 
Niftrik  et  al.  2004).  Is  it  possible  that  this  sort  of  relationship  oecurs  in  the  aquifer?  Are 
these  deep  arehaea  themselves  denitrifiers?  Further  work  must  be  done  in  this  area  before 
the  significance  of  anammox  processes  in  this  groundwater  aquifer  can  be  assessed.  This 
work  should  inelude  moleeular  data,  sueh  as  funetional  genes,  or  enriehment  eultures  that 
implieate  speeific  anammox  organisms  and  rule  out  mineral-pore  water  interaetions. 

Implications 

The  Waquoit  Bay  groundwater  aquifer  is  souree  of  nitrogen  to  the  eoastal  oeean. 
The  seale  of  nitrogen  introduced  to  through  this  avenue  is  a  function  of  the  nitrogen 
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loading  in  the  watershed,  the  flux  of  water  from  the  watershed  to  the  coastal  ocean  and 
any  transformation  processes  that  can  attenuate  the  dissolved  nitrogen  en  route.  Above 
and  in  the  previous  chapters  we  demonstrate  geochemical  signals  of  nitrogen  removal  in 
the  upper  sediments  associated  with  groundwater  denitrification  coinciding  with 
potential,  archaea-mediated  nitrification.  Thus  these  aquifer  sediments,  especially  the 
surface  sediments,  may  represent  an  area  where  anthropogenic  nitrogen  may  be 
ameliorated.  It  is  hypothesized  that  all  discharging  water  must  transit  through  such  a 
salinity  zone  before  entering  the  surface  waters  of  the  bay  and  if  these  zones  function 
similar  to  the  surface  salt  zone  significant  nitrogen  could  be  removed  from  the 
discharging  waters. 

Nitrogen  transiting  through  the  freshwater  plume  remains  relatively  unchanged 
over  the  course  of  the  transect,  due  to  the  relatively  slow  rate  of  reactions  compared  to 
the  advection  rate  of  water.  If  this  water  were  to  exit  into  the  bay  in  this  state  it  would 
represent  new  nitrogen  entering  the  bay.  However,  previous  work  has  estimated  that  65% 
of  the  nitrogen  flux  in  the  freshwater  plume  is  lost  over  a  distance  of  12  m  (14  day 
transit)  near  the  CCC  line,  much  of  this  was  reported  in  mixing  zones  between  water 
parcels  (Kroeger  and  Charette  2008).  At  our  vertical  mixing  between  parcels  appears 
limited,  as  evidenced  by  steep  gradient  especially  between  the  freshwater  plume  and  the 
deep  salt  wedge.  However,  it  is  reasonable,  and  predicted  that  much  of  the  water  flowing 
into  the  bay  will  transit  through  a  mixing  zone.  We  have  demonstrated  that  mixing  zones 
such  as  the  upper  meter  can  harbor  populations  that  can  transform  inorganic  nitrogen 
species  at  environmentally  significant  rates. 
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Future  Directions 


More  detailed  incubation  studies  are  needed  to  constrain  the  distribution  and  rate 
of  denitrification  in  the  aquifer  sediments.  The  next  critical  study  location  ought  to  be  the 
discharge  zone.  As  demonstrated  in  the  upper  meter  of  the  sediments  in  this  study,  the 
last  centimeters  or  meters  of  transit  prior  to  discharge  may  be  the  most  active  in  terms  of 
nitrogen  removal  from  groundwater  in  this  system.  However,  precious  little  is  known 
about  the  microbiology  and  geochemistry  of  this  zone  other  than  the  results  presented 
here.  Further  study  in  this  zone  will  yield  a  better  understanding  of  how  much  nitrogen 
can  be  attenuated  prior  to  discharge.  A  second  focus  for  future  studies  is  to  look  at  the 
geomicrobiology  in  other  areas  of  the  bay.  For  example,  how  does  denitrification  within 
each  water  parcel  east  of  the  CCC  line  compare  with  the  current  study?  Are  the  areas 
east  of  CCC  line  more  efficient  at  removing  nitrogen  due  to  the  presence  of  a  less 
oxidized  freshwater  plume?  Knowledge  of  this  type  may  influence  policy  makers  in 
prioritizing  such  issues  such  as  sewer  development  in  local  communities,  land  use  and 
water  rights.  Furthermore,  detailed  study  of  the  archaeal  ammonium-oxidizers,  especially 
those  presence  in  near  the  deep  salt  wedge  could  lead  to  a  better  understanding  of  the 
plasticity  of  their  metabolisms  and  their  overall  effect  on  nitrogen  biogeochemical 
cycling  in  a  variety  of  environments. 
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Figure  1.  (A)  Head  of  the  bay,  Waquoit  Bay,  Cape  Cod,  MA.  Contours  represent  10  ft. 
intervals  and  hatched  area  denotes  wetlands.  (B)  The  current  and  previous  study  location 
are  depicted.  Cape  Cod  Commission  wells  (CCCl,  CCC2),  transect  from  Charette  et  al. 
2005,  2006  (A  to  A’)  and  Kroeger  et  al.  (2008)  (A  to  A’,  B  to  B’  and  C  to  C’).  Mulligan 
and  Charette  (2006)  extended  the  B  to  B’  transect  east  and  west. 
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A  B  C  D 


Figure  2.  Schematic  of  SGD  at  Waquoit  Bay.  Three  major  water  parcels  are  evident,  a 
surface  salt  zone,  a  freshwater  plume  and  a  deep  salt  wedge.  Locations  A,  B,  C  and  D 
represent  area  of  piezometer  studies.  Sediment  cores  were  obtained  from  location  B. 
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